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SUMMARY 


Phenomena associated with long range propagation of sound over irregular topog- 
raphy motivated the research work which was described in this thesis. Specifically, 
the goal of the work was to analyze the diffraction effects which would occur near 
the tops of hills and ridges. From this particular goal, the research work evolved 
into a study of the diffraction of a high frequency plane wave due to its grazing of a 
two-dimensional curved surface of finite impedance. Laboratory scale models were 
constructed and measurements were made of the field on, above, and behind either of 
two curved surfaces possessing distinctly different impedances; that is, one was soft 
while the other was hard. The experimental technique consisted of simultaneously 
measuring the pressure at a reference point and at a field point due to a transient 
pulse generated by an electric spark. The pressure waveforms were digitized and 
processed. As described in the thesis, the ratio of the discrete Fourier transforms 
of the two wavefor ms provided an estimate of the insertion loss between them. The 
results of the measurements were compared with the predictions of a theory which 
was derived by Pierce using the method of Matched Asymptotic Expansions (MAE). 
The predictions relied upon the experimental evaluation of the impedance of each 
surface at grazing angles of incidence. This evaluation was achieved by a fairly stan- 
dard technique involving empirical models of various generic types of surfaces. An 
example was shown of the important role that the structural intricacies of a surface 
play in the determination of an appropriate model. The comparison between the 
measurements and predictions clearly indicated that the theory gives an excellent 
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description of the field anywhere near a curved surface. Further, with a simple 
modification, the theory was also shown to give nearly as good of a description of 
the field surrounding a curved surface even at distances far behind the surface yet 
near the line of sight. 
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CHAPTER I 

INTRODUCTION 

The principal goal of the research discussed in this thesis was to gather and analyze 
experimental data that would be useful in understanding the effects of irregular to- 
pography on long range sound propagation in the atmosphere. At the outset, it was 
assumed that the most interesting and influential effects on such propagation occur 
near the tops of hills or ridges located in the terrain. In the vicinity of these tops, 
rays associated with the propagating sound wave are incident at angles near grazing, 
particularly when the source of the wave is distant enough that the incident field is a 
nearly vertical plane wave. A real example of such circumstances occurs when a low 
flying aircraft produces an acoustic disturbance which propagates over distant hilly 
terrain. It can be shown that an analogous situation involving flat ground occurs 
when wind or temperature gradients exist in the atmosphere . 1 Curved ray paths 
and possibly multiple shadow zones are present in both fields. In either situation, 
the interaction between the impinging sound field and the surface topography is 
expected to be complex. In light of this expectation, a prototype laboratory scale 
model of a single, cylindrical ridge was constructed of plywood. The expectation 
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was that a knowledge of the basic effects of a single ridge on the incident field could 
be used to construct the field over more complex topography. Therefore, experi- 
ments were designed to measure the insertion loss on, above, and behind the single 
scale model ridge. A spark generator provided an acoustical point source which 
was located many wavelengths away from the ridge. The frequency content of a 
generated waveform was such that kR » 1, where R was the radius of curvature 
of the scale model ridge. The large value of kR was in keeping with the typically 
large scales associated with hills or ridges in real terrain. 

In a related vein, the duration of a typical waveform was approximately 40 fis. 
This short duration made it possible to identify the desired direct, reflected, and 
diffracted waveforms separate from other spurious waveforms (reflected from walls 
and objects in the lab room) . The insertion loss was calculated as the ratio of the 
discrete Fourier transform of a field waveform to that of a reference waveform, after 
appropriate corrections for spreading and absorption losses. In this context, a field 
waveform typically referred to a measured waveform which had been reflected or 
diffracted from a ridge, while a reference waveform referred to a waveform measured 
prior to any interaction with a ridge. 

An important parameter in such propagation problems is the specific acoustic 
impedance of the surface of the topography. Two types of surfaces were used in 
the scale model experiments; they were a plywood curved surface and a carpet- 
on-plywood curved surface. It was found that these two surfaces could reasonably 
represent a hard and soft ground, respectively. An important part of the experi- 
ments was the quantitative evaluation of the specific acoustic impedance of the two 
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surfaces used on the model ridge. The evaluation was made by assuming that the 
impedance of the surface when curved was equivalent to that when the surface was 
flat. Then, measurements were made of the field over flat table tops which had 

c 

surfaces nearly identical to those of the model ridge. The theory for sound prop- 
agation over a flat impedance plane was employed along with empirical models of 
hard and soft surfaces to deduce the impedance from measurements of the pressure 
at locations just above the surface. In particular, a model developed by Delany 
and Bazley for the impedance of fibrous absorbent materials and an asymptotic 
version of a five parameter model developed by Attenborough were used to describe 
the caxpet-on-plywood and plywood surfaces, respectively, in terms of an effective 
flow resistance parameter, o . The results indicated that values in the range of 
0.8-1.6 x 10 6 kg/(m 3 -s) adaquately described the carpet-on-plywood surface while 
a value of 120 x 10® kg/(m 3 -s) sufficed for the plywood surface. 

Subsequent to the measurements over the flat table top surfaces, measurements 
of the field on, above, and behind the model curved surfaces were made. Specifically, 
the pressure was measured on the curved surfaces as well as along vertical lines 
located 0, 63, 120, 200, and 400 cm downstream from the apices of the model curved 
surfaces. These distances were chosen knowing that the radius of curvature, R, of 
the model surfaces was approximately 2.5 m. The spark source was located roughly 
R meters in front of the model surfaces. The generated acoustic pulse was typically 
40 fis in duration with a peak amplitude of approximately 130 Pa at 1 meter. A 
discrete Fourier transform of such a captured pulse revealed a broad frequency 


band such that kR » 1. 
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The results of the experiments were compared to the computational results 
of the asymptotic theory for diffraction over curved surfaces of finite impedance. 
The theory was based upon the work of V.A. Fock, who investigated the diffracted 
field of an electromagnetic wave propagating over a curved body. A.D. Pierce re- 
developed the theory in an acoustic context, and employed the method of Matched 
Asymptotic Expansions to produce a solution which matched the impedance bound- 
ary condition at the surface and which asymptotically approached the geometrical 
(ray) acoustics solution at locations not too far above the surface. Highlights of 
the solution include an application of the method of steepest descent to produce a 
approximate solution in terms of the field behind a rigid knife-edge plus two slowly 
varying so-called background terms. The background terms were shown to vanish at 
large distances downstream from the apex, thus leaving the functional form of the 
solution equivalent to that of the solution for the field far bahind a rigid knife-edge. 
However, it was shown that, at such distances, a dimensionless height involved in 
the solution had to be replaced by another more appealing dimensionless height in 
order for the solution to yield a good comparison with experimental data. This 
result did not invalidate the theory because the original solution was derived under 
the assumption that the listener was close to the ridge. In fact, the theory com- 
pared excellently with experimental data at distances from the apex which were 
less than R . Further, an application of Cauchy’s Residue theorem resulted in an 
approximate form of the solution in terms of an infinite residue series often referred 
to as a creeping wave series. Though thought to be valid only in the shadow zone 
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of a curved surface, the series solution was shown to give excellent results well into 
the lower bright zone. 

This thesis was assembled so that the reader would be first introduced to some 
history and then to the theory. Subsequently, the laboratory, and the preliminary 
and main experiments, and their results, were described in detail. The main body 
of the thesis was concluded by a short summary section of conclusions and recom- 
mendations. Several appendices, which contain certain details of the work, were 
included as well. 
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CHAPTER II 

A GLIMPSE AT THE PAST 


This thesis contains the results of a study of the propagation of sound over a 
curved surface of finite impedance. Of particular interest were the diffraction effects 
inherent in such a propagation. It so happens that in the area of wave propagation 
over curved surfaces, a majority of the work done previously has dealt with electro- 
magnetic waves. For this reason, much of the discussion of previous work, whether 
theoretical or experimental, involved electromagnetic waves. However, the similar 
wave mechanics involved in either propagation made all of the discussed work useful 
and pertinent. In fact, the foundation of the current theory of interest was the re- 
sult of a careful translation between electromagnetic and acoustic variables from a 
theory derived for long range radio wave propagation over the surface of the earth. 

2.1 ELEMENTS AND EVOLUTION OF THE THEORY 

When the literature search was being conducted, it was not unusual to find refer- 
ences to some fairly old work by men of considerable standing such as Sommerfeld, 
Kirchhoff, Fresnel and even Newton, amongst others. Much of this work pertained 
to light and its propagation. In this way, a certain, although surely limited, view of 
the history of thought on diffraction, and wave propagation in general, was obtained. 
This history is intrinsically interesting. It also highlights some of the key concepts 
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involved and it lends continuity to the current and more recent work. Therefore, 
an attempt is made in this chapter to sketch the evolution of thought on the nature 
of light and its ability to diffract. The sketch demonstrates some of the interlaced 
traditions of light and sound propagation. 

The current theory of interest is an example of the interlaced tradition in 
that the concepts involved were developed in the area of radio wave propagation. 
Furthermore, and in particular, the presented theory results from an asymptotic 
solution of the wave equation; but the result is explained in terms of ray concepts. 
In high frequency cases, wave concepts and ray concepts are quite compatible. 
However, this compatibility was not always well understood or accepted, even by 
Newton hi ms elf. There is evidence that Newton noted the ray-like behavior of 
light in co mm on surroundings as a part (but far from all) of his argument for the 
corpuscular (particle) nature of light. The ensuing struggle to resolve the wave- 
particle dispute led to an understanding of ray and wave mechanics which, as is 
well known, also have applications in the realm of acoustics. 

2.1.1 Some Early Thoughts and Observations of Diffraction 

An early, and perhaps the first, “careful” observation of the optical phenomenon 
known as diffraction was made by Francesco Maria Grimaldi (c. 1650). 2 | He noted 
that light does not always travel in straight lines but can bend around comers. This 
observation was subsequently made by Newton, Hooke, and Huygens among others 
before the turn of the 18th century. Although all concerned parties agreed that light 
does actually bend around corners or edges, they did not agree on the explanation. 
The scientific community at the time was split between the corpuscular theory 

f The year(a) of a discovery, publication or other event is usually enclosed in parentheses 
and placed nearby the reference to the event. Hopefully, this approach has enhanced 
the readability of the text, and will lend the reader a sense of the chronology without 
the strain of searching through the bibliography. 
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of light which was championed by Newton and the wave theory of light which 
was promoted by Huygens. Each theory put forth explanations for interference, 
refraction and diffraction effects. Newton argued that if light propagated as a wave, 
like sound does, then “it would bend into the Shadow”. He noted that “Rays which 
pass very near to the edges of any Body, are bent a little by the action of the 
Body,” but, he continued, that “So soon as the Ray is past the Body, it goes right 
on .” 3 From a practical point of view, this observation was no doubt correct for 
the range of bodies he observed. However, Newton was misled by the difference in 
the quantitative size of the diffraction effect between light and sound when each 
interacts with ordinary objects. This misunderstanding is clearly present in the 
following remark which he used as part of his argument against the wave theory of 
light. 


For Pression or Motion can-not be propagated in a Fluid 
in right lines, beyond an Obstacle which stops part of 
the Motion, but will bend and spread every way into the 
quiescent Medium which lies beyond the Obstacle.^... 
The Waves on the Surface of stagnating Water, passing 
by the sides of a broad Obstacle which stops part of 
them, bend afterwards and dilate themselves gradually 
into the quiet Water behind the Obstacle. The Waves, 
Pulses or Vibrations of the Air, wherein Sounds consist, 
bend manifestly, though not so much as the Waves of 
Water. For a Bell or a Cannon may be heard beyond a 
Hill which intercepts the sight of the sounding Body, and 
Sounds are propagated as readily through crooked Pipes 
as through streight ones. But Light is never known to 
follow crooked Passages nor to bend into the Shadow . 4 


Apparently, the notion of mechanical waves in air and water was well accepted 
but the notion of light waves was in contest. Obviously, accepted analogies between 
light and sound were quite limited at the time. Further, Newton was able to observe 
several phenomena such as the heating of bodies by light which he contended could 
be explained only from a corpuscular point of view. “If Light consisted only in Pres- 
sion propagated without actual Motion, it would not be able to agitate and heat the 
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Bodies which refract and reflect it.” Newton’s Opticks Dover ed. Impressed by these 
observations, Newton held to his corpuscular theory of light, and his prodigious 
reputation ensured this theory a large and faithful following. 

A major step toward the resolution of this split occurred in 1802 when Thomas 
Young announced his principle of interference. Nearly two decades later, but un- 
aware of Young’s announcement, Augustin Jean Fresnel announced in 1818 his own 
principle of interference. Fresnel was able to explain the relationship between the 
wavelength and the size of diffraction and interference effects. Fresnel relied heavily 
upon Huygens’ principle and explained “diffraction as the mutual interference of 
secondary waves emitted by those portions of the original wavefront which have not 
been obstructed by the diffracting obstacle.” 5 Shortly after Fresnel’s announce- 
ment, Hamilton began work in 1824 on the formulation of the geometrical theory 
of optics (or “ray theory”). His work was based upon Fermat’s Principle of Least 
Time; that is, that the optical distance between any two points is stationary. As 
a result of this and other work, the relationship between wave propagation and 
ray optics became fairly well understood by the time of the publication in 1865 of 
Maxwell’s Theory of Electromagnetism. Maxwell’s theory crowned the wave theory 
of light and, in so doing, made clear the analogy between light and sound. 

Further investigations into the relationship between ray optics and wave propa- 
gation were made by Kirchhoff. In 1882, he published a “proof” which showed that, 
as the wavelength approaches zero, the results given by the wave theory approach 
those given by geometrical optics. The strength of this proof was mitigated by some 
shortcomings in the Kirchhoff-Huygens principle upon which it was based. Perhaps 
a better and more widely accepted account of the relationship between geometri- 
cal optics and wave theory was given by Sommerfeld and Runge (1911). Later, 
Luneburg (1944) and Kline (1954) showed that an asymptotic expansion in inverse 
powers of the wave number can be found for the solution of the wave equation in 
the illuminated region of a given scatterer. The leading term of this expansion is 
found to be the geometrical optics solution. 6 Thus, the relationship between the 
concept of waves and that of rays came to be understood in terms of ray theory 
being a high frequency approximation of wave theory. 
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2.1.2 Long Range Radio Wave Propagation and the Shadow Zone 

Whereas the diffraction of light around common objects is limited to very small 
scales and effects, the propagation of radio waves around the surface of the earth 
involves scales on which these effects are more readily observed. In the early 1900 s, 
a great deal of interest arose in phenomena associated with the transmission of 
radio waves beyond the line of sight and into the geometrical shadow (generica Uy 
termed the shadow zone ) of the earth. Generally, the model used to describe this 
situation was that of a vertical electric dipole in the presence of a metallic sphere. 
The radius of the sphere was assumed to be much greater than any wavelength of 
interest. Unfortunately, under such an assumption the eigenfunction expansion of 
the solution in the shadow zone is slowly convergent. However, Watson (1918), who 
drew upon the work of Poincare (1910) and Nicholson (1910), managed to express 
the expansion in terms of a residue series. From this residue series, an associated 
contour integral was constructed. Watson showed by deforming the contour in a 
particular way that a new, faster converging series can be derived. This process is 
known as Watson’s Transformation. White (1922) extended the work of Watson to 
include cases when the observer is outside the shadow zone. In the late 1930’s, van 
der Pol and Bremmer (1936-38) published several notable papers on the propagation 
of radio waves over the surface of the earth. (With the same theory, they also 
addressed the optical phenomenon known as a rainbow.) Their work included the 
use of what they termed “the third order approximation” in the saddle point method 
for the solution of diffraction integrals of the type dz. This approach yields 

a solution involving Airy functions which would be central in the work of V.A. Fock 


soon to follow. 
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2.1.3 The Local Field in the Penumbra 

Shortly after the work of van der Pol and Bremmer, in 1945, V.A. Fock (Fok) 
published a paper in which he put forth the “Principle of the Local Field in the 
Penumbra Region.” Fock considered the case of a high frequency wave incident 
upon a perfectly conducting convex surface of continuously varying curvature. By 
way of a physical argument, Fock claimed that the surface current distribution in 
and near the shadow boundary depends upon the local curvature and upon the 
incident field. This dependence was expressed in the form of a “universal” function 



"^rr da 

w;(a) 


( 2 . 1 ) 


where 1) £ is a representative distance of a given point on the surface from the 
geometric shadow boundary, and 2) T is the contour which runs from infinity to the 
origin along the path arg(z) = 2 tt/ 3 and from the origin to infinity along the real 
axis. u>; (a) is simply related to the derivative of the Airy function and is a bounded 
solution to w"(a) = aw(a). The function G(£) decreases exponentially with £ for 
large positive values of £ in the shadow zone and asymptotically approaches 2 for 
large negative values of £ in the illuminated region. The characteristic width of the 
penumbra was given as d = [AR 2 /*] 1 ' 3 , where Rc is the radius of curvature of the 
surface at the grazing point of the incident wave and in the plane of incidence. 

In the following year, 1946, Fock derived the same result directly from the 
Helmholtz equation. In addition, Fock extended the solution to the region in the 
“neighborhood” of the surface and included bodies which are “good” conductors 
in the sense of the Leontovich boundary condition. 7 (In acoustic terms, the Leon- 
tovich boundary condition (used below) is analogous to the impedance boundary 
condition associated with locally reacting surfaces.) In short, Fock began with the 


Helmholtz equation 


V 2 tf + k 2 V = 0 


( 2 . 2 ) 
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but then wrote 

^ = c <fc *^* (2.3) 

and assumed 1) that 'i'* varies slowly over a distance of one wavelength, and 2) 
that variations in '£* normal to the surface are much greater than variations in \t* 
tangential to the surface. Then, he transformed the local coordinates to a parabolic 
system (£,rj) where r) = 0 represents the local surface. After further manipulation, 
the solution was shown to be 

= g - * 4 '' e* 4 */ 3 V (2-4) 


where V is the solution to 


with the boundary condition 


d 2 V dV , I/ n 

__ + ,_ + „v = o 


+ qV = 0 at rj = 0 . 

on 


(2.5) 


( 2 . 6 ) 


Fock used the Sommerfeld radiation condition as the exterior boundary condition. 
The parameter q is a dimensionless surface admittance or conductance. One iden- 
tifies a diffusion-like operator d 2 / dr] 2 . According to Fock, Malyughinetz first in- 
troduced the notion of transverse diffusion which arises due to this operator. 8 As 
its name implies, the transverse diffusion is in the t] direction; that is, at any point 
the diffusion is along lines of constant (. Therefore, it contributes to the transport 
of energy into the shadow zone. According to Fock, the transverse diffusion is in- 
significant except near the line of sight and near the surface. (Interestingly, (2.5) is 
of the same form as the SchrSdinger equation.)! A solution for V was found in the 
form of a contour integral 

4= / U («-,)- 1 i B j ~ (« - «)] *"■ « da (2.7) 

l,q) 2^ J c i V ’ rv[ (a) - qw x (a) 1 

where (o) is a second, linearly independent Airy function. 9 At the surface (rj — 
0), and for a perfectly conducting surface ( q = 0), the function V reduces to (2.1) 
derived by Fock the previous year. 

t A decimal number contained in parantheeee it a reference to a displayed equation of 
the tame number. Numbered flguret, tablet, etc. are epecifleally referred to at tuch. 
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2.1.4 Geometrical Theory of Diffraction 

Approximately one decade after the publication of Fock’s work, Keller (1956) pub- 
lished the “Geometrical Theory of Diffraction.” Keller’s theory was based upon a 
heuristic application of Fermat’s principle. In order to account for diffraction effects, 
Keller added a modified form of the principle which states that “the diffracted rays 
connecting the points P and Q are those curves which render stationary the Fermat 
integral among all curves connecting P and Q, and having a point or arc on the 
diffracting obstacle.” 10 (Note: The original wording was "... arc on the cylinder.” 
since this was the geometry originally examined.) This statement can be aplied to 
any body of arbitrary shape. The resulting solution is in terms of unknown diffrac- 
tion coefficients”, which are analogous to the reflection coefficient, and which are 
found by comparison to well known exact solutions of various canonical problems. 
These coefficients are functions of the incident field, the wavelength, and the surface 
properties and geometry at the point of diffraction. These coefficients follow the 
spirit of Fock’s notion of the local field in the penumbra. Keller’s theory formalized 
the concept of “creeping” waves, so-called by Franz and Depperman (1952), for an 
arbitrarily shaped body. The “creeping” wave hypothesis is that, from the point of 
diffraction, a portion of the incident wave proceeds to travel into the shadow zone 
along a ray path coincident with a geodesic of the surface. This wave continuously 
sheds rays tangential to the geodesic ray path in accord with the modified Fermat s 
principle of diffraction. Keller’s theory is a high frequency approximation although 
good results have been reported at wavelengths on the order of the dimensions of 
the diffracting obstacle. 11 This theory is not valid near the geometric shadow 
boundary or on the surface of the diffracting obstacle; that is, the ray approach is 


invalid at caustics. 
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2.1.5 Recent Developments 

In 1959, Goodrich examined Fock’s theory in light of the work of Keller, Franz and 
Depperman, and Kazarinoff and Ritt (1959). Goodrich showed that by changing 
Fock’s choice of coordinate transformation that Fock’s theory will yield properly the 
creeping wave behavior discussed above. Also, Goodrich noted that Fock’s theory 
is essentially a 2-dimensional one in that it does not account for ray convergence on 
the surface of a 3-dimensional body. This difficulty was overcome in the derivation 
of Hong (1967). 

An analysis of previously published diffraction theories was made by Logan 
(1959). Subsequently, Logan and Yee (1962) presented a thorough mathematical 
treatment of the interrelationships between the relevant theories for diffraction from 
a convex body. More recently, Ivanov (1971) examined plane wave diffraction from 
an ideally reflecting cylinder. He proceeded to match an asymptotic solution valid in 
the shadow zone to an asymptotic solution valid in the illuminated region. Ivanov’s 
matched solution was continuous but not necessarily smooth at the shadow bound- 
ary. Later, Pathak (1979) found his own uniform, asymptotic solution. His solution 
is for a perfectly conducting surface and is valid everywhere except at the surface. 
Pathak’s solution is given in terms of well tabulated functions similar to those of 
Fock. 

The current theory of interest is that presented by Pierce (1986). His solu- 
tion employs the method of Matched Asymptotic Expansions to derive a uniform 
solution that 1) is valid on much of the surface, and 2) reduces asymptotically to 
the geometrical acoustics solution in the illuminated region. The theory also ac- 
counts for locally reacting surfaces of finite impedance. A summary of this theory 
is presented in chapter 3. 
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2.2 RECENT EXPERIMENTAL STUDIES 

A survey of related experimental studies is presented now. This survey includes the 
results of experimental studies of the diffraction of waves incident upon cones and 
other bodies of revolution as well as experiments on wave propagation over and on 
curved surfaces. As stated above, a significant portion of the relevant experimental 
work has involved electromagnetic wave propagation. An effort is made to include 
this work although an emphasis is placed upon work performed in the realm of 
acoustics. To the best of the author’s knowledge, the survey includes most of 
the relevant experimental work from the past two decades. The objective here is 
to show the extent of the experimental evidence which supports (or opposes) the 
theory discussed above. As is discussed below, the majority of experimental results 
have been found to be in good or excellent agreement with the theory. Typically, 
the greatest difficulties are encountered in regions near the line of sight. Further, 
it appears that there is a large amount of experimental data on scattering from 
screens and wedges compared to that on scattering from smoothly curved bodies. 
This discrepancy is due to interest in the design and performance of barriers. In this 
area, a legitimate question concerns the impact of having a curved barrier top versus 
a pointed or flat barrier top. Since this question will be addressed later in the thesis, 
a brief discussion of results involving screens and wedges is included. However, most 
of the discussion is intentionally tilted towards cases involving curved surfaces. 


2.2.1 Edgy Barriers: Screens and Wedges 

A renewed interest in experimental and theoretical work related to the performance 
of sound barriers occurred during the early 1960’s due to the public alarm over the 
rise in environmental noise pollution. This rise in noise pollution was attributable 
in large part to the proliferation of airplanes and automobiles, and the location of 
factories and facilities in urban and suburban centers. In 1968, Maekawa gathered 
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experimental data on the attenuation of sound over a wide barrier (screen). He was 
concerned principally with sound from sources near the barrier. Maekawa used this 
data along with the Kirchhoff-Fresnel solution for sound propagation over a semi- 
infinite screen to develop a chart method for predicting the attenuation of sound 
over a finite screen. Maekawa’s work included direct and ground reflected rays 
although the ground was idealized as rigid. Scholes, Salvidge and Sargent (1970) 
carried out a series of full experiments which measured the barrier performance 
of a screen for multiple receiver positions near and far from the screen. They also 
considered the effects of favorable and unfavorable wind conditions. In most of their 
experiments, the source was relatively close to the screen. The results for zero wind 
were in good agreement with those predicted by Maekawa’s method. 

Jonasson (1971,1972) and Ambaud and Bergassoli (1972) performed experi- 
ments on noise reduction from wedged barriers. Jonasson considered the effect of 
the surface impedance of the ground and wedge on barrier performance. The re- 
sults of his analytical model compare well with experiment for situations when the 
source and receiver are deep in the shadow zone and the closed angle of the wedge is 
greater than jt/ 2. Kurze (1973) published a literature survey in which the effects of 
various components of barrier performance were outlined. He analyzed the accuracy 
of the then current prediction methods with respect to various source and receiver 
positions, ground and barrier surface impedances, and wind conditions. Kurze ex- 
amined both wedge and screen type barriers. More recent experiments have been 
carried out by Bremhorst and Medwin (1978), and Thomasson (1978). Thomasson 
improved upon the work of Jonasson (1972) and Lindblad (1970) by using an in- 
tegral equation approach involving Green’s functions. In addition to the improved 
mathematical rigor of the solution, the results for finite screens of low height were 
good when compared with experiment. He also introduced an approximate method 
for “solution of the impedance boundary condition.” Bremhorst and Medwin ex- 
amined the diffraction of a point source pulse by an infinite screen and wedge. They 
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used the Biot-Tolstoy theory (1957) as a building block for the prediction of finite 
screen and wedge diffraction in the shadow zone. 

A rather voluminous study of highway noise propagation over barriers was 
published by Hayek et al. (1978). Their theoretical analysis included barrier ge- 
ometries such as screens, wedges, cylinder-topped wedges, and trapezoids. They 
also examined the influence of a finite impedance ground and barrier surface. The 
results indicated that some of the then current prediction procedures were insen- 
sitive to parameters which their analysis showed to be important in determining 
barrier performance. Some of these parameters were screen thickness, barrier sur- 
face im pedance, and interior wedge angle. The authors concluded their study by 
recommending further studies, particularly experimental and involving curved bar- 
riers. More recent results were presented by Nicolas et al. (1982) and deJong et 
al. (1982). deJong et al. measured the sound transmission loss for cases of finite 
and discontinuous impedances along the surface of the barrier and ground. After 
extending Jonasson’s method, they found good agreement between the calculated 
and measured results for the above cases. Nicolas et al. measured the diffraction of 
direct and reflected rays over a semi-infinite screen. Finite ground impedance effects 
were included. In cases where the source or receiver is near the barrier, they found 
a large discrepancy (sometimes on the order of 12 dB) between their experimental 
results and those predicted by the first-order Macdonald solution (Bowman et al.) 
as well as those predicted by the classic Kirchhoff- Fresnel theory. Much better 
agreement was found when their experimental results were compared with results 
based upon the method of Pierce and Hadden (1981). This method is derived from 
a reformulation of the problem of a point source in the vicinity of a rigid wedge. The 
method reduces quickly to various limiting cases of the source and listener positions 
and the closed wedge angle. The method was shown to be valid for arbitrary source 
and listener positions. An efficient numerical solution is found through the use of 
Laguerre integration. A comparison between the numerical results of the method 
and the experimental results of Ambaud and Bergassoli shows good agreement. 
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2.2.2 Bodies Born of Revolution 

Cones and other bodies of revolution offer simple 3-dimensional geometries 
which 1) are easy to find or construct and, thus, allow for experimental testing of 
theories such as the Geometrical Theory of Diffraction (GTD), and 2) readily admit 
numerical solutions due to symmetry considerations. For this reason, many recent 
papers have reported the use of cones, spheres, and ellipsoids in their model ex- 
periments. Keller (1961) made a comparison between theoretical and experimental 
results for the backscattering cross-section produced by the electromagnetic irradi- 
ation of a perfectly conducting cone. The experimental results were obtained by 
Keys and Primich (1959) for six different cones with half-angles ranging from 4 to 20 
degrees. The data is in good agreement with predictions based upon Keller’s Geo- 
metrical Theory of Diffraction. Bechtel (1965) carried out a series of experiments to 
measure the radar cross-section (RCS) of cones with half-angles of 4, 15, and 90 de- 
grees. The size of the base, a, was also varied between 0.98 and 2.87 wavelengths in 
diameter. Again, Bechtel compared his results with predictions based upon Keller s 
GTD. For ka on the order of 8 or 9, Bechtel found good agreement between the 
theory and experiment except for cases when the cone is observed within 30-40 
degrees of nose-on with a vertically polarized wave (i.e. wavefront parallel to the 
axis of the cone). For smaller cone bases on the order of ka = 3, the predicted 
RCS matched the observed RCS to within 5 dB but the predicted shape or form of 
the RCS was not well observed, particularly within 30 degrees of nose-on. Another 
set of reported data on cone diffraction was gathered by Bargeliotes et al. (1975). 
Mittra and Safavi-naini (1978) compared these experimental results with theory 
which considered the potential field as being produced by an array of point sources 
on the surface of a *cone. The array represented the surface currents as derived 
asymptotically by Fock and others. The comparison showed good agreement when 
the diffracted rays had traveled a significant angular distance around the cone. No 
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results for near grazing angles of incidence were found in their report. An extensive 
study of acoustic scattering from various bodies was made by Lang (1980). He 
measured forward and back scattering as well as creeping wave type phenomena 
from smooth, rigid objects such as baffles, cylinders, cones, and cone variants. 

Experimental data on high frequency diffraction and/or reflection from other 
bodies of revolution has been reported by a number of sources. Neubauer (1967) 
measured the travel time from source to receiver of a sonic pulse diffracted around 
a circular aluminum cylinder. The measurements were made in water. He posi- 
tioned a narrow-beam source such that the central ray arrived tangential to the 
cylinder surface and normal to the axis of the cylinder. An array of five equally 
spaced hydrophones measured the amplitude and arrival time of the pulse. The 
hydrophones were deep in the shadow zone. Neubauer found that his results were 
in excellent agreement with those predicted by the “creeping wave” theory of Franz 
and Keller. Neubauer also reported that a slit of less than a wavelength in width 
and positioned along the grazing line allowed a . .large part of the wave to pass. 
Further, he placed a baffle against the cylinder surface to block passage through 
the water and, subsequently, found little or no transmission. Thus, he concluded 
that the wave was transmitted on the water side of the cylinder boundary. Foxwell 
(1970), and Blake and Wilson (1977) carried out related experiments. Foxwell mea- 
sured the diffracted field on the shadowed surface of a rigid sphere. A point source 
located on the surface formed the illuminated pole. At high wavenumbers, Foxwell 
observed interference effects predicted by creeping wave theory. Blake and Wilson 
performed an analogous experiment using a highly eccentric prolate spheroid. Their 
sound source was located at a distance from the surface equal to the length of the 
spheroid and along the major axis. Measurements revealed an illuminated spot at 
the antipole and the existence of the creeping wave interference pattern near the 
antipole. In addition, Blake and Wilson reported that 1) the measured levels at 


20 


the antipole “decrease roughly as the reciprocal of frequency” , 2) that the shadow 
zone measurements are “well approximated by Keller’s Geometric Theory of Dif- 
fraction” , and 3) that deviations on the order of 3 dB from the theory “are observed 
at high frequencies and at coordinates off the antipole.”. Lang (1980) also examined 
acoustic scattering from the surface of a thin prolate spheroid and used the results 
to deduce characteristics of acoustic scattering from any smoothly curved surface. 
Lang concluded that the “backscattered pressure from smoothly curved bodies is 
determined almost exclusively by specular effects, even at wavelengths that are rel- 
atively large compared to the appropriate dimensions of the scatterers.” He also 
observed “evidence of surface fields very close to smoothly curved scatterers” that 
he “identified with Franz-type creeping waves.” 

More recently, Almgren (1986) measured the insertion loss above and behind 
convex and concave cylinders. A distant spark source provided a near planar inci- 
dent wave. Almgren was interested in examining the analogous relationship between 
sound propagation over a curved surface and that in a medium with a linear ver- 
tical sound speed gradient. He found that “it is reasonable to simulate the effect 
of refraction due to a sound speed gradient . . . with a curved ground scale model. 
Only small errors on the order of 3 dB or less were found for all measurements 
including those at grazing incidence. These errors were relative to the theory as 
presented in the works of Pekeris (1946), Pridmore and Brown (1962), and Pierce 
(1981). (Almgren also examined the impact of the acoustic boundary layer in cases 
involving a very hard surface. His data shows that the boundary layer makes only 
a small contribution to the total insertion loss; however, he did conclude that the 
boundary layer effect is important when comparing scale model results to those at 
full scale.) Berry and Daigle (1987) also were interested in the refraction of sound 
in a stratified atmosphere. Like Almgren, they chose to examine the analogous case 
of diffraction over a convex curved surface. A tone burst mechanism was used to 
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produce a narrow band point source. Their measurements were confined to listener 
positions deep in the shadow zone and in the penumbra region while the point 
source was placed either on or slightly above the surface of the ridge. Measure- 
ments were made over a frequency range of 0.3 to 10 kHz. For listener positions 
below the shadow boundary, they compared their results to predictions based upon 
an extended version of the creeping wave theory. For listener positions above the 
shadow boundary, they compared their results to predictions based upon geomet- 
rical theory. At the shadow boundary, they reported that the two theories agreed 
to within 1/2 dB. However, the experimental results differed from the theoretical 
results by as much as 5 dB in cases where the listener was either on or near the 
shadow boundary. Comparisons outside of this region were good. 

2.3 In a Nutshell 

Based upon the surveyed literature, it seems clear that although the qualitative 
nature of the sound field produced by propagation over a curved surface is well 
understood, problems still exist in terms of accurately predicting the magnitude 
and phase of the field, particularly near caustics such as the line of sight. The 
remainder of this thesis involves comparing newly acquired experimental data to 
the results of a refreshed and improved new version of the theory. The objective is 
to improve upon the accuracy of previous predictions and to sort out the important 
parameters and aspects of the scattered sound field near, around, and far from a 
curved surface. 
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CHAPTER III 


THEORY OF SOUND PROPAGATION 
OVER A CURVED SURFACE 


The objective of this chapter is to present all of the important ideas and as- 
pects of the current theory of interest as developed by Pierce. 12 Many of the 
mathematical details have been omitted but they can be found in Appendix A. 

3.1 Prototype Problem 

Consider the two-dimensional prototype problem shown in Fig.( 3.1). A plane wave 
of constant frequency, / = u/2n, (e~ iut time dependence) and with complex pres- 
sure amplitude Pit ik * is incident upon a locally reacting surface whose radius of 
curvature R is not necessarily constant, but is nevertheless everywhere large com- 
pared with Ijk. The surface is characterized by a finite impedance, Z,, which, 
for simplicity, is not a function of position. A cartesian coordinate system (x,y) is 
situated such that its origin is the point where a ray associated with the incident 
plane wave just grazes the surface, and such that the x-axis is coincident with this 
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ray. In this context, the origin is referred to as the apex of the surface since the rest 
of the surface lies in the half-space y < 0. 

The interest here is in obtaining a solution for the total acoustic field above and 
behind ( x > 0) the curved surface. From physical intuition, one can discern three 
basic zones surrounding the surface. As shown, these zones are labeled the bright 
zone , the penumbra, and the shadow zone. These terms derive from the analogous 
problem of a beam of light traveling past an obstacle. The term “penumbra” literally 
means a region of partial shadow between regions of complete shadow and complete 
illumination. 13 Although the sound pressure in the shadow zone is not exactly 
zero, it is small enough such that the term penumbra is still technically accurate. 
In this context, the penumbra refers to that region straddling either side of the line 
of sight. As in the optical analog, the width of the penumbra is flared as one moves 
further and further away from the apex. The penumbra is also often referred to as 
the transition zone. 

3.2 Geometric Acoustics Solution 

As a first step, one can apply the concepts of ray acoustics to derive a solution at 
points located in the bright zone (Fig. 3.2). This is reasonable since the propagation 
is considered to be relatively high frequency such that kR 1. The so-called 
geometrical acoustics solution is described as a superposition of an incident and a 
reflected ray, such that 

p = Pie ikx + P4A(0)/A(£)] 1/2 »e* fc *'e <kf (3.1) 

where p is pressure, P, is the magnitude of the incident plane wave, 9? is the reflection 
coefficient and A(i) denotes ray tube area after propagation a distance £ from the 
reflection point. The reflection point (x 0 >!/o)> the local angle of incidence 0, , the 
local curvature radius R, and the reflected ray path length £ can all be determined 
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Incident 
Plane Wave 



Figure 3.1 Depiction of the prototype problem for an investigation of long 
range sound propagation over irregular terrain. It is believed that the 
do minan t, effects on such propagation are due to interactions with the 
tops of hills or ridges. A generic hill top is approximated as a curved 
(parabolic) surface of finite impedance. The incident wave is plane. 
Three regions above and behind the curved surface can be identified. 
These are the bright zone, the penumbra, and the shadow zone. 
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for given listener coordinates (x, y) using the law of mirrors and the mathematical 
description of the surface. 

For the purposes of this study, consider that portion of the surface in the 
vicinity of the apex. One can approximate the local shape of the surface by a 
parabola y = -x 1 /2 R, where R is the local radius of curvature. In those cases 
where the reflected ray originates from this portion of the surface, the geometric 
acoustics solution can be written as 


P~P,C jl+l ZQ J l Q _l x+ £c_ R \ ] 

(3.2) 

where 


Q = [(4/9)x 2 + (2/3 )fly] 1/2 

(3.3) 

0 = (2 k/R 2 )[Q 3 - (8/27)x 3 - (2/3)i?xy] 

(3.4) 

Scaling parameters 


L x =R(kR )~ 1/ 3 L y = R{kR)~ 2/3 

(3.5) 


can be introduced such that, when p e~ ikx (as shown above) is expressed in terms 
of x/L x and y/Z,„, the resulting expression is independent of k and R, but not of 
their product. 

It is well known that geometric acoustic solutions such as the one derived above 
are invalid near caustics. In the prototype problem, the line of sight is one such 
caustic. Therefore, one does not expect good results from (3.2) at points on or near 
the line of sight. However, one recognizes that a solution (as yet unknown) valid 
within the line of sight caustic should have a reasonably smooth transition to the 
geometric acoustics solution in the region above and adjacent to the line of sight. 
This reasoning suggests that the geometric acoustics solution could be used as a 
radiation (boundary) condition when solving the Helmholtz equation in the caustic 
region. More specifically, since an impedance boundary condition is assumed for 
the surface, one can pose the problem of determining the field near and on the apex 
in terms of the Helmholtz equation with an impedance (mixed) boundary condition 
at the surface and a geometric acoustics radiation condition at points well above 
the apex. 



26 



Figure 3.2 Geometric description of the acoustic field at a point in the 
bright zone above a curved surface. The reflection occurs according to 
the law of mirrors. The effects at the point of reflection are described 
by the plane wave reflection coefficient and by the divergence of rays 
striking adjacent infinitesimal segments of the surface. The solution is 
a superposition of a reflected and direct ray. 
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3.3 Matched Asymptotic Expansions and the Field near the Surface 

To begin, as was done in deriving the geometric acoustics solution, the local shape 
of the surface near the apex is approximated by the parabola y = — x 2 /2 R. Then 
one transforms coordinates from the cartesian set (x,y) to a parabolic cylinder set 
(u,v) such that the new equation of the surface is v = 0. The relation between the 
two coordinate sets is described by 

x = u(l + [v/R]) 

(3.6) 

y = V (l +[v/2R])-u 2 /(2R) 

The scaling par am eters L x and L y can be introduced such that a set of scaled and 
non-dimensional parabolic cylinder coordinates (£,rj) emerge. They are defined by 


_ u fkR\ 

" R V 2 ) 


2v (kR \ 2/3 


V = ~R V~2" 


One can use the basic rules of orthogonal curvilinear coordinate transformations to 


express the Helmholtz equation 


(V 2 +Jk 2 )p = 0 


in terms of the £ and tj. Further, if one writes 

p = P < e <fctt e <4,/3 G(e,r?) 


and introduces a small parameter 


= (- 

\kR 


(3.10) 


then, after discarding higher order terms in the parameter c, one is left with the 


parabolic equation 


d 2 G ,3G, _ 

^r + , a£ + ’> G = 0 - 


This is the same differential equation found by Fock (see (2.5)). 


(3.11) 
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The impedance boundary condition at the surface is, in general, of the form 

Vp • n + % —=—p — 0. (3-12) 

z* 

where n is the unit vector normal to the surface. In the scaled parabolic coordinate 
system, the impedance condition takes the form 

^ + qG = 0 at r) = 0 (3.13) 

dv 


where 


1/3 


*-•(“) T (3H) 

is a nondimensional scaled admittance. One notes that since Re{Z, } > 0 for a 
passive surface, the phase of q is such that 0 < arg(g) < ir. 

The radiation condition is formulated from an asymptotic analysis of (3.2) 
when £ < 0 and rj » 1. The analysis suggests that the radiation condition be 


-«•/#«<«* j 1+ ^ 




— 2 £ + (£ 2 +317) 1/2 

3{t 2 +Zv) 1/7 

-2£+(g 2 +37 ? ) 1 / 2 +3i< 7 

2e + (e 2 +3r?) 


r 

7 + - ^ ) } 

— Ziq) J 


(3.15) 


where 

6 = -j£l + ^{£ s +(£ 2 +3>)) 2/2 } (3-16) 

This asymptotic form of (3.2) is required to match the solution of (3.11-3.14) at 
points where £ < 0 and rj 1. It can be shown that £tj — £ 3 /3 + 0 satisfies the 
eikonal equation. This fact suggests that (3.15) is of a reasonable form. 

Fock solved the boundary value problem posed above by the method of sepa- 
ration of variables. An alternative but similar approach is to introduce a function 
T(r/,a) and assume a solution for G in terms of the Fourier transform of T. That 


G(e,*l)« / e <a «T(» h a)da 

Jc 


is, write 


(3.17) 
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where the contour C is independent of £ and tj. The introduction of this transform 
into (3.11) separates the variables and results in the requirement that 


<PT 

dr/ 2 


+ (»? - a)T = 0 


(3.18) 


which is the Airy differential equation. This differential equation is the fundamental 
member of a class of differential equations known as turning point problems. As the 
name implies, the character of the solution of such problems changes when a so- 
called turning point is crossed. In (3.18), the turning point is r) = a. One notes 
that for T) > a, T will have a sinusoidal character, while for r) < a, T will have an 
exponential character. Two independent solutions of (3.18) are 


v(a — r/) = 7r 1/,2 Ai(a — tj) 

to, (a - r?) = «" / ‘2ir 1/, Ai (V ! ' ,3 (a - ij)) 


(3.19) 


where v(z) and Wy ( z ) are so-called Fock functions and Ai(z) is the Airy function of 
complex argument. The Airy function is an entire function and can be expressed 
as a contour integral of the form 


Ai (z) = —[ e i[,,/3+M,] ds (3.20) 

2 * Jc Al 

where the integration contour C Ai can be any path originating at infinity in the 
sector 27 t / 3 < arg(s) < n and ending at infinity in the sector 0 < arg(s ) < */3. 

The boundary condition (3.13) is easily satisfied by substituting in the assumed 
solution (3.17), letting T(rj, a) = v(a — tj) + c 0 (<*)«;! (a — tj), and solving for c 0 (o). 
The task of satisfying the radiation condition (3.15) is, at first glance, rather im- 
posing. However, physical reasoning requires that (3.17) approach the first term of 
the radiation condition when the listener is in the illuminated region (i.e. when £ is 
large and negative). The result of such a requirement can be shown to satisfy the 
entire radiation condition. (See appendix A.) 
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The resulting complete solution is 


« «■■»•«) = 7? /: - ”) - ^^ r (a - 4 '"' (3 - 21) 


which is also referred to as the Fock-vasi-der-Pol-Bremmer function and is trivially 
related to what Logan (1959) calls “Fock’s form of the van der Pol-Bremmer dif- 
fraction formula.” The two ter ms in the integrand are associated with the incident 
wave and the scattered wave, respectively. This association is borne out by the 
facts that the second term alone is a function of the surface property, q, and the 
first term is integrable as 





(3.22) 


A comparison of this result with (3.15) indicates that (3.21) is associated with the 
incident wave. 


3.4 Alternate Forms of the Fock-van-der-Pol-Bremmer function 

Although (3.21) is a general solution to the prototype problem, the integral itself 
is not easily computed, either analytically or numerically. The reason for this is 
that the integrand is often rapidly varying along the entire length of the negative 
real axis. However, alternate or approximate forms of (3.21) can be found for the 
penumbra, shadow zone, and along the surface such that better or more efficient 
numerical results can be obtained. The presentation of computational results is 
deferred until chapters 7 and 8. There, all of the theoretical predictions of the 
insertion loss at various field positions will be presented along with appropriate 
experimental results. In the remainder of this chapter, the various alternate forms 
of the Fock-vsni-der-Pol-Bremmer function, which were used during computations, 
are presented and discussed. 
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3.5 Field near the Apex and on the Surface 

At this point, it is helpful to introduce a third Fock function w 2 which is a linear 
combination of tvi and v of the form 


w 3 (z) = [zt-"* 13 ). (3.23) 


Then, (3.22) can be reexpressed without approximation as 

«*■•> ■ t- /." I* - « - 

where fj = r?e _,2,r/3 , q = qt i2lt/ 3 , and s = se <2,r/3 . This form is preferable to (3.21) 
because the integrands become vanishingly small for values of s greater than s « 8. 


Although several integration schemes were tried, satisfactory results were found 
using the extended Simpson’s rule. (In fact, in cases where the integrand is rapidly 
oscillating, Simpson’s rule can be more accurate than other more sophisticated, 
higher order methods.) For |£| < 2, a step size of As = 1/20 was used. For larger 
values of |£|, a step size of As = l/10|£| was used. 

On the surface, one can simplify (3.24) by multiplying the numerator and 
denominator of the first term of the integrand by the denominator of the second 
term and combining the two terms. Then one can make use of the Wronskian 


relations W(t>,ti>i) = W(v,tu 2 ) = 1 to yield 

f- e -t*</2 e - «{\/3/2 


1 r°° 

G(£’ °» 9 ) = ^ J o (s) g» 2*/3 

1 f°° 

\f* Jo W 'l ( 5 ) - 9 


qw 2 (s) 


- ds 


ds 


(3.25) 


qw i (5) 

Again, the extended Simpson’s rule was used to estimate the integral. In both cases, 
the integrands themselves were evaluated with the help of an algorithm developed 
by Pierce for the computation of Ai(z) where z is arbitrary complex. 14 This 
algorithm is presented in some detail in appendix B. 
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3.6 Field in the Penumbra 

In the penumbra, one can use the method of steepest descent to derive an approx- 
imate expression for G. This expression can be written as the sum of two parts. 
The first of these parts contains an integral of type commonly found in problems 
involving Fresnel diffraction, and is represented by the symbol $. For reasons which 
will be soon clear, $ is referred to as the knife edge diffraction term. The second 
part consists of two integrals which are collectively referred to as the background 
(the term used by Fock), and represented by the symbol Thus, if one writes 

G = $ + tf, (3.26) 


then it cam be shown that 


vW„ »;(*)-»««'/»»,(») v ’ 

1 f°° i. t V{s)-qv{s) , w 

_— / (\ w ' {s-v)ds 

V* Jo w i ( 5 ) “ 9 w i ( s ) 


(3.27) 


and 


$ = 


c -.«*/ 3 c «n h{Y) - ~j=e 
v2 


A D (Y ) 


(3.28) 


where H(Y) is the Heaviside step function. The argument Y, which represents a 
dimensionless height, is positive in the bright zone and negative in the shadow zone, 


and is written as 

r-GT^-n 

The diffraction integral, A D (V), is 

1 f +eo e~ u * du 

i */ 2 ) l ' 3 Y-e~<*'*u 


(3.29) 


(3.30) 


and is commonly found in asymptotic solutions of the wave equation in cases when 
at least part of the boundary has a sharp edge. 15 For purposes of computation, 
A d (Y) can be expressed as the sum of two auxiliary Fresnel functions each of which 
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has a simple rational approximation valid for positive values of Y . The fact that 
A d (V) is an odd function of Y enables one to use the rational approximations for 
both positive and negative values of F. 16 

As an example of another problem in which the diffraction integral, A D , arises, 

A 

consider the problem of a point source of strength 5 radiating over a rigid wedge of 
exterior angle /?, where 7 r < (3 < 2tt (Fig. 3.3). 17 When both kr and kr, are large 
compared to 1, when the source angle <f>, is between ir and /? but close to neither, 
and when the listener is near the line of sight, then an approximate expression for 
the complex pressure amplitude is 



H(X) - 



(3.31) 


where 

xL ) 

where L = r, + r. When the incident wave is plane, and the listener is near the 
line of sight, another dimensionless height can be derived from X. This derivation 
is accomplished by setting r « x, y « and letting r, — *■ oo. Then 


A 4> and A<f> = <f> — ( <j> , — ir) (3.32) 



/ ic \ 1/2 

x *y [— ) = * ( 3 - 33 ) 


Note that the wedge solution is not a function of the wedge angle /?, and is valid 
for any thin wedge or knife edge. One recognizes the similarity between (3.31) and 
(3.28); in fact, one can show that they are identical except for the difference in the 
dimensionless height variables X and Y. Thus, one can infer that similar physical 
mechanisms axe involved, and so it is appropriate to refer to $ as a knife edge 
diffraction term. Interestingly, $ is not a function of the acoustic impedance and is 
virtually independent of the shape of the diffracting body. (In fact, there is a slight 
dependence through the scaled coordinates.) This independence is in accord with 
physical intuition in that when the source and receiver are far from the surface and 
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Figure 3.3 Geometry associated with the diffraction of sound over a wedge. 
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near the line of sight, one expects the curved surface to be indistinguishable from a 
thin wedge or knife edge. One could also apply the concepts of the the Geometrical 
Theory of Diffraction (GTD) and reason that as the source and receiver move further 
and further apart along the line of sight, less and less of the surface near the apex 
is involved in the diffraction process. In the limit of such reasoning, the diffracted 
rays impinge upon the apex of the surface and depart from there without traveling 
any appreciable distance along the surface itself. In this way, the far field effect of 
diffraction from a curved surface is like that from a knife edge. 

According to the above reasoning, all curved surfaces eventually diffract as a 
knife edge in the limit of the source and listener being very far from the surface. 
Therefore, at listener locations in the penumbra and far from the apex, it is rea- 
sonable to propose (3.31) as the appropriate solution. It has already been proposed 
that a valid solution near the apex is G = $ + ^ (3.26-3.28). If these solutions 
are accepted, then several questions immediately come to mind. Namely, is there 
a smooth transition from (3.26) to (3.31) at points neither near nor far from the 
apex? If not, then clearly there is some mid-range region in which neither (3.31) nor 
(3.26) are valid or at least in which they are incompatible, and so the next obvious 
question is in what regions are they valid? Furthermore, in regions where neither 
is valid, how do the solutions behave and do they still yield reasonable results? To 
address these questions, one need first examine certain properties of (3.26-3.28). 

Computations show that the magnitude of the background, 'F, is significant 
at points close to the apex but less so at points far behind the apex, say x » R. 
At such points, one might be tempted to let G « $. Such an approximation is 
appealing because the work involved in evaluating 'F is much greater than that 
involved in evaluating $. Thus, a good deal of computational effort would be saved 
by neglecting 'F, when possible, and approximating G as $. 

At first glance, this result is also very appealing because the functional form of 
$ is identical to that in (3.31), and because (3.31) is valid at distances far from the 
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apex. Therefore, one might conclude (erroneously) that a smooth transition exists 
between (3.26) and (3.31), and that $ can be used as the extended solution to (3.26) 
at points in the penumbra far downstream from the apex. However, an inspection 
and comparison of $ and (3.31) reveals that although the form of $ is identical 
to that of (3.31), an important difference exists between the dimensionless height 
variables, Y and X. The difference between Y and X is a result of the different 
geometries involved in the respective problems. From symmetry considerations, the 
zero height line for a knife-edge must be the line of sight. In (3.32), one sees that 
X = 0 when A <j> = 0, which is the line of sight. However, Y = 0 is not coincident 
with the line of sight. Instead, relative to the (x,y) coordinate system, Y = 0 is an 
upward sloping curve (parabola) originating at the apex, which is the point (0, 0) 
in both the (x,y) and (£,»?) coordinate systems. Therefore, at points far from the 
apex, Y = 0 represents points of cartesian height, y » 0; that is, far from the apex, 
points of Y = 0 lie well above the line of sight. Thus, the divergence of Y = 0 from 
X = 0 at large x or r indicates that $, and thus also (3.26), will not asymptotically 
approach (3.31), the knife-edge diffraction solution. That the form of the knife-edge 
diffraction solution should fall out of the integral solution (3.21) would seem to be 
a consequence of the fact that near the apex, Y X. 

An interesting case of (3.26) is found when a value of q is such that ^ « 0 when 
Y « X (i.e. when x is small). In such cases, the field behind the curved surface will 
be nearly the same as if the surface were a rigid knife-edge. That such q's should 
exist is supported by again reasoning as in the GTD. For instance, when the curved 
surface is soft or highly absorbent (large q), one expects the insertion loss behind 
the surface to be larger than that behind a knife-edge because all rays impinging 
upon the surface and all creeping waves traveling along the surface will be highly 
attenuated.! In other words, there is a relatively large amount of absorbent surface 

f Iniert ion lo.. i* defined as that loss, expressed in dB, which occurs at a field point due 
to the insertion of an obstruction into the field. In the case of a single curved surface, 
the insertion loss is defined by IL = 20 log 10 |p/pi| • 
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area near the line of sight. Conversely, when the curved surface is very hard, one 
expects the insertion loss behind the surface to be smaller than that of a knife-edge 
because the creeping waves will propagate into the penumbra and shadow zone with 
very little attenuation. Geometrically, one senses that more of the incident wave is 
likely to scatter forward in the case of a curved surface. If this is true, then more of 
the incident energy is also spread over the area behind the curved surface. The fact 
that the area behind a curved surface is smaller than that behind a knife edge is 
important when one considers that no energy may cross an infinitely hard or rigid 
surface boundary. Thus, one expects a stronger overall field behind a rigid curved 
surface than behind a rigid knife-edge. Such reasoning indicates that a moderate 
surface admittance q (not too hard or too soft) might yield results very similar to 
that of a rigid knife edge even when the listener is near the apex. 


3.7 Field in the Shadow Zone 


Well below the line of sight, in the shadow zone, another approach to the solution of 
(3.21) is available. After some simple algebra, (3.21) can be rewritten in the more 
compact form of 


G{£,ri,q) 


1 r Mo=,*?) -<7X(Q»*7) ' 
n/W-oo L wj (a) - gw* (a) . 


e' ai da 


(3.34) 


where 


<p{a , tj) = w[ (a)v(a - rj) - v'(a)w i (a - tj) (3.35) 

and 


x(“>»?) = w l (a)v{a - rj) - v(a)wi(a - tj) 


(3.36) 
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In the limit rj = 0, x(«>°) = 0 identically, while the Wronskian relation W(v, u> x ) = 
1 yields ^>(a,0) = 1. So, for sufficiently small rj, the behavior of the integrand of G 
is dictated by e <0< /(«"i (a) - qw i (a) )• Further, at large a, 

< 3 ' 37 ) 

These facts allow one to close the integration contour along a semi-circle at infinity 
in the upper half of complex a plane without effecting the value of the integral. 
Then, the Cauchy residue theorem can be applied to evaluate it. The integrand has 
poles located at the roots of 

w[ (a) — qwi (a) = 0 (3.38) 


The fact that the Airy function and its associated Fock functions are entire ensures 
that the above roots are the only singularities of the integrand in the finite complex 
plane. Since all of the singularities are simple poles, the residue of the integrand is 


^(ttn) ~ ?v(o!n) 

<(«») ~ ?<(««) 


ti>i (ci£ n - r?)e’ a * t 


(3.39) 


where a n are the roots of (3.38). This expression may be simplified by substituting 
w"(a n ) = a n Wi (a n ) and w[ ( a n ) = qw x (a n ). Then, one can write down the residue 
series solution for the field in the shadow zone as 


v'(cO -gv(<*„) w 1 {a n - rj) ia ^ 
^ / \ C 

0£« 


(3.40) 


= -2* 1/3 iy' . , i 

v ' a n - q 2 «>i( a n) 

n = I 

Further, from the aforementioned Wronskian relation, and from the definition 
of the a n roots as solutions of tu'^otn) = gtt^a,,), it can be shown that t/(a„) — 
qv(a n ) = —l/wi (<* n ). Thus, a more compact expression for the above residue series 


solution is 
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G{Z , »?,?) = 2tt 1/2 i i ^2 7 ^ 3 7~ 

“ [«» -? 2 ]u; 2 («n) 


(3.41) 


In this way, each of the a„ poles is identified with a single creeping wave 
and the field in the shadow zone can be depicted as a superposition of an infinite 
number of creeping waves. Each creeping wave originates at the grazing point of the 
surface, travels along a geodesic of the surface, and continuously sheds rays into the 
penumbra and shadow zone. The shed rays leave the surface along tangent lines. 
Fig. (3.4) shows the location in the complex plane of the first ten a n poles for the 
plywood and carpet-on-plywood surfaces. The real and imaginary parts of a n are 
positive in accord with the notion that the creeping wave attenuates as it travels 
along the surface and as it sheds rays. When the listener is deep within the shadow 
zone, the series (3.41) is rapidly converging, and in some cases it suffices to keep 
only the first two or three terms. However, as the listener approaches the line of 
sight, more and more terms need to be kept. In general, as £ becomes small, many 
terms must be kept because the reduction in magnitude of succesive e' a ' ( terms 
becomes small. Also, as mentioned above, the solution relies upon a sufficiently 
small value of rj. It has been suggested elsewhere that £ > r? 1/2 is a necessary 
criterion for the application of this solution. The computational results based upon 


this solution are discussed in chapter 8. 

One clear advantage of this approach is that one need calculate the roots, a„ , 
only once for a particular q. This implies that only the functions ti/j (a n — r?) and 
e* a *' 7 need to be recalculated at different spatial points. This attribute offers the 
potential for significant computational savings over other forms of the solution. 
Thus, (3.41) is the equation of choice when the listener is located in the shadow 


zone. 
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Roots for CW Series 
o Carpet— on— Plywood 
x Plywood 


Freq = 10 kHz 


Real Part 


Figure 3.4 Location of the first ten roots of w[ (a n ) — qw x (o n ) = 0 for the 
cases of the plywood and carpet-on-plywood surfaces. 
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3.8 Evaluation of the Zeros, z„, involved in the Creeping Wave Series 

The only remaining difficulty is to calculate the poles a n which, of course, are the 
roots or zeros of (3.38). In order to do this, it is helpful to rewrite (3.38) in terms 
of the Airy function itself and then to note certain properties of the new equation. 
Letting a n = z n c~ i7 *t 3 , the equation for the zeros, z „ , is 

Ai' ( z n ) - pAi(z„ ) = 0 ( 3 - 42 ) 


where p = qe ~ i2 " / 3 . It has been stated that 0 < arg(g) < it. Thus, — 2 tt/ 3 < 
arg(p) < jt/ 3. In the following discussion, no restriction on the phase of p is 
assumed. However, in the end, the above restricted phase range of p proved to be 
important because, for such a range of p, all of the a n zeros were found to lie in 
the first quadrant. 18 

The Newton-Raphson method can be used to solve numerically for the z„ zeros. 
As will be shown later, the zeros of Ai(z) and Ai* (z) are well known. Let z n be the 
n th zero of Ai'(z). Then to find a root of (3.42) for small p, let z n = z' n + A n , and 
write 

/(An) = Ai'« + A„) -pAi« + An) = 0 (3.43) 


where A„ is a small and generally complex number. Application of the Newton- 
Raphson method of solving for A„ results in 


Ai* (z„ + AjJ — pAi(z^ + A|J 
(z; + Aj> — p 2 )Ai(zJ, + A|J 


(3.44) 


where A 0 = 0. When p is not small, the z„ zeros still can be found by using 
(3.44) recursively. That is, when p is not small, then a sequence of intermediate 
zeros, z„, which are zeros of (3.42), can be found for a sequence of intermediate 
values of p. If the sequence of intermediate p values converges upon the true value 
of p, then the sequence of intermediate zeros, z„ should converge upon z n . For 
example, a sequence of intermediate zeros could be found which correspond to a 
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sequence of intermediate p values, all of which have the same phase as the true 
value of p. The method consists of picking an initial p whose phase is that of the 
true p but whose magnitude is small enough that the Newton-Raphson method 
converges upon an associated intermediate zero, z n . The new z n can be substituted 
for z' n in (3.44). Another intermediate p, whose phase is that of the true p, but 
whose magnitude was chosen to be incrementally larger than that of the previous 
intermediate p, can be inserted into (3.44). Then, the Newton-Raphson method will 
yield a new intermediate zero, z n . In this way, for a given phase of p, the Newton- 
Raphson method can be employed to march along the n root (zero) trajectories 
from 0 — * |p|.f The same method can be used when |p| is large except that z' n 
should be replaced by in (3.44), where z° Me the zeros of Ai(z). The difficulty 
here is that, when |p| is not large or small, the method requires a relatively large 
number of evaluations of the Airy function and its derivative. These evaluations 
require a considerable amount of computation. Another approach is to find a z* 
such that, if the Newton-Raphson method, beginning at z‘, were applied to (3.42), 
the method would always converge upon z n and do so in relatively few steps. Such 
a number as z* is referred to as an approximate zero. The remainder of this chapter 
is concerned with finding approximate zeros of (3.42). 

To begin, an idea of the nature of the zeros z n can be obtained by examining 
certain properties of the Airy function and its derivatives. Toward this end, one 
notes that Ai(z) = Ai(z) (an example of Schwarz’s principle of symmetry), where 
z denotes the complex conjugate of z. Thus, it is also true that z„(p) = z n (p). 
So, for arbitrary complex p, one need only consider 0 < arg(p) < i r. One also 
notes that Ai(z) and its derivative are real when z is real. This is easily seen from 
the fact that the power series expansion of Ai(z) has all real coefficients and is 

f In this context, a root trajectory ii a path in the complex plane traced out by the 
continuous sequence of roots associated with a continuous variation of |p( from sero to 
infinity for a given phase of p. 
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absolutely convergent for all z (See appendix B). Thus, when z n is real, p must be 
real. Conversely, since z n (p) = z n {p), and also since p = p when p is real, then 
z n (p) = z„ (p) when p is real. This can only be true if z n (p) is real. Since the 
statement and its converse are true, then the contrapositive and its converse are 
true, namely, that when p is complex, z n is complex and vice versa. 

Another useful fact is that all of the zeros of Ai(z) and Ai*(z), labeled z° and 
z' n , respectively, are real and negative. These zeros are interlaced along the negative 
real axis such that 

0 > z\ > z° > z' 2 > z\ > z' 3 > . . . (3.45) 

The first ten of these along with an asymptotic approximation for the rest can be 
found in various mathematics handbooks. 19 One also notes that when p = 0, 
z n = z' n . Further, since 

Ai"(0) = 0 Ai'(0) < 0 Ai(0) > 0 (3.46) 


one can conclude that, for real p, as p — ► +oo, the zero z n travels from z' n to 
z°, and that as p — * — oo, the zero z n travels from z‘ n to z°_ x . In the case of 
n = 1 and p -*• -oo, the zero z n travels from z[ to +oo where Ai(x) and Ai' (x) 
asymptotically approach zero. Thus, the entire real axis is divided into segments 
of z n root trajectories, each segment occupying that part of the real axis between 
successive points in the sequence 


(oo, 


z* ^ z* z ^ z 1 




(3.47) 


The segments correspond alternately to negative and positive values of the real 
parameter p. 

When p is complex, it was shown that the zeros z n are always complex. This 
implies that the complex root trajectories never cross the real axis. However, the 
complex root trajectories converge upon the real axis when |p| — ► 0 and |p| — ► oo. 
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What are the complex root trajectories? To get at that question consider each 
root z n to be a function of the parameter p ; that is, write z n = z n (p) . Then 
differentiate (3.42) with respect to p, and employ Ai "(z) = zAi(z), to show that 

dzn _ 1 

dp z n p 2 



Immediately, one notes that z n = p 2 is not a valid root. This result is fortunate 
since otherwise (3.41) could be unbounded. (This result can be proven formally 
with the help of the Argument Principle of complex analysis.) The shape of the 
trajectories can be ascertained through inspection of (3.48). Let p = |p|e‘ 9 . Since 
z n = z' n when |p| = 0, consider how the trajectories leave z' n when |p| is small. In 
this case, one can write 


dz n e i{e - n) 

d \p\ w ki 


, when |p| < | z' n 


(3.49) 


So as 0 = arg(p) goes from 0 — ► 7r, arg(dz n /d|p|) approximately goes from —n — ► 0, 
for |p| < |<|. For instance, when 0 = 0 , dz n /d\p\ is negative and real, and when 
0 = it, dz n /d\p\ is positive and real, as was deduced previously. Similarly, since 
z n = z° when |p| — ► oo, consider how the trajectories approach z° when |p| is large; 
in particular, when |p| |z°|. In this case, 


dZn _ 

d\p\ ~ \p\ 2 


, when |p| > z® 


(3.50) 


So as 0 = arg(p) goes from 0 -*• n, arg(dz„/5|p|) approximately goes from tt -* 0, 
for |p| » |z°|. 

With similar reasoning, one can sketch out the approximate trajectory for any 
0 as |p| goes from 0 — * oo. With the exception of those trajectories emanating from 
z^ and heading toward positive infinity, the result is that the trajectories are the 
shape of smooth convex curves connecting adjacent points in the sequence (3.47). 
In fact, when |0| or \n — 0\ is less than 27r/3, the curves are roughly circular arcs. 
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The trajectories all lie in the second or third quadrants, and near the negative 
real axis when |0| is not near x/2. Of course, each of these trajectories has an 
image trajectory located in the other quadrant associated with p. Furthermore, the 
trajectories which move the furthest from the negative real axis before returning are 
those trajectories associated with 9 « 7t/2 (or —x/2). Somewhat curiously, when 
one considers all of the trajectories which originate at z' n , 9 = x/2 (or —x/2) is not 
the transition value of 9 below which the trajectories travel to z° , and above which 
the trajectories travel to z°_ x . The actual transition value is somewhat greater than 
9 = x/2 (or less than 9 = —x/2) such that p 2 lies in the third (or second) quadrant 
near the negative real axis; that is, p 2 is potentially near to z n when 9 = x/2 + 6 
(or when 9 = —x/2 — S), where 6 is a relatively small positive number. For such 9 , 
dz n /dp does not exist as seen from (3.48). 

A power series representation of z n (p) can be constructed for the case of |p 2 | < 

| z n | by substituting 

oo 

z„ = £a fc p fc (3.51) 

Jb = 0 

where Oq = z' n . One finds that a 1 = 1 jz‘ n and the remaining coefficients are 


a k+l 


(fc ~ l)flfc-l ~ Ol Ofc 

Oo(fc + 1) 


yi (m + 1) a m + ! a* - m 
— ; (* + !) °o 

m = 1 


fc = 1,2,... (3.52) 


The first two terms of (3.51) comprise the often stated approximation of z n for small 
p. A complementary power series representation of z„ ( p ~ 1 ) can be constructed for 
the case of Ip 2 1 > \z n |. The series can be written as 


o 1 V- bt 


<= 3 


The coefficients are found from 


1 * -3 

b t = ( m + l)bt- m - 3 b m + 1 


m ~ 0 


(3.53) 


i = 3,4,... 


(3.54) 
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The series (3.53) does not converge except for very large |p|; and although (3.51) 
converges, it encounters difficulty for moderate values of |p|. However, the utility 
of these series is in their ability to generate approximate zeros z* (p) for z n (p). 

Consider the problem of finding z x (p) when arg(p) = 8.1° (as in the case 
of the plywood surface) and when arg(p) = 46.5° (as in the case of the carpet- 
on-plywood surface). Fig. (3.5) shows the root trajectories as calculated by the 
Newton-Raphson method at small increments A|p| = 0.02 from |p| = 0 — + 10. 
(This same trajectory can be obtained by solving (3.48) with a standard Runge- 
Kutta method for ordinary differential equations.) Also shown are the approximate 
root trajectories obtained by successively including more and more terms in (3.51) 
and (3.53). The figure can be deceiving in that the point for a given |p| on an 
approximate root trajectory tends to be further along the trajectory than is the 
actual root. However, the trick is that the approximate root trajectories provide 
the location of a point or points from which the Newton-Raphson method can then 
converge to the actual root. Convergence tests were run where the number of series 
terms kept was varied along with the value of arg(p). The best results of these tests 
were obtained when 6 terms were kept in (3.51) and when 4 terms were kept in 
(3.53). Thus, a reasonable conclusion was to limit the series to a finite number of 
terms and calculate approximate roots from the abbreviated series 

e 

z\ = ^a fc p k when |p 2 | < \z n \ (3.55) 

* = o 

and 

<=*H + ; + £jr when Ip 5 I > \ z ~ I ( 3 - 56 ) 

p tt p 

The convergence properties tend to hold true for a wide range of arg(p) and for any 
n of z n . Thus, a relatively fast way of finding z n when p is not small or large is to 
approximate the z n by (3.55-3.56) and then apply the Newton-Raphson method to 
polish the zero. 
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Real Part 


Figure 3.5 Example root trajectories when arg(g) is equal to that for the 
plywood and carpet-on-plywood surfaces, respectively. The trajecto- 
ries rolling over to the right are associated with the plywood surface. 
The discrete symbols, (□), represent the actual root trajectories found 
through a recursive application of the Newton-Raphson method. The 
solid lines represent the approximate root trajectories predicted through 
the calculation of the first few terms of the series solutions (3.51) and 
(3.53). J indicates how many of the first j series terms were used. 
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Exceptional cases may exist in which the optimum or required number of terms 
kept in the series (3.55-3.56) is different. As the accuracy of the < approximations 
worsens, the propensity of the Newton-Raphson method to diverge increases. Safe 
guards should always be included to assure that the Newton-Raphson method is con- 
verging. When the Newton-Raphson method diverges from a particular z ' n , one can 
almost always either increase or decrease the |p| used to generate the z n such that 
a new z‘ n is generated for which the Newton-Raphson method does converge. One 
can then increment |p| and return to the |p| of interest with the Newton-Raphson 
method. In this way, even if (3.55-3.56) do not generate a viable approximate zero 
for a particular |p|, they are likely to generate a viable approximate zero for a \p\ 

nearby. 

A word of caution is in order concerning the Z\ zeros when arg(p) ~ 27r/3 or 
greater but less than or « 4tt/3. In these cases, the magnitude of the Zj zeros 
grows without bound as |p| oo. Although it seems reasonable from the previous 
discussion of the general characteristics of all zeros of (3.44) that z x zeros do exist 
for all p, it may become increasingly difficult to find these zeros for the phase range 
of p stated above. This is particularly evident when one realizes that the distance 
between z x (p) and z a (p+ A p) grows as p ^ oo. Hence, the sensitivity of the method 
to errors in the starting point increases as |p| increases. This behavior is limited to 
the zeros (i.e. the zeros of the trajectories that have an endpoint at z' ). 

3.9 Summary 

In s umm ation, when the listener is well into the bright zone, the well known 
geometric acoustics solution 


p = Ptc““ + 


(3.1) 
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is valid. In cases where the point of reflection is not too far from the apex, an 
approximate forn of (3.1) was shown to be 


P<t 


ikx 


{ 


1 + 


r gj-jx i i/’ r g-Jx-j^R i 1 

3 Q J iQ-lx+j^Rl / 


(3.2) 


where 


Q = [(4/9)x 3 + (2/3) J Ry] 1 / 2 (3.3) 

0 = (2 k/R 2 )[Q 3 - (8/27)x 3 - ( 2/3)Rxy ] (3.4) 


This solution was used as the outer boundary condition to solve an asymptotic form 
of the Helmholtz equation. The resulting solution, valid on and near to the curved 
surface, was termed the Fock-van-der-Pol-Bremmer function 


- T. L - ’> - roSs- 1 * - 4“ *■ <“> 


and was considered the general solution to the problem of diffraction of a plane wave 
over a curved surface of finite impedance. Several alternate forms of this equation 
exist. In the penumbra, another form was shown to be 


G = $ + tf, 


(3.26) 


where 


= 


1 /■“ 

V* Jo 

1 f°° 

V* Jo 




i Jr/I 


t>'(s) — qe i2w J*v(s) 
tu' (s) — qt i7 "l z w 2 {s) 


w 2 (s — rje l2n ^ 3 )ds 


i.t v'(s)-qv{s) 

c ~rH ( s ~ ^) ds 

ty((s) - ?tvi(s) 


(3.27) 


and 


$ = 




H(Y) - ^e'i 


"‘'’A C (Y) 


(3.28) 
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Thus, when the listener is near the apex, the solution can be interpreted as the sum 
of a knife-edge diffraction term and two so-called background terms. The back- 
ground terms eventually vanish as the listener moves away from the apex. How- 
ever, at such locations, it was argued that the dimensionless height Y is probably 
no longer appropriate. Instead, an alternate dimensionless height, such as 

1/3 


* 


-(£) 


( 3 . 33 ) 


may be substituted for Y to yield the well known knife-edge diffraction formula for 
listener locations far from the apex. In the shadow zone, an appropriate form was 
shown to be 

u, ‘ (a ;~'! ) / -* (3.4i) 

v ’ “ [a n -<7 2 W(a n ) 

Each individual term of the series can be interpreted as a creeping wave such that 
the field in the shadow zone is seen as a superposition of such waves. Typically, 
only a few of these terms are needed to describe the field in the deep shadow. The 
series is believed to converge when £ 2 > r). 
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CHAPTER IV 


EXPERIMENTAL FACILITY AND PROCEDURE 


Scale model experiments were conducted as part of the study of propagation of 
sound over curved surfaces and its relation to the more complex problem of long 
range outdoor sound propagation over undulating terrain. This chapter contains 
a description of the laboratory and experimental methods used in the study. In 
general ter ms , the laboratory consisted of a large table on the center of which rested 
one or two hollow cylindrical ridges. Side by side, the ridges spanned the width of 
the table. A spark generator, which was located near one end of the table, was 
used as an acoustic source. Near this same end of the table, a reference microphone 
was given a fixed (permanent) position. This microphone was positioned such that 
it could receive the incident acoustic pulse prior to the arrival of any reflections. 
Another microphone was used to make field measurements of the acoustic pulse after 
it had been reflected or diffracted by a ridge. The microphones were connected via 
a preamplifier, power supply, and amplifier to an A/D data collection board and 
associated personal computer. Data reduction was performed by software stored in 
the personal computer. Of interest was the influence of the ridges upon the acoustic 
signal in terms of insertion loss. The insertion loss was inferred from the ratio of 
the digital Fourier transform of a pulse received by the field microphone to that of 
a pulse received by the reference microphone. The details of this procedure and the 


laboratory follow. 
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4.1 Laboratory Design 

The laboratory was housed within a standard university small laboratory room 
of dimensions 8.2 m by 6.1m by 4.3 m and located on the campus of the Georgia 
Institute of Technology. The walls of the room were cinderblock, the floor was 
tiled, and some of the walls were lined with shelves. In no sense did this room 

approximate the ideal of an anechoic chamber. 

Four tables were made to be used for the scale model experiments. Each 
table was 1.2 m wide by 2.4 m long and 0.9 m high with a table top made of 3/4 in 
(« 2 cm) CDX plywood. Each small table frame was constructed of two-by-six (5 cm 
by 15 cm) yellow pine grade #1 planks; the table legs were constructed of four-by- 
four (10 cm by 10 cm) yellow pine frade #2 beams. The fasteners holding the table 

together were machine bolts and wood screws. 

The table frame was made by running two 2.4 m length two-by-sixes parallel 
to each other, 1.2 m apart. Then five equally spaced two-by-sixes were mounted 
in between the first two. Next, a shelf was cut into the four-by-fours so that they 
would fit into the corners of the frame and still leave some of the frame resting on 
the shelf. Each four-by-four was bolted into the frame with three machine bolts. 
Finally, the CDX plywood was placed on top of the frame and secured with wood 
screws around the perimeter. The four small tables were bolted together, forming 

one large table 4.9 m long by 2.4 m wide. 

Two identical curved surfaces were constructed to be mounted on the large 
table and used as laboratory scale models of topographical ridges. Each surface 
was formed by bending a 1/4 in CDX plywood sheet, measuring approximately 4 ft 
wide by 8 ft long, over two exterior plywood ribs. The ribs were cut such that their 
top edges were circular arcs 2.5 m in radius. These ribs were placed parallel to 
each other and equidistant (« 17 in) from the lengthwise centerline of the plywood 
sheet. Thus, two cylindrical ridges approximately 32 cm high and 235 cm along the 
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base were formed. Each ridge was mounted upon a separate sheet of 3/4 in CDX 
plywood which was the same size as the top of one of the four smaller tables. A 
ridge (s) were placed upon the table such that the lengthwise centerline of the table 
top and the centerline arc of the ridge(s) were coplanar. This plane was termed the 
centerline plane. 

Large rectangular strips of a 1/4 in thick commercial carpet were cut so as to 
fit atop the plywood curved surface. These strips were later placed on top of the 
plywood curved surface(s) to provide a second type of surface impedance distinctly 
different from that of the bare plywood surface. The carpet had a hemp mesh 
backing interspersed with glue which held the carpet fibers in place. 

Two Bruel & Kjaer 4136 quarter-inch condenser microphone cartridges were 
used for making the sound pressure measurements. Each microphone cartridge was 
attached to a Bruel & Kjaer 2615 cathode follower (pre-amplifier) via a Bruel & Kjaer 
UA0035 adaptor. The cathode follower and the microphone cartridge were powered 
by a Bruel & Kjaer 2801 microphone power supply. This assembly was referred to 
as a microphone. Each microphone cartridge had a dynamic range rating of up 
to 180 dB, and a rated sensitivity of « 1 mV/Pa. Further, the diaphragm of each 
cartridge was approximately 6.35 mm in diameter, and was covered by a protective 
grid housing which screwed on to the tip of each cartridge. More details about these 
microphones can be found in the manuals available from Bruel <k Kjaer. 

One of these microphones, referred to as the reference microphone, cartridge 
serial #1200301 with open circuit sensitivity of 1.26 mV/Pa, was positioned such 
that the sound arriving at it during the time interval associated with the direct wave 
arrival was as close to free field (no reflected arrivals during the time interval of the 
direct arrival) as possible and such that the microphone casing and suspension did 
not impede or interfere with propagation of the acoustic wave to the field micro- 
phone. Specifically, the reference microphone was suspended at a height roughly 



54 


three times the height of the ridge (32.4 cm) by three separate lines of 1/16 in nylon 
cord. The other microphone, referred to as the Geld microphone, cartridge serial 
#1200305 with open circuit sensitivity of 1.32 mV/Pa, was mounted upon a Fisher 
Scientific Flexifraune portable stand by way of a burette clamp. The microphone 
was held by the clamp such that the microphone diaphragm was roughly 9 cm for- 
ward of the clamp itself. When needed, a 6 ft long, 5/8 in aluminum tube acted as a 
horizontal extension of the clamp, the microphone held within and protruding out 
of the tube. 

Each microphone was connected in line with its own low current amplifier. Each 
amplifier consisted of a Motorola LF351N FET operational amplifier microchip, 
which had a high voltage slew rate of 13 V/fis and a flat response up to 100 kHz. 
The am plifier was designed to produce an amplified, non-inverted signal. A set of 
nine possible gains was achieved by the inclusion of an array of feedback resistors 
controlled by an external multi-position switch. Gains from unity to 100 were pos- 
sible. In all of the subsequent experiments the op amp position switch was set to 
either #5 or #6. These two settings corresponded to measured gains of approxi- 
mately 5.8 and 11.1, respectively. Each amplifier was powered by two parallel sets 
of four 9 V batteries. The transient level of the parallel sets of batteries was guar- 
anteed by parallel 0.1 /xF capacitors. The amplifier was enclosed by an aluminum 
box which was grounded to the batteries. The box served as a precautionary shield 
against electromagnetic noise generated by the spark source or other instruments. 

The amplified analog signal of each microphone was converted to digital form 
by an integrated hardware and software system produced by RC Electronics Inc. 
and called “Computerscope ISC-16.” This system consisted of a 16 channel A/D 
board, which was inserted into the IBM PC, an external instrument interface, and 
the scope driver software. The system was capable of recording an input voltage 
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signal with a peak to peak range of 20 volts centered at zero and of resolving it to 
12 bit accuracy, or equivalently to approximately 1 part in 4000. 

Thus, the data acquisition system was composed of two microphones, two am- 
plifiers, an analog-to-digital converter, and an IBM personal computer (Fig. 4.1). 
The system was capable of gathering data at an aggregate rate of 500 kHz and pos- 
sessed an aggregate memory buffer of 64 kilobytes. Software purchased from RC 
Electronics Inc. enabled the digitized data to be downloaded into standard ASCII 
coded files which were subsequently processed by the PC. 

The sound source was an electric spark generator (Fig. 4.2). The design of the 
spark generator was inspired by a photo of an apparatus designed by Dr. Mendel 
Kleiner at the Chalmers Institute of Technology. The important features of the 
spark generator were as follows. A gap of approximately 1 mm separated two stain- 
less steel rods or electrodes. (Originally, copper electrodes were used; however, 
problems were encountered with the tips either exploding or ablating away.) The 
electrodes were pointed at each other forming a spark axis. Each electrode was ap- 
proximately 1.5 mm in diameter. A third electrode was positioned midway through 
the gap and pointed perpendicular to the other two electrodes. All three electrodes 
were held in place by a narrow, wedge-shaped plexiglass block. (The shape was 
chosen to minimize reflections from the block.) Each of the first two electrodes was 
connected by an insulated copper cable to one point of a 1/xF, 10 kV (nominal) 
capacitor. One point of the capacitor was grounded while the other was connected 
to a 10 DCkV variable power supply through a 10 W resistor, nominally rated at 
10 Mfl. The third electrode was connected to a standard automobile ignition coil 
by way of a lOkfl resistor. 

The spark source was oriented such that its axis was perpendicular to the 
centerline plane and such that the center of the spark gap lay m that plane. In 
this way, the output of the source was expected to be circularly symmetric in the 
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Reference 
Spark Microphone 



Figure 4.1 Schematic of the interior of the laboratory room used in the 
study. Shown is the sequence of data processing associated with a single 
firing of the spark source. Reference and field microphones were sta- 
tioned in front and behind of the ridge, respectively. A triggering device 
activated the data capturing system simultaneous to the firing of the 
spark source. The analog signals of the two microphones were sampled 
at a rate of 250 kHz, digitized to 12 bit precision, stored in a 64 kilobyte 
memory buffer, and low pass filtered. Interesting portions of the total 
data field were then input to a discrete Fourier transform program. The 
output of this program was interpreted as the Fourier components of 
the original waveform portions. 
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Figure 4.2 Photo of the spark source electrodes and the wedge-shaped 
plexiglass block which held them. A gap of approximately 1 mm sep- 
arated two of the 1.5 mm diameter stainless steel electrodes. A third 
electrode delivered a low current, high voltage electric pulse to the gap. 
The subsequent ionization of the air in the gap enabled the charge stored 
in the 1 fi F capacitor to travel to ground. 
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centerline plane, and also symmetric about that plane. The nature of a typical 
waveform produced by the spark source is discussed in detail in the next chapter. 
A typical waveform had a peak acoustic pressure on the order of 130 Pa at 1 m from 
the source, a pulse duration of approximately 40 /is, and which was broadband from 
5kHz to at least 50kHz (Fig. 4.3). 

The spark discharge was initiated by a manual trigger attached to the ignition 
coil. The coil delivered a low current, high voltage (« 10-15 kV) electrical discharge 
through the tip of the third electrode. This discharge ionized the gap separating the 
first two electrodes, which were the power electrodes. Once the gap was ionized, its 
electrical resistance plummetted, and the charge stored in the capacitor traversed 
the gap to ground. During the transition to ground, some of the electrical energy 
was converted into heat. The subsequent conduction of this heat resulted in an 
expansion of the air around the gap, and, thus, an acoustic pulse as well. An 
antenna situated near the spark source transmitted the spark’s electromagnetic 
pulse to an external TTL port associated with the A/D board. Data collection was 
triggered by the reception of this pulse. 

A model wedge was also constructed with the anticipation of comparing the 
experimental data associated with a curved ridge to that associated with a wedge. 
Basically, the wedge consisted of a plywood body, made of two planks of approx- 
imate dimensions 3.5 ft x 12 in x 1/4 in, with a stainless steel tip formed from a 
3 ft x 1 ft sheet. The planks were joined such that they shared a common lengthwise 
edge and such that there was a 34° angle between the planks. This position was 
fixed by attaching the unjoined lengthwise edges to a plywood base. The stain- 
less steel sheet was creased along its lengthwise centerline, then folded there, and 
placed atop the shared edge of the planks. In this way, the crease formed the apex 
of the wedge. The entire wedge received two coats of spray paint. The wedge was 
approximately 13 in high with a base of 8 in. 
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Figure 4.3 A typical waveform at a distance of lm from the spark gap, 
and the results of a discrete Fourier transform of the waveform data. In 
general, such waveforms had a duration of 40 /xs, and a peak acoustic 
pressure of 130 Pa. The transformed data revealed a broad frequency 
band of Fourier components of the original waveform. 
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4.2 General Procedure 

The general experimental procedure was similar for all of the experiments reported 
here. Most of the experiments used two microphones, one at a reference point in 
the free field, and the other at a desired field point. The reference microphone was 
oriented such that the normal to its diaphragm was in the centerline plane and, in 
addition, was pointed directly at the spark source gap. A microphone alignment 
or mounting of this type was referred to as normal. The field microphone was also 
aligned normally, although the normal to its diaphragm did not always point directly 
at the source. (A transverse alignment was when the normal to the microphone 
diaphragm is perpendicular to the centerline plane.) The two corresponding analog 
voltage signals, which were output from the microphones, were amplified, sampled 
at intervals of 4.0 *ts, and digitized to 12 bit precision. A constant was automatically 
added to each data sample such that the mean of the data sample was approximately 
zero. Since the voltage increment registered by a microphone was opposite in sign to 
that of the corresponding pressure increment, the sign of the shifted digitized data 
was reversed. Then the data sets were filtered by a lowpass digital filter with an 
upper cutoff frequency of 75 kHz. A waveform portion of interest and representative 
of what would be received if there were no undesired reflections contaminating the 
data was selected. Data points outside these portions were then replaced by zeros. 
Then, two abbreviated time windows were manufactured such that each waveform 
portion of interest was centered in its new window. These windows had the same 
duration (usually, 500 x 4.0 ms), but usually not the same time beginning. In this way, 
two “windows” containing a reference microphone waveform and a field microphone 
waveform, respectively, were generated. The replacement of the extra data points 
by zeros was in accord with the expectation that an acoustic pulse has negligible 
residual effect on the ambient pressure, and that the acoustic pressure before and 
after each waveform would ideally be zero. The digital Fourier transform of each 
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window was computed and interpreted as if the value of the acoustic pressure were 
identically zero before the start and after the end of each window. 

The Fourier transform algorithm used in the subsequent calculations was taken 
from a listing of a Fortran program developed by Norman Brenner and Charles 
Rader of the MIT Lincoln Laboratory circa 1967 and which, according to the com- 
ment statements included in the program listing, had some relationship to what was 
described in a paper “Fast Fourier tranforms for fun and profit” by W . Gentleman 
and G. Sande, presented at the 1966 Fall Joint Computer Conference. As used in 
the data reduction process of a generic experiment, the program returned a set of 
N numbers <j(n) from an input set v(n) according to the relation 


N 


m = 1 


Thus, if the input set was described by 
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From these relations, the value of the Fourier transform was identified once one 
knew the time t 0 at the start of each record, the number N samples, the sampling 
interval At, and the apparent transduction constant K (volts per pascal). In all 
cases, the actual measured time domain data occupied less than 20% of the 500x4/is 
total time domain duration and, as mentioned above, was centered in the 500 point 
window. For this reason, no adjustment, such as obtained through the application of 
a H amm ing window, was made to the finite time window. For time windows such 
as these, one can show that such an adjustment has a negligible effect upon the 
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Microphone Calibration History 
(with B&K 4220 Pistonphone) 

Date 

Reference Microphone 

Field Microphone 

Op Amp 
Setting 

K 

(mV/Pa) 

Op Amp 
Setting 

K 

(mV/Pa) 

12/22/86 

5 

3.015 

5 

4.776 

2/10/87 

5 

3.000 

5 

4.865 

6/22/88 

5 

2.980 

6 

9.425 

9/10/88 

6 

5.542 

6 

9.297 

9/14/88* 

5 

2.890 

5 

7.724 

9/14/88* 

5 

3.027 

6 

15.120 

1/6/89* 

5 

3.015 

6 

16.144 

8/12/89 

5 

3.135 

6 

15.433 


Table 4.1 Values of the transduction constant, K, at various dates. K 
was deduced through the use of a B& K pistonphone. The “Op Amp 
Setting” number refers to one of the nine possible gains of each operational 
amplifier. The gain is included in the transduction constant. The asterisk 
denotes the installation date of a new preamp for the field microphone. 


transformed data. The transformed data was interpreted as the discrete frequency 
distribution of the transient signal. Such frequency data was obtained at 500 Hz 
intervals from 0 to 125 kHz. 

The transduction constant was determined through the use of a Bruel k Kjaer 
4220 pistonphone. The pistonphone was calibrated by Bruel k Kjaer in July of 
1986 such that it produced a nominally pure tone of 250 Hz and 123.9 dB re 20/iPa 
(or ss 31.34 Pa rms) at the coupler termination. A transduction constant was de- 
termined by capturing a minimum of 40 cycles of the pure tone with each micro- 
phone, calculating an RMS voltage over these 40 cycles, and dividing this voltage 
by 31.34 Pa rms. Table 4.1 shows the value of the transduction constant for each 
microphone transmission line at various dates. 
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In many of the experiments, the objective was to obtain a measure of the 
insertion loss at a given field point, where the insertion loss was defined as 


IL = -20 log! 0 


pM 

Pi{m) 


(4.5) 


where p(m) and Pi (m) were equal to A(m)q(m) for the field and reference mi- 
crophones, respectively. Here, A(pi) included the adjustment for the frequency 
response of each microphone as well as adjustments for absorption in air and spher- 
ical spreading. Absorption in air was calculated from the ANSI standard formula 
where the pressure decrease due to absorption is written as 


P = Poe 


(4.6) 


Equivalently, the change in sound pressure level due to absorption is 


AIL, - —20 log 10 


_P 

Po 


8. 69 as, in dB. 


(4.7) 


The variable, s, is the propagation distance, in meters, while the parameter a is 
defined to be the attenuation coefficient due to atmospheric absorption and has 
units of (Np m" *). (The neper, Np, is the natural logarithm analog to the decibel, 
dB, such that 8.69 dB = lNp.) a can be written as the sum of four terms 


a = a et + a rot + a vibiQ + a„ i 6 i j/ 


(4.8) 


where a el represents the classical absorption, a rot represents the molecular absorp- 
tion due to rotational relaxation, and ct a n,,o an d &vib,N represent the molecular 
absorption due to vibrational relaxation of the oxygen and nitrogen constituents of 
air, respectively. With this notation, one can write 

1.84 x !<.-/» (*f) (£) ' 


(4.9) 
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and 


- 2(aA )"“-‘ (0 i +(///,,: 


(4.10) 


where t = O for oxygen and * = N for nitrogen. Here, 

1/2 


c = 343.23 


S) ■ * 


in m/s 


(4.11) 


and (aA) maXii is a relaxation parameter described physically as the maximum ab- 
sorption, in Np, in a distance of one wavelength due to relaxation. Its standard 
expression is 

where 0; is the characteristic vibrational temperature and Xi is the fractional molar 
concentration. For dry air, these parameters have the values 


0 O = 2239.1 K 0 N = 3352.0 K 
X 0 = 0.209 X N = 0.781 


(4.13) 


Finally, / r> , is the relaxation frequency, in Hz. Expressions for this parameter are 


fr,0 

fr,N 



24 + 4.41 x 10 4 h 


0.05 + h 
0.391 + h 


9 + 350he _6142|(r/rorl/ ’' 11 


(4.14) 


In all of the above expressions, p , o — 101.325 kPa, T 0 — 293.15 K, and p,, T, /, 
and h are atmospheric pressure (in Pa), temperature (in K), acoustic frequency (in 
Hz), and molar concentration of water vapor (in %), respectively. Typically, the 
large variation in the relaxation frequencies due to variations in h account for the 
largest variations in a, while the direct effect of temperature on oc is much smaller. 
However, h is dependent on T and in this way temperature has an added effect. It 
so happens that variations in pressure can be handled through a scaling relationship 
involving / and h. The standard lists the accuracy of the above formulae for a at 
±5% for 2 < / < 15 kHz and 10% < h T < 100%, where h r is the relative humidity. 
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In the more extended frequency range of 4 x 10" 4 < //p. < 10 Hz Pa 1 , the 
standard lists the accuracy of a at < ±10%. As examples of the predicted losses 
due to absorption, when h r = 50%, T = 20° C and s = 100m, AIL, 4.13 dB 
at / = 5 kHz while AIL, = 14.69 dB at / = 10 kHz. These numbers indicate the 
nearly quadratic growth of AIL, with / for a given s. For the circumstances of the 
laboratory, the propagation distances were restricted to 4 m or less. 

In addition to the adjustment for absorption, a correction of 6 dB per dou- 
bling of propagation distance was made for spherical spreading. The losses due to 
absorption were assumed to be uncoupled from the losses due to spherical spread- 
ing. Further, the data was corrected for the frequency response of each microphone 
cartridge. Each of these frequency responses was measured by Bruel ic Kjaer techni- 
cians. In regard to the ambient conditions in the room, casual observations revealed 
that the they were relatively uniform and prone to only small variations with time. 
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CHAPTER V 


PRELIMINARY EXPERIMENTS AND ANALYSIS 


A set of preliminary experiments were performed in order to establish l) the nature 
of a typical incident waveform, and the effects of the finite size of the acoustic 
source and microphones, 2) the ambient speed of sound and the speed of sound in 
the plywood table top, 3) the directionality of the source, 4) the rate of propagation 
losses due to spreading and absorption, and 5) the size of nonlinear effects. The 

O 

objective was to gain a qualitative and quantitative understanding of the source 
waveform, of the free field propagation effects, and the effects of the microphone 
on the recorded shape of the waveform. Although all of these effects distorted the 
measured waveforms from the actual waveforms, the concern was with the effect 
of these distortions upon the measured insertion loss. In this chapter, the term 
“insertion loss” retains its usual definition but is used somewhat loosely in that it 
means the magnitude of the discrete Fourier transform of any field waveform divided 
by that of its associated reference waveform, in dB. Even if no object was “inserted” 
between the reference and field microphones, the above ratio, in dB, was referred to 
as an insertion loss. For example, if a single waveform was recorded at two different 
free field locations, then one could still calculate an insertion loss between the two 
locations, even though a zero insertion loss might have been predicted. 
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As is discussed and shown below, the measured incident waveforms contained 
a certain anomalous feature not found in a standard N-wave or short duration 
sine pulse. It was unclear whether this feature was due to the spark source or 
the microphone or both. With regard to the source, the relationship between the 
shape of the incident waveform and the discharge voltage, or the gap separation 
distance, or the orientation of the spark gap was unclear. As for the microphones, 
it was unclear what effect their orientation had on the recorded waveforms and, 
more importantly, on the measured insertion loss. 

5.1 Incident Waveforms 

Figs.(5.1-5.2) show typical waveforms and their amplitude spectra (the magnitude 
of the discrete Fourier transform), at a distance of approximately lm from the 
source, which resulted from spark discharges at voltage levels of 1.5, 2, 2.5, 3, 
3.5 and 4kV. In the time waveforms, two trailing troughs were usually found to 
follow the initial crest. Apparently, the size of the second of these troughs was 
not a consistent proportion of the size of either the first trough or the initial crest. 
Further, the peak pressure recorded for a 4 kV discharge was approximately 67% 
larger than the peak pressure recorded for a 2 kV discharge. All of the peak pressures 
were well within the dynamic range of the microphone. Further, if one accepts the 
notion that the diffraction effects due to the finite size of the microphone were 
proportional to the peak incident pressure (i.e. the shape of a waveform would not 
be significantly altered by these diffraction effects), then the shape of a waveform 
was apparently also a function of the discharge voltage level. The actual physics 
involved in a spark discharge was no doubt complex, and the resulting waveform was 
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probably a function of gap separation distance, and electrode material, geometry, 
and alignment as well. 

The concern was to pick a suitable discharge voltage given the spark source 
design parameters discussed in chapter 4. That is, it was desired to produce a 
waveform which was loud enough to keep the signal-to-noise ratio large at large 
R, but which was not so loud as to introduce nonlinear propagation effects, and 
which was of such a shape as to produce a relatively smooth and broad band fre- 
quency spectrum. An examination of the discrete Fourier transform of waveforms 
measured at various distances from the source showed that the third criterion was 
well satisfied by the waveforms produced by the 2kV or 2.5 kV discharges. Further 
tests indicated that sufficiently loud signals could be produced with these discharge 
voltages without unwanted nonlinear effects. In fact, a gap voltage of 2.2 kV was 
used in most of the experiments. (Interestingly, a gap voltage of close to 4.5 kV 
was necessary to initiate an untriggered discharge.) In terms of available electrical 
energy, 2.2 kV corresponded to approximately 2.4 J. Typical peak acoustic pressures 
at a distance of 1 m were on the order of 130 Pascals. A typical rise time and pulse 
duration was on the order of 10 /zs and 40 /zs, respectively. The resulting sound 
source was broadband with the peak level residing near 25 kHz. 

5.2 Effects of the Finite Source Size 

Consideration can be given to the possibility that the finite size of the source con- 
tributed to the anomalous shape of the waveform. Almgren performed experiments 
with a similar spark source and analyzed the effect of his finite sized source. He 
assumed that his spark source behaved in a fashion similar to that of short line 
source consisting of a continuous smeax of identical point sources all sounding in 
unison. For such a distributed source, it was possible to approximate the sound 
pressure at a point by taking the integral over the source volume of the Green’s 
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Figure 5.1 Typical measured waveforms produced by the spark source 
when the gap voltage was 1.5, 2, 2.5, 3, 3.5, and 4kV, respectively. The 
measurements were made with normal microphones with the protective 
grid in place. The duration of the initial rise and fall cycle was ap- 
proximately. 40 /xs in all cases. The peak acoustic pressure increased as 
the voltage was increased. The waveform at 2 kV had a peak acoustic 
pressure of approximately 130 Pa at 1 meter, and had no noticable trail- 
ing ripples. Most of the subsequent experiments involved waveforms 
generated with 2.2 kV across the spark gap. 
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Figure 5.2 Spectra produced from the waveforms shown in Fig. (5.1). Here, 
a spectra was defined as the magnitude of the discrete Fourier transform 
of a waveform. For a gap voltage of 2 kV, the peak spectrum magnitude 
was roughly 210 x 10“ 5 Pa/Hz and occurred at a frequency of approxi- 
mately 25 kHz. 
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function multiplied by the source strength distribution. 31 Let the length of the 
line source (e.g. the length of the air gap between the electrodes) be denoted by d B , 
and let the distance to a listener be denoted by R. Almgren reported that when the 
line source was oriented such that its axis was perpendicular to the line connecting 
the source and listener (as in the experimental arrangement described in chapter 
4), then the performance of the line source relative to that of a point source, in free 
field conditions, can be represented by 

jP,ho,t lin - e - = - [C{x) + iS(x)] (5.1) 

-Ppoint iourcc 



/ = 40 kHz, 5.3) yields 

•Pihort lin e 
> 

pointtource 

Thus, it appeared safe to consider the spark source as a point source, and it also 
appeared unlikely that the finite size of the spark source was responsible for the 
anomalous shape of the waveform. In addition, the thin (60 mil diameter, or « 
1.5 mm) pointed electrodes were not expected to contribute any significant reflected 



or diffracted rays to the incident field. 

Wright noted that deviations in the spark path from a straight line could cause 
distortions in the incident waveform. 32 Such flutter in the spark was difficult to 
quantify but given the dimensions of the spark source is seemed unlikely that this 
effect was significant at the frequencies of interest. No further attempt to analyze 


this effect was made. 



72 




Figure 5.3 Typical measured waveforms at 1 meter produced by the spark 
source when the gap voltage was 2.2 kV and the microphones were ori- 
ented normally and transversely, with and without the protective grids. 
Clearly, the overall sensitivity of a normal microphone is significantly 
higher than that of a transverse microphone. Also, the peak pressure 
measured with a normal microphone was greater when the protective 
grid was in place than it was when the grid had been removed. The 
same effect was not seen in the transverse microphone measurements. 
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FREQUENCY (kHz) 


Figure 5.4 Spectra produced from the waveforms shown in Fig.(5.3). 
Again, a spectra was defined as the magnitude of the discrete Fourier 
transform of a waveform. The spectrum of the waveform measured nor- 
mally and with a grid in place was significantly larger and broader than 
any of the others. The spectrum of the waveform measured transversely 
was the narrowest with only modest levels present past 25 kHz. The 
spectrum of the waveform measured transversely with the grid did not 
exhibit the same fall off at 25 kHz. 
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5.3 Effects of the Finite Microphone Size 

The fact that the microphones were not ideal omnidirectional, point receivers intro- 
duced distortions in the measured incident waveforms (Fig. 5.3— 5.4). Specifically, 
the sensitivity of the microphones varied with orientation and frequency. The cause 
of the distortions was attributed to the way in which a microphone averaged the 
pressure over its diaphragm face, and to diffraction effects emanating from the 
circular edge of the diaphragm. 

5.3.1 Microphone Orientation 

There appeared to be two sensible ways in which to orient each microphone; as 
stated in chapter 4, they were the normal and transverse orientations. A clear 
advantage of the normal orientation was the increased sensitivity due to near pres- 
sure doubling at the diaphragm face. This was important when the received signal 
was weak, particularly when the signal-to-noise ratio was not too large such as was 
found in the shadow zone. However, the sensitivity of a normally oriented micro- 
phone was also highly directional relative to a transversely oriented microphone. 
That is, the sensitivity of a normally oriented microphone was a function of the 
angle of incidence of a ray located in the centerline plane. (Due to circular sym- 
metry, the same was not true for a transversely oriented microphone.) Thus, errors 
in the normal measurements were present when a microphone measured direct and 
reflected or diffracted rays which arrived at varied angles of incidence. In order to 
reduce these errors, an effort was made to approximate the direct, reflected, and 
diffracted ray paths prior to a measurement and subsequently adjust the orientation 
of the microphone such that all rays arrived at the microphone at similar angles of 
incidence. The dimensions involved in most of the experiments were such that these 
angles of incidence were small, on the order of 15° or less. Thus, it was concluded 
that differences between the angles of incidence were even smaller. As measured 
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by Bruel & Kjaer, the nominal variation in sensitivity between 0° and 30° angles 
of incidence was not more than 2dB for frequencies less than 60 kHz. Therefore, 
the errors introduced by splitting the relevant arrival angles should have been much 
smaller than 2 dB. In this way, the directional effects of the microphones for normal 
orientation were minimized. 


5.3.2 Pressure Averaging over the Microphone Diaphragm 


Another concern was the effect of the finite surface area of a diaphragm and the 
fact that the recorded pressure was, in some sense, an average of the pressure at 
each point of the surface. Almgren considered a transversely mounted microphone 
and neglected diffraction effects from the edges and from the protective grid. In his 
case, and under the assumption that the microphone diameter, D , was much smaller 
than the distance from the source to the microphone, the relationship between the 
averaged pressure and the midpoint pressure is known to be 


5- flr. membrane average 

r — p 

•* ff , p oint 


2 7T 

sF 


f*(“) 


(5.5) 


where S is the membrane surface area, and J x is the first order Bessel function 
which has its first non-trivial zero at kD = 7.664. For D = 6.35 mm, this zero 
corresponds to a frequency of / = 66.4 kHz. Table 5.1 outlines the relationship 
between the ratio of the averaged pressure to the point pressure and the frequency. 

The reduction in P with frequency is attributable to the variable phase of the 
incident wave at different points of the diaphragm. The Bessel function, Jx{x), 
results from the circular geometry. From this analysis, it was clear that a trans- 
versely mounted microphone would tend to suffer significant changes in sensitivity 
at frequencies higher than 20 kHz. 

The sensitivity of a normally mounted microphone was also a function of fre- 
quency due to the same mechanism. However, when D <Zi R> the effect is small 
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Frequency 

(kH*) 

P 

1 

1.000 

10 

0.959 

20 

0.842 

40 

0.466 

60 

0.088 

66.4 

0. 


Table 5.2 The effect of the circular geometry and finite size of the 
microphone diaphragm upon the sensitivity of a transversely mounted mi- 
crophone. P is the ratio of the averaged pressure on the diaphragm to the 
pressure at the diaphragm center point. 


except at very high frequencies. This is true because the pathlength difference be- 
tween an incident ray which strikes the microphone center and that which strikes the 
microphone edge is very small. Thus, it seemed likely that a normally mounted mi- 
crophone would, in this sense, act more like a point receiver than would a transverely 
mounted microphone. 

Further, the diffusive effect of the protective grid would tend to reduce the 
loss of sensitivity with increasing frequency due to pressure averaging. Also, with 
regard to the insertion loss, it was thought that the pressure averaging effects at the 
field and reference microphones would tend to cancel each other out when the ratio 
of their discrete Fourier transforms was calculated. In sum, the pressure averaging 
effect due to the finite surface area of a microphone diaphragm should have had 
little, if any, influence upon the measured insertion loss, unless the frequency was 
such that P was near zero. Then, errors associated with the ratios of small numbers 
could have appeared. 
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5.3.3 Diffraction at the Edge of the Diaphragm 

Another effect which contributed to the anomalous shape of some of the wave- 
forms in Fig.(4. 1-4.2) was the diffraction of the incident wave at the edge of the 
microphone’s circular housing. Using the concepts of the Geometrical Theory of 
Diffraction, diffracted rays travel from the edge of the microphone housing inward 
over the diaphragm. Previous experiments using short pubes and baffled micro- 
phones have shown the existence of these waves. 33 Given the dimensions of the 
microphones used here (microphone diameter of approximately 6.35 mm), and given 
a sound speed of 343 m/s, the diffracted rays from the edge of the circular housing 
would arrive at the center of the microphone at a retarded time of approximately 
9.3 (is. This retarded arrival time was much smaller than the duration of the total 
waveform. The arrival of the diffracted wave upon the active portion of the micro- 
phone diaphragm could be delayed by constructing a relatively large baffle around 
the microphone. Then, the diffracted wave would arrive after and separate from 
the incident wave and, thus, could be gated out, as was done for other spurious 
reflections. However, the presence of a large baffle around each microphone would 
introduce other difficulties, and although the diffracted wave clearly affected the 
recorded shape of the incident waveform, it was hypothesized that the effect of the 
diffracted wave would cancel out when calculating the insertion loss. 

To test this hypothesis, it was decided to measure the insertion loss with nor- 
mally and transversely mounted microphones located at typical points of interest. 
The diffraction effects for the normal and transverse mountings should have been 
quite different. Therefore, if the diffracted wave (or the diaphragm averaging effect, 
for that matter) influenced the insertion loss to any appreciable degree, then, it was 
reasoned, the measured insertion losses for the two microphone mountings should 
be different. Measurements were made with normal and transverse microphones, 
with and without the protective grids, at the ridge apex, at a point 120 cm behind 
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the apex and on the line of sight, and at a point 8 cm above the second point. These 
three different field points were typical of the measurements described in chapter 
8. Although the calibration factors pertinent to normally mounted microphones 
should have been approximately correct for transversely mounted microphones, a 
separate calibration of the transverse microphones was made by placing them side 
by side equidistant from the source and comparing the RMS voltage values of each 
microphone’s measured waveform. The measured shapes were practically identical 
and the scaling ratio between microphones turned out to be roughly the same as 
that deduced for the normal microphones using a Bruel & Kjaer pistonphone. The 
results of the measurements were adjusted by this small correction. The measured 
insertion losses lay within a band of approximately ±ldB. As will be shown, this 
number is within the experimental error range deduced for the normally mounted 
microphones. The result was not surprising since the magnitude of the diffracted 
wave should have scaled with the magnitude of the incident wave, and, given that 
the microphones were similarly aligned, the scaled effect should have been the same 
at both the field and reference microphones, thus, yielding no effect upon the in- 
sertion loss. Therefore, it was concluded that the diffracted edge waves played an 
insignificant role in the value of the measured insertion loss. 

In sum, the results indicated that a normal microphone had the greatest average 
sensitivity and that with minor precautions presented only small errors due to 
orientation. 

5.4 Source Directivity 

A set of waveform measurements were made at seven different field points all lo- 
cated approximately 1 m from the source. Each field point was separated from its 
neighboring points by an angular distance of 6° such that the total set of points 
spanned an angular range of -12° to 24° from the horizontal. The results are shown 
in Fig.(5.5). Over a broad band of frequencies and over the extent of the angular 
range, the ratio of the discrete Fourier transform of the field waveform to that of 
the reference waveform deviated from unity by less than 1 dB. 
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Figure 5.5 Insertion loss measured at seven distinct field points all located 
approximately 1 m from the source but at different angles of inclination 
from a horizontal line which passes through the spark gap. Results over 
a broad frequency range are shown. 




TRAVEL DISTANCE (m) 


Figure 5.6 Arrival times of direct pulses measured at ten different distances 
from the spark gap. The dashed line represents a least squares best fit 
of the data. The inverse slope of this line gave a sound speed in air of 
345.6m/s. 
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5.5 Speed of Sound 

Free field waveforms were recorded at 10 different field points, all of which were 
located approximately 40 cm above the table top and which spanned the length of 
the table. The field points were located at 25 cm intervals from R « 175 cm to 
R « 400 cm, where R was the distance from source to microphone. A reference 
waveform was also recorded for each of the 10 field waveforms. The sound speed 
was estimated by plotting the arrival time of each waveform versus the its associated 
value of R. The plotted data and a least squares best fit line are shown in Fig. ( 5.6). 
The slope of the best fit line gave a sound speed estimate of 345.6 m/s. This estimate 
was very close to the sound speed predicted by 

c.et = (1 + 0.16/i)c dry = (1 + 0.16h) (331.3 + 0.61T C ) = 346.5 m/s 

when T e = 22.2° C and h « 0.03, where h is the fraction of molecules in air that 
are water. The given value of h corresponds roughly to a relative humidity of 50%. 
The calculation was made assuming p A i r = 1.20 kg/m 3 and P* m b = 740mm Hg. 

The speed of sound in the plywood table top was estimated by a similar tech- 
nique. Two high frequency (1 MHz) transducers were mounted on both sides of 
a small patch of plywood identical to that comprising the table top. A suitable 
couplant was smeared between each transducer and the plywood in order to ensure 
transmission. Two cycles of a 40 kHz sine wave were amplified and used to drive one 
of the transducers while the other was used as a receiver. Transmitted signals were 
recorded for a single and double thickness of the plywood. The data indicated that 
the second thickness of plywood resulted in a delay of roughly 18 n s. The plywood 
thickness was 3/4 in (« 1.9 cm). Thus, a sound speed estimate of 1050 m/s was 
obtained. A substantial amount of unwanted noise, an anomalous signal compo- 
nent in the data, and the high attenuation in the plywood made the 18 ps reading 
somewhat questionable. Therefore, a second attempt at measuring the sound speed 
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in the plywood, was made by measuring the transmission time across a 1 in thick 
piece of plexiglass and comparing this measurement to one made across the plexi- 
glass and plywood in series. This attempt resulted in an estimate of 1300 m/s for 
the sound speed in the plywood. Both of the sound speed estimates seemed rea- 
sonable for the plywood since the speed of sound in many woods is on the order of 
3000-4000 m/s with the grain, although the speed can be substantially lower across 
the grain. That the plywood sound speed should be lower than a typical wood is 
supported by the notion that the glued layers comprising the plywood sheet contain 
a significant amount of trapped air. Also, to the extent that plywood has a grain 
pattern, the measurements were certainly made across the grain. As an example of 
the effect of trapped air on the sound speed, consider that the listed sound speed in 
cork is 500 m/s, which is substantially lower than the sound speed in most woods. 
However, all of the above figures are significantly higher than the sound speed in 
air. Thus, it seemed reasonable and safe to conclude that the sound speed in the 
plywood is greater than the sound speed in air, regardless of inhomogeneities or 
anisotropy. 

5.6 Free Field Propagation Losses 

The same set of 10 field measurements were used to deduce the rate of free field prop- 
agation loss. The measured free field insertion loss, for given discrete frequencies, 
is shown in Fig. (5.7). The data was normalized with respect to spherical spreading 
losses of 6dB per distance doubling, and with respect to absorption losses in air 
as predicted by the ANSI standard. All of the data lie within 2 dB of zero, with 
the lower frequency data being within 1 dB of zero and being consistently closer 
to zero than the higher frequency data. This is in keeping with the accuracy of 
the standard which is expected to decrease with increasing frequency; in fact, the 
standard is expected to underpredict the actual losses at high frequencies. 
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Figure 5.7 Insertion loss measured at ten different free field locations. The 
data was adjusted for spherical spreading and absorption, as prescribed 
by the ANSI standard. ,-s 
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5.7 Nonlinear Propagation Effects 


The short duration, finite amplitude pulse used in these experiments lends itself to 
the existence of nonlinear wave steepening. The concern was with the extent of such 


effects which occurred between the position of the reference microphone and that 
of the field microphone for any particular measurement. To estimate the extent 
of these nonlinear effects, consider the propagation of a simple, finite amplitude 
wave. 34 If one assu mes that the particle velocity is a single-valued function of 
pressure, and that 


dp 

dp. 


= c 


then the mass and momentum conservation equations yield 


(5.6) 


S+<« + .)£-o 


(5.7) 


where u is the particle velocity, c is the isentropic speed of sound, and where both 
are functions of p. As is well known, this equation indicates that points of a given 
pressure move with the s am e velocity and that points of different pressure move with 
different speeds. Further, if (c + u) a p, then points of greater pressure move faster 
than points of lesser pressure, and wave steepening occurs. (5.7) can be reexpressed 
as 

g + fc+MfS-o M) 

where 0 = (i + l)/2 is the coefficient of nonlinearity for an ideal gas, and c 0 is 
the ambient sound speed. Now, if p(0,t) = g(t), then the approximate Earnshaw 
solution to (5.8) is 

p(x,t)=9W (5 ' 9) 

Letting g(rjj) = p 0 sinu >xj>, the shock formation distance, f, is defined as the earliest 
value of x for which g(ip) is a multi-valued function of t — x/c 0 . This distance is 
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-.3 
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(5.10) 
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where (i — 1.2 and. p « 1.2 kg/m 3 in air. Letting c 0 = 345.6 m/s and p — 120 Pa, 
then 

£ tv 5.5 m at / = 10 kHz (5.11) 

Of course, £ oc 1 //, so the shock formation distance is smaller for higher frequencies. 
The value of 5.5 m is on the order of the dimensions of the laboratory; however, it 
should also be a very conservative value since it is based upon a waveform which 
propagates without dissipation or spreading. The spark waveform is, in fact, spher- 
ical not planar. Thus, the peak pressure is attenuated at least by 6dB per distance 
doubling. This fact extends the shock formation distance although it is well known 
that spherical spreading cannot on its own prevent the formation of a shock. An 
estimate of the shock formation distance of a spherical wave in a lossless medium 


can be had from the formula 


^oniet *0 GXp 


/>c 3 //?r 0 


. . (5.12) 

l(dp(r„t)/dt) m . x J 

An estimate of the denominator can be taken from the plot of a typical waveform 

at r 0 = lm. From Fig. (3.3) or (4.1), a reasonable value seemed to be 

/ dp(r 0 . f ) \ 120 Pa 
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dt /max 


8/iS 


(5.13) 


Given this value, r oniCt ^ 15.7 m, nearly three times as far as the plane (simple) 
wave value. As mentioned above, the key concern was the extent of nonlinear 
effects between distances of 1 m and 4 m. The above analysis indicated that they 
were probably small but not necessarily negligible. An attempt was made to observe 
nonlinear effects by comparing the discrete Fourier transforms of free field waveforms 


captured at 1, 2, 3, and 4 m. After corrections were made for spherical spreading and 
absorption losses, no clear indication of an energy redistribution from low to high 
frequencies was apparent. The data shown in Fig.(5.7) was some of the data used 
to make the observation. In addition, the shapes of the waveforms were observed 
for signs of wave steeepening. Again, no significant changes of shape were apparent. 
For these reasons, it was concluded that nonlinear effects play a negligible role in 
the waveform propagation induced by the spark source. Consequently, the ambient 
speed of sound was denoted by c throughout this thesis. 



CHAPTER VI 


SPECIFIC ACOUSTIC IMPEDANCE 

OF THE 

CURVED SURFACE 


Recall from chapter 3 that the Fock-van-der-Pol-Brennner function, G, was a func- 
tion of the scaled parabolic coordinates, (£,*/), and of the scaled admittance pa- 
rameter, < 7 , which characterized the boundary condition at the surface. A value 
of the parameter q was required before G could be computed and compared with 
experimental results. In this chapter, a method used to evaluate q for the two lab- 
oratory scale surfaces (chapter 4) is described. As will be shown in this and the 
following two chapters, the results of this method yielded excellent comparisons to 
experimental data. The results also tend to indicate that the method should work 
reasonably well for a wide range of real surface types. 

The task of evaluating q is equivalent to that of evaluating the specific acoustic 

impedance of each surface since 



where R is the radius of the curved surface, k and pc are the wavenumber and char- 
acteristic impedance of air, respectively, and Z, is the specific acoustic impedance 
of the curved surface. The specific acoustic impedance is defined by 
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where So is the nominal location of an infinitesimally small area on the surface. 
A direct measurement of the acoustic impedance of a surface is not a simple task. 
The difficulty becomes apparent when one considers that, to obtain Z, directly, the 
pressure and normal velocity of the fluid at the surface would need to be measured 
simultaneously. To complicate matters, Z, is not an intrinsic property of the surface 
but depends upon the incident wave as well. In particular, Z, = Z,[0 lt u) where u 
is the frequency of the incident wave and is the angle between the propagation 
direction and the inward normal of the surface. To simplify matters, a number of 
assumptions concerning the surface and the incident wave were made so that Z, 
could be deduced from measurements of the pressure, p, at points off of the surface. f 

6.1 Analytic Methods of Evaluating Surface Impedance, Z 3 

A reasonable and simple first approach was to assume the case of a plane wave 
incident upon locally reacting surface. A locally reacting surface is one on which the 
motion of any particular portion of the surface is independent of the motion of the 
rest of the surface. As shown below, an implication of the local reaction assumption 
is that the surface impedance is independent of the angle of incidence, 6 x . These 
assumptions result in a well known and simple expression for the pressure at any 
point in the field. However, there are limitations associated with these assumptions, 
and in certain cases, particularly when 0^, « tt/2, it is necessary instead to consider 
the reflection of a spherical wave from a plane boundary. 

Consider two semi-infinite homogeneous half spaces separated by a plane im- 
pedance boundary. Let a plane wave of constant angular frequency u> be located in 
one of the half spaces and be incident upon the boundary at angle . (The half 

j* No attempt was made to use an impedance tube to make the required meaiurement 
because kR was made large (intentionally). At high frequencies, standard impedance 
tube techniques are inapplicable. 
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space containing the incident wave is referred to as the upper half space and its 
properties are denoted by a subscripted numeral 1. Its image space is referred to 
as the lower half space and its properties are denoted by a subscripted numeral 2.) 
Then, the acoustic pressure at the receiver can be written as 


p = pe 


— tw t 


(6.3) 


with 

p = Pi (6-4) 

where P< is the magnitude of the incident plane wave, k x is the wavenumber vector, 
F/ and f R are position vectors associated with the incident and reflected waves, and 
J? p (fl 1 ,w) is the plane wave pressure-amplitude reflection coefficient. This reflection 
coefficient can be expressed as 


R p 


Z 2 cos 9i — Pi Ci 
Z 2 cos &i + piCi 


(6.5) 


where Z 2 is the specific acoustic impedance of the lower half space, and p x c x is the 
characteristic impedance of the upper half space. If one lets Z 2 = P?c 2 1 cos 0 2 , then 
Z p can be reexpressed as 


A special version of (6.6) 
(Snel’s) Law 


£ _ P2C2 cosfli - P1C1 cos #2 6 ) 

p p 2 c 7 cos Si -I- pi Ci cos 0 2 

is found in cases where k x 2 < k 2 2 . In such cases, Snell’s 


ki sin 0i = k 2 sin 0 2 


(6.7) 


indicates that 0 2 is close to zero, independent of the value of 9 X . Then, (6.6) reduces 


Pj C<2 COS 0 1 p\ Cj 

p 2 C 2 COS 0! + PiC x 


to 


( 6 . 8 ) 
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Thus, the reflection coefficient is still a function of Oi, but Z 7 « pic 7 Is a constant 
impedance independent of 0 X . In such cases, the boundary is referred to as locally 
reacting; that is, when the propagation in the lower half space is negligible due to 
a slow sound speed, Ci ^ c 7 . Another case when the boundary can be considered 
as locally reacting is when the attenuation in the lower half space is large enough 
that the transmitted wave can be neglected. 

In principle, one could measure the phase and amplitude of the pressure at 
any point above the surface, and use (6.4) to deduce R p , and subsequently the 
surface impedance, Z 7 . However, in practice, neither a true incident plane wave 
nor true plane wave reflection at the boundary is attainable. The second of these is 
particularly true at large angles of incidence. In fact, a plane wave cannot propagate 
parallel to a flat surface of finite impedance. 36 . 

One may account for the lack of a true incident plane wave by writing the 

pressure as 

p (j ( e — 

P<~\k .R, ) +Kp \k l R 2 
which is valid for a point source located such that Ri and R 7 are the distance 
from the source to the receiver and the distance from a mirror image source to 
the receiver, respectively. As kR gets large, this expression approaches (6.4) (i.e. 
the incident spherical wave becomes approximately planar when the point source 
is moved far away). Thus, in practice, a reasonable approximation for a plane 
incident wave can be had. Still assuming plane wave reflection at a locally reacting 
boundary, then the impedance can be deduced from (6.9) and (6.5). At frequencies 
such that the errors in the measurement of source and receiver positions are small 
compared to a wavelength, and at moderate angles of incidence, this method should 
yield reasonable results. However, if one assumes that Z 7 is finite, then as the angle 
of incidence 0 -*■ tt/ 2, (6.5) indicates that £ p -► -1, and thus, from (6.9), p -» 0 
at the boundary. This result is referred to as the grazing incidence paradox since 
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it predicts that no energy is absorbed by the surface regardless of the value of 
the surface impedance, although it seems clear that some of the energy of a wave 
propagating parallel to a plane surface would go into causing motion of the surface 
in a direction perpendicular to the assumed direction of the wave*” 27 f 

This contradiction can be avoided by considering spherical wave reflection from 
the plane boundary. Work in this regard goes at least as far back as Sommerfeld 
[1909,1926] who worked on problems of long range radio wave propagation. Other 
contributions were made by van der Pol [1935] and Norton [1936], who first outlined 
the form of the solution used currently. As in the case of a rigid boundary, they 
interpreted the solution in terms of a source point and its image below the plane 
boundary. However, in their interpretation, this image point lies at the center of 
a roughly elliptical region from which secondary waves are emitted. 28 The size 
of this region is an inverse function of the magnitude of the imaginary part of the 
wave number in the medium below the boundary. When this region is small, due to 
a sufficiently large attenuation in the lower medium, the combination of the image 
source and the secondary waves was referred to by van der Pol as a diffuse image 
source. Clearly, one characteristic of a diffuse image source is that its behavior 
approaches that of a classic image source as the boundary impedance approaches 
that of a perfectly rigid boundary. Rudnick [1947] followed the work of these men 
to present the analogous solution for acoustic propagation over an impedance plane. 

A complete and up-to-date treatment of the derivation and solution was pre- 
sented by Brekhovskikh [1980]. The derivation made use of the fact that a spherical 
wave can be expressed as an integral of homogeneous and inhomogeneous plane 

| Ingard noted that the gracing incidence paradox could be overcome by eimply eetting 
$ = p Q co # n $ = p iCi/z 3 . With iuch an expreeeion for the impedance, and when 
$ = t/ 2, the gracing incidence paradox is found when n — 0. But when n = 1, the 
reflection coefficient is (/9 0 ” O/fao + 0» * nd th ® surface pressure is non-sero. In fact, 
when n> 1, the reflection coefficient is 1, and the surface pressure is double the incident. 
See U. Ingard, Ph.D. dissertation, 1950. 
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waves. Subsequently, the method of steepest descent was used to derive an asymp- 
totic solution for large kR. The details of this procedure are rather involved; there- 
fore, the end result is simply stated here as 


£ 

A 

S 


gilt i R i 

ki Ri 


+ TIT K 1 - *p) F M + *.l 

K± ilj 


( 6 . 10 ) 


where 


F(to) = 1 + iy/irwe w erfc (— ty/w) 


( 6 . 11 ) 


to 


2k x Rz Z\ 


K-5-") 


( 6 . 12 ) 

(1 -Rp)*Zi'~ 

and where erfc(z) is the complex complementary error function and 0 is the angle 
of incidence. This result is sometimes referred to as the Weyl-van der Pol equation. 
The ray acoustic solution has reappeared on the right hand side of (6.12), along 


with the pl an e wave reflection coefficient, R p . The new term, F, was referred to 
by Rudnick as the boundary loss factor while others referred to it as the ground 
wave f un ction, to has been commonly referred to as the numerical distance .f The 
function F is such that as w — ► 0, F — * 1 and as w ► oo, F * 0. Examination of 
(6.12) reveals that these two cases represent the case of a rigid surface and the case 


of plane wave reflection, respectively, 
of to in that 

to a 


This behavior is consistent with the definition 


k^Rz Z\ 

(i - *,)’ 


(6.13) 


Thus, it is clear that i) when the source is very far from the receiver or the source 
is very far above the boundary, to ^ oo and the plane wave result is found, and ii) 
when Zz -* oo, to -*• 0 and the rigid surface result is found. (One should be careful 
to note that when Zz — y oo, R p — * 1 and to is indeterminate. However, in this 


The numerical dUtance is related to the distance (in the complex plane) between poles 
of the plane wave reflection coefficient and the saddle point found as part of the method 
of steepest descent used in the solution. 
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case (6.10) gives the correct rigid surface result directly after plugging in R p = 1 
and knowing that F{w) is bounded.) Finally, the aforementioned grazing incidence 
paradox, found in cases of plane wave reflection, is overcome since when 0 -*■ sr/ 2 
and — ► —1, 

e iklRl / \ 

— 2F[w) — — (6-14) 

ill 

In this form, it is clear that F is an attenuation function since for F = 1, (6.14) 
yields the rigid boundary solution. Interestingly, for certain values of the numerical 
distance, particularly when arg(tx>) — * tt/2, |.F| can be greater than unity. Then, 
surface pressures greater than those at a rigid surface are predicted. 

Similar work by Wenzel showed that for certain types of surface impedances, 
the solution for the total field contains a surface wave term. A surface wave is 
characterized by a magnitude which decreases exponentially with height above the 
impedance plane and by a propagation speed parallel to the boundary which is less 
than the free space propagation speed. Using the notation of Wenzel, this surface 
wave was expressed as 

^e~' T ' l Hj ) 1) [(fc 2 + 7 2 ) 1/2 rj where 7 = (6.15) 

where h is the distance between the source and its image, and r is the horizontal 
distance between the source and receiver. It is worth noting that the asymptotic 
behavior of the Hankel function is 



In this form, the magnitude of the surface wave is seen to fall off exponentially 
with distance along the surface. However, if the imaginary part of 7, which is 
proportional to the real part of Z 2 , is small enough, then the magnitude of the 
surface wave falls off as r -1 ^ 2 . Since the other terms in the total solution fall off 
as r _1 , the magnitude of the surface wave may dominate the solution. In such 
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cases, a so-called negative excess attenuation can occur near the surface. (Excess 
attenuation is defined as the attenuation beyond that due to spreading or pressure 
doubling at a surface.) Another aspect of the solution, which is apparent from 
(6.16), is that the phase speed of the surface wave is slower than that in free space 

when the imaginary part of is positive. 

A surface wave term is not always present in the solution. Rather, its existence 
is determined by the following relationship between Z and 9, where the subscripts 
have been omitted, and Z and 9 are understood to be Zi and 9 1 . 

Im |t'(l + 0cos9 - (1 - /? 2 ) l/2 sin0)j <0 (6-17) 

where 9 is the angle of incidence and (3 = P\Cx / Z . 20 When the surface wave does 
exist its contribution is typically weak enough such that a positive excess attenuation 
results in the field over the impedance plane. However, as discussed above, when 
the real part of the impedance is small and the imaginary part is relatively large, 
the surface wave contribution can be significant. It is the surface wave component 
which is responsible for those instances when \F\ > 1. 

6.2 Proposed Experimental Method 

As stated before, the interest was to obtain an estimate of the surface impedance at 
grazing singles of incidence so as to evaluate the scaled admittance parameter, q. It 
was assumed that the impedance at a point on a curved surface was equivalent to 
that of a point on a flat surface, given that the surface material properties and the 
angles of incidence were the same. This assumption seemed reasonable from a ray 
acoustics point of view and, since high frequency propagation was assumed, it was 
adopted. Then, using the above formulae, the impedance of a plane surface could 
be calculated or inferred from measurements of the pressure above that surface. A 
method for doing this was outlined by Chessel and successfully applied by Embleton 
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tl a/.. 30 The method consisted of comparing measured sound pressure levels with 
predicted sound pressure levels for the field of a point source above an impedance 
plane. Of course, to predict sound pressure levels using the above formulae, an 
estimate of the impedance for each given angle of incidence was needed. However, 
it was assumed a priori that the surface was locally reacting, with the intent of 
determining a posteriori the validity of the assumption from the results of its ap- 
plication. Thus, it was decided to use a semi-empirical, single parameter model 
put forth by Delany and Bazley for fibrous absorbent materials. 31 The model 
estimates the impedance independent of the angle of incidence, and so it inherently 
assumes local reaction at the surface whenever it is used for cases of oblique angles 
of incidence. The resulting impedance estimate was used in (6.10) to predict the 
SPL over a wide range of frequencies, and the resulting curve was compared to that 
found from experimental measurements. When a curve of the predicted SPL, for 
some estimated impedance, was found which matched the experimental data, then 
the actual impedance was assumed to be equivalent to the estimate used in the 
prediction. 

The single parameter model was obtained from an analysis of measured impe- 
dance data for a wide range of fibrous absorbent materials. All of the measured 
data was acquired through the use of an impedance tube and, thus, the data rep- 
resented normal incidence surface impedances. The model related the real and 
imaginary parts of the impedance, Z 2 , to the ratio of frequency divided by specific 
flow resistance per unit thickness (henceforth referred to as the flow resistance). 
Relationships to the same ratio were also found for the complex propagation co- 
efficient in the absorbing material. Expressed as power laws, these relationships 
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were 


Re /-^-l = -^_ = l + 0.057 (Px//a)- 0 - 75 

ypiCx) pic x 

(5.18a) 

Im{ — 1 = — = 0.086(pi / /a) - 073 
L Pi c i ) pi c i 

(5.186) 

Re {£} = = 1 + 0.096( ( . 1 //<t)-» ' 0 

(5.18c) 

Imjti) = Sl =0.1920., 

(5.18d) 


where / was the frequency, p 1 c l was the characteristic impedance of air, Z 2 was the 
specific acoustic impedance of the reflecting surface, k 2 was the complex wavenum- 
ber in the absorbent material or surface, and o was the flow resistance parameter in 
kg/(m 3 -s). The use of [pif/o) as the scaling parameter was suggested in the the- 
ory of sound propagation in a porous medium as presented by Zwikker and Kosten 
[1949]. 

Attenborough pointed out that a price paid for the single parameter model was 
that it fixed the functional dependence of the impedance to frequency as well as 
the relative magnitude of the real and imaginary parts of the impedance. These 
constraints could become problematic for uses of the model outside the range of 
/ and a associated with the experimental data upon which the model was based. 
Delany and Bazley stated that the model “may be used with confidence within 
the interpolating range 12 < Pl f/a < 1200 (where Pl = 1.2 kg/m 3 ) but caution 
against its use otherwise. In particular, they warned against its use for low values 
of Pl f f 0 \ that is, for acoustically hard surfaces. Indeed, Attenborough reported 
some difficulty when using this model to predict highway noise propagating from a 
line of traffic over a barren or lightly vegetated ground. He also showed that a five 
parameter model involving flow resistance, porosity, tortuosity, steady flow shape 
factor, and dynamic shape factor yielded “superior” results compared to those of 
the single parameter model when hard but porous surfaces such as sand or soil 
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were considered. In the five parameter model, Attenborough showed, for large flow 
resistivities and low frequencies, that Z 3 and Aj 3 were proportional to 

JL) (6.18) 

aCl ) 

where ft is the porosity, and that the real and imaginary parts of Z 3 were equal. t 
This result showed that i) the scaling exponent was no longer the same as that 
in the single parameter model, and ii) it was likely that what Delany and Bazley 
called the flow resistance, a, was, in fact, better termed the effective flow resistance 
since it probably included the influence of other factors such as the porosity, Cl. 
Attenborough also stated that “for high flow resistivities the empirical formulae give 
severe overestimates of attenuation constants.” Furthermore, at high frequencies 

Z * t ^2 ~ r (6.19) 

Pi ^1 ^ ^1 

where r was the tortuosity. However, in another study by Chessel, the single pa- 
rameter model found good use outside of its stated range of validity, and, due to its 
simplicity, was worth using as an approximate relation between a characteristic sur- 
face property, a, and the surface impedance, Z 2 . The five parameter model, while 
widely applicable and promising in accuracy, required the measurement or assess- 
ment of both the porosity and flow resistivity, and the acceptance of an approximate 
relationship between the dynamic and static shape factors. 


f Konstantinov found the same functional dependence between Z a and f when considering 
the effect of thermal and viscous boundary layers on reflection from a hard boundary. 
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6.3 Actual Experimental Procedures and Results 

Before the experimental method outlined above was used, an attempt was made to 
evaluate the impedance, Z 2 , by assuming specular reflection and directly applying 
(6.8). The experimental setup was such that the spark source and receiver were 
located high above the table surface (Fig. 6.1), such that the direct and reflected 
pulses arrived separately at the receiver. The angle of incidence was simply stated 
as 

tan (|-#)- (6-20) 

An advantage of this approach was that it required only one microphone, and so no 
account was needed of frequency responses amongst different microphones. In order 
to mitigate the effects of the directional response of the microphone, the microphone 
was oriented such that the direct and reflected rays impinged upon the diaphragm 
face at equal angles of incidence; that is, such that 2^ = tt/2 + (a — 6), where 


tan a = 



( 6 . 21 ) 


Another advantage of the single microphone capability was that the entire 500 kHz 
sampling rate of the A/D converter was available to collect data. This rate was 
sufficient to assure accurate measurement of the arrival time of the pulses to less 
than 5 % of the period, T, at 20 kHz. An obvious disadvantage of this method was 
that in order to keep the arrivals of the direct and reflected pulses distinct at large 
angles of incidence, the distance L had to be large enough so that A Rf c > T. This 
requirement tended to decrease the signal-to-noise ratio since, in practice, the pulses 
were spreading spherically. On the other hand, the maximum size of L was dictated 
by the dimensions of the laboratory. This limitation on L y along with the need for 
separate arrivals of the pulses, prohibited measurements at angles of incidence less 


than 10° . 
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Source 



Figure 6.1 Sketch illustrating the method used to acquire initial estimates 
of the surface impedance. Depicted are the direct (Ri) and reflected 
(R 2 ) ray paths between a point source and a listener. The reflection 
occurs according to the law of mirrors but with the reflected ray ex- 
periencing a change of magnitude and phase during the reflection. In 
the actual experiments, a spark source produced a short transient. Be- 
cause {R 2 - Ri) was sufficiently great, the two arrivals were separated 
in time. With appropriate corrections for spherical spreading and ab- 
sorption, the incident and reflected waveforms at the point of reflection 
could be deduced. Then, the ratio of their discrete Fourier transforms 
yielded the reflection coefficient, and with additional calculation, the 
surface impedance. 
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The incident and reflected waveforms were corrected for spherical spreading 
and absorption, and the reflection coefficient was calculated as the ratio of the 
magnitudes of the discrete Fourier transforms of these waveforms. Results obtained 
at angles of incidence of 10° -60° for the bare plywood surface indicated that the real 
part of the plywood impedance was relatively small and positive but independent 
of frequency. However, the imaginary part of the plywood impedance decreased 
roughly as 1 j f but was significantly larger than the real part over the range of 
frequencies 5-25 kHz. Near 25-30 kHz, the imaginary part tended to go through 
an inflection and change sign. Thereafter, its magnitude tended to grow at least 
quadratically. Thus, the magnitude of the impedance tended to go through a mini- 
mum near 25 kHz. Finally, the magnitude of the imaginary part tended to increase 
as the angle of incidence increased. Thus, these results indicated that the surface 
was not locally reacting. 

A troubling aspect of this data was that a relatively large scatter was found 
in the measured reflection coefficient phase data while the magnitude data was 
repeatable. At high frequencies, this result was not surprising since the error in 
the measured positions of source and receiver scaled inversely with the wavelength. 
This error was believed to be the source of most of the trouble with this technique, 
certainly at moderate angles of incidence. As discussed above, for large angles of 
incidence, the analytical model itself was questionable. 

Measurements such as those described above were made over the carpet-on- 
plywood surface as well. The results at moderate angles of incidence, 10 6 -60°, 
indicated that the real and imaginary parts of the impedance were approximately 
the same with both falling off roughly as 1 //. Again, a large scatter was found 
in the phase data. Also, the finite thickness of the carpet and the presence of the 
hemp and glue carpet backing made an estimation of the effective point of reflection 


difficult. 
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When it was not possible to acquire distinct arrivals of the direct and reflected 
pulses, the impedance was calculated using a method relying upon the standing 
wave ratio (SWR). In this method, the direct and reflected pulses were transformed 
together so as to produce an oscillating interference pattern. The magnitude and 
phase of the surface impedance can be deduced from the magnitudes and locations 
of adjacent nodes and antinodes, respectively. 33 However, this method also relied 
upon the accurate measurement of the source and receiver positions, and so difficulty 
was encountered in obtaining repeatable estimates of both the magnitude and phase. 

The method finally settled upon was that developed by Chessel, used by Em- 
bleton, et a/., and described in the previous section. As mentioned above, Ches- 
sers method considered the reflection of a spherical wave from a plane impedance 
boundary. 33 This consideration was important when the angle of incidence was 
large, as was true for the model experiments conducted in this study. Also, use of 
the s emi -empirical, single parameter model of Delany and Bazley allowed for the 
calculation of the impedance, Z 2 , from the magnitude of the pressure at a point 
above an impedance plane; that is, no phase information was required. Thus, the 
difficulty of accurately and precisely measuring the source and receiver positions 
was relieved. 

Typical results obtained for the carpet-on-plywood surface at a number of 
different large angles of incidence are shown in Figs. (6.2- 6.3). The data shown in 
Fig.(5.2) is that for inclination angles of 1.1° , 5.0° , 9.8° , and 15.3° from the horizon. 
Fig.(5.3) shows data obtained at an angle of incidence of 2.2° from the horizon. 
Other inclination angles examined were 1.3°, 1.7°, 3.6°, and 7.0°. The height 
of the receiver was 0.3 cm in all cases. The horizontal distance ranged between 
60 cm and 120 cm, while the source height was adjusted to provide the desired 
angle of inclination. The solid line represents the results predicted by (6.9) when 
the impedance is estimated by the single parameter model with an effective flow 
resistance of o — 1600 x 10 3 kg/(m 3, s). 
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Figure 6.2 Typical results obtained for the determination of the effective 
flow resistance, a , of the carpet-on-plywood surface. Shown are the 
measured SPL (0) and the predicted SPL (— ) versus frequency curves 
for sound propagating over the flat table top covered with a thin com- 
mercial carpet. In this context, the SPL refers to the sound level, in dB, 
at the point of reflection relative to that of the free field extrapolated 
to that point. The predicted SPL curves were computed using a value 
for the effective flow resistance of o = 1.6 x 10® kg/(m 3, s). The listener 
height was 0.3 cm while the horizontal distance between the listener and 
source was either 100 or 120 cm. The source heights were such that the 
shown angles of incidence were produced. 
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Figure 6.3 Typical result obtained for the determination of a variable ef- 
fective flow resistance, a, for the ca rpet-on-plywood surface. The curves 
and data represent SPL quantities similar to those shown in Fig. (5.2). 
The solid curve represents the predicted SPL curve computed using a 
value of o - 1.6 x 10® kg/(m 3 -s) for the effective flow resistance. The 
dashed curve represents the predicted SPL computed using an effective 
flow resistance which varied linearly from a — 1.6 x 10® kg/ (m 3 *s) at fre- 
quencies below / = 14.5 kHz to a = 0.8 x 10® kg/(m 3 *s) at frequencies 
above / = 17 kHz. The listener height was 0.3 cm, the source height was 
2.0 cm, and the horizontal distance between the listener and source was 


60 cm. 
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The experimental data fit the prediction very well up to approximately 14 kHz, 
but it overestimated the experimental data for frequencies greater than 14 kHz. 
As shown in Fig.(5.3), a improved fit was found when a was varied linearly from 
1600 x 10 3 to 800 x 10 3 kg/m 3 s over the frequency range 14 < / < 17.5 kHz, and then 
maintained at a = 800 x 10 3 kg/m 3 s past 17.5 kHz. It was expected that the best fit 
value of a would change with frequency because the relative importance of porosity 
and other factors, which do not scale the same as flow resistivity, was expected to 
change with frequency. In fact, as noted previously, porosity and tortuosity should 
have become increasingly important at higher frequencies. Furthermore, an effective 
flow resistance parameter might be expected to decrease with increased frequency 
simply on grounds that as frequency increased, the relevant length scales decreased, 

and the effective size of a pore increased. 

Somewhat surprising was that a good match between theory and experiment 
occurred at all. Embleton et a 1. pointed out that the formulae of Delany and Bazley 
are invalid for surfaces consisting of a thin porous layer on top of a hard backing. 
In fact, Delany and Bazley advised adjusting the impedance for such cases by the 
formula 

2 7 = Z 2 coth (6.22) 

where Z 2 = Ri + iX 2 , 7 = - ia 2 is the propagation constant with & 2 and a 2 

as defined in (6.17), and £ is the effective thickness of the layer. As mentioned 
above, and in chapter 4, the carpet-on-plywood surface consisted of a 1/4 in thick 
commercial carpet laid atop a plywood table. The carpet fabric was not tightly 
woven. However, the carpet fabric was attached to a backing of glue interspersed 
amongst a woven grid of thin hemp cord. The glue layer was riddled with holes 
on the order of 1 mm in size. This glue layer may have acted as a soft absorber or 
scatterer atop the relatively hard plywood surface. Thus, it was unclear whether 
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the above adjustment for a rigid backing was appropriate for the carpet-on-plywood 
surface. 

Many trial calculations using the correction (6.22) were made. The possible 
values of £ were alternately chosen to be the total thickness of the carpet, one-half 
the total thickness of the carpet, and the thickness of the glue-hemp backing. The 
reasoning was that, if the carpet behaved as a single layer at all, then the effective 
layer thickness might be the total thickness, or the mean thickness, or perhaps, if 
the carpet fabric was assumed to play no role in reflection or absorption, then the 
thickness of the glue-hemp backing. The various numerical results generated using 
(6.22) provided no better fit with experiment than that generated without (6.22). 
In fact, in most cases the fit was worse. It was found that when £ was set to the total 
thickness of the carpet, 6mm, that the signal attenuation tended to decrease with 
increased frequency when the frequency was greater than approximately 15 kHz. 
This decrease in attenuation was not found in the experimental data shown. Similar 
behavior was found for £ = 3 mm. When £ was set to the thickness of the hemp-glue 
thickness, 1 mm, a negative excess attenuation was predicted at frequencies below 
5 kHz; that is, the predicted pressure was greater than the 6dB pressure doubling 
predicted for a rigid surface. This result was reasonable since at low frequencies, 
ki was quite small and the surface was effectively very hard. Thus, no apparent 
improvement in the fit with experiment was found from the use of the thin layer 
correction, and since an excellent fit was found without the correction, Z 7 , not Z 2 , 
was taken as the surface impedance for the carpet-on-plywood surface. 

In the case of the bare plywood surface, initial predictions, using the single 
parameter model with a = 80 x 10 6 kg/(m 3 -s) to estimate the impedance, did not 
fit the experimental data particularly well, although it was the best fit obtainable 
(Fig. 6.4- 6.6). Although the prediction fit the experimental data at certain points 
in the frequency range, the predicted shape of the curve certainly did not follow 
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Figure 6.4 Result obtained for the determination of an effective flow re- 
sistance, <t, for the plywood surface. The curves and data represent 
SPL quantities similar to those shown in Fig. (5. 2— 5.3). The solid curve 
represents the predicted SPL curve computed using the two parame- 
ter im pedance model. In that model, an effective flow resistance of 
= 120 x 10 6 kg/(m 3 *s) was used along with an effective surface layer 
thickness of d t = 7.5 x 10" 6 m. The dashed curve represents the pre- 
dicted SPL computed using the single parameter model with an effective 
flow resistance of o = 80 x 10® kg/(m 3 -s). The listener height was 0.3 cm, 
the source height was 2.0 cm, and the horizontal distance between the 
listener and source was 100 cm. 
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FREQUENCY (kHz) 


Figure 6.5 Result obtained for the determination of an effective flow re- 
sistance, o - , for the plywood surface. The curves and data represent 
SPL quantities similar to those shown in Fig. (5.2-5 .3). The solid curve 
represents the predicted SPL curve computed using the two parame- 
ter impedance model. In that model, an effective flow resistance of 
— j20 x 10 6 kg/(m 3, s) was used along with an effective surface layer 
thickness of d e = 7.5 x 10" 5 m. The dashed curve represents the pre- 
dicted SPL computed using the single parameter model with an effective 
flow resistance of a = 80 x 10® kg/ (m 3 -s) . The listener height was 0.3 cm, 
the source height was 2.0 cm, and the horizontal distance between the 
listener and source was 80 cm. 
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FREQUENCY (kHz) 


Figure 6.6 Result obtained for the determination of an effective flow re- 
sistance, <r, for the plywood surface. The curves and data represent 
gPL quantities similar to those shown in Fig. (5.2— 5.3). The solid curve 
represents the predicted SPL curve computed using the two parame- 
ter impedance model. In that model, an effective flow resistance of 
a e = 120 x 10 6 kg/(m 3 -s) was used along with an effective surface layer 
thickness of d e = 7.5 x 10“ 5 m. The dashed curve represents the pre- 
dicted SPL computed using the single parameter model with an effective 
flow resistance of o = 80 x 10“ kg/ (m 3 *s) . The listener height was 0.3 cm, 
the source height was 9.0 cm, and the horizontal distance between the 
listener and source was 100 cm. 
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that of the data. In two of the figures, a very shallow dip in the data below 7 kHz 
was followed by a small rise until approximately 16 kHz. Past that frequency, the 
data tended to fall off more rapidly until roughly 35 kHz at which point another 
rise and fall began. (Due to the limited spectrum of the free field incident wave, 
the data past 35 kHz was not to be trusted.) A major difficulty with this data was 
in how to interpret the small rise in SPL between 7 and 16 kHz. In this range, 
the SPL of the free field incident wave was large enough so that numerical error 
associated with division by small numbers was negligible. Also, it did not seem to 
be an interference effect because the pathlength difference between the direct and 
reflected ray paths was on the order of 0.5% of a wavelength at these frequencies. 
An interesting fact about the small rise is that it began at the 6dB line. Therefore, 
the rise appeared to be a negative excess attenuation. As mentioned earlier, this 
effect can occur in special cases when the surface wave component of the solution 
dominates at a point in the field. A surface wave was shown to exist whenever (6.16) 
was satisfied. It can be shown for frequencies around 10-15 kHz, that the impedance 
estimate obtained from the single parameter model, when a = 80 x 10® kg/(m 3 s), 
satisfied the requirement for a surface wave . However, for such an estimate, the 
magnitude of the surface wave can be shown to have been small. This result was 
not surprising because the impedance model of Delany and Bazley fixed the ratio 
of the real and imaginary parts of Z 2 such that they were of comparable value even 
though X 7 was usually slightly greater than as required by the surface wave 
existence condition (6.16). This fact was important because the magnitude of the 
surface wave was shown to fall off exponentially with R 2 , the resistive component 
of impedance. Thus, for a surface of large Z 2 , the single parameter model automat- 
ically predicted comparably large resistive and reactive components. Therefore, the 
model could not have simultaneously predicted a large impedance and still allowed 
for a significant surface wave contribution. In the case of a plywood table top, a 
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fairly large value of Z 2 was expected, and, as shown in the figure, the results based 
upon the single parameter model did not predict the negative excess attenuation 
found in the measured data. 

For hard surfaces, such as the plywood surface, a suggested alternative to 
the single parameter model was found from Attenborough’s five parameter model, 
mentioned briefly above. His idea was to model a hard surface as a thin layer of 
thickness, <f, atop a rigid backing. Then for sufficiently small kd, the impedance Z 2 
was shown to modelled by 

Zi , . .38.99 

— — - 0.82 x 10" 5 a, d e + t-—- (6-23) 

Pi Ci /“« 

where d e is an effective thickness, in meters, and o e is an effective flow resistance, 
in kg/(m 3 -s). 34 In this model, the real part can be small even when the reactive 
part is large. Such a formulation was suggested in the work of Donato, who showed 
for a rigidly backed thin surface of constant porosity that at low frequencies the 
impedance was nearly purely imaginary and scaled as l/f. Donato also showed 
that this result was identical to that found for a surface with a porosity which 
decreased rapidly with depth into the surface (although Attenborough disputed 
Donato’s conclusion). 35 This two parameter model was in accord with physical 
intuition in a number of regards. First, for large d t , the model gives a real part of 
impedance which is much larger than the imaginary part. This result agrees with 
the notion that deep pores admit more flow and, thus, higher losses due to flow 
resistance, than do shallow pores which tend to act more like capacitors which do 
not attenuate the field. Second, at higher frequencies, the model gives a real part 
which is again much larger than the imaginary part. This result reflected the notion 
that, at higher frequencies, the width and depth of the pores was effectively larger, 
and the pores again admit more flow. Thus, the effect at higher frequencies was 
to make the surface seem softer and more like an absorber. An examination of the 
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plywood table top indicated that the two parameter model might give a reasonable 
estimation of the plywood impedance. 

The method for deducing the parameters was similar to that used before except 
that two parameters, a, and d c , had to be recursively estimated and plugged into 
the formula. Since the single parameter model yielded a « 80 x 10 6 (MKS), a e was 
initially held to be on the order of o. As before, the two parameters were varied until 
the predicted SPL approximated that of the experimental SPL. Figs.(5.4-5.6) show 
the predicted SPL based upon the single parameter and the two parameter models, 
respectively. Best fit values of o e = 120 x 10 6 kg/(m 3 *s), and d e = 7.5 x 10" 6 m 
were found for the three sets of experimental data. Note that the measured negative 
excess attenuation was better approximated by the theory when the two parameter 
model was used. The best fit values were chosen based upon the fits with all three 
data sets. Even though the results do not support the local reaction assumption, 
they do indicate that the two parameter model was a much better approximation 
for the plywood surface because it predicted the qualitative behavior of the data 
even though it had difficulty matching the data itself over a range of incident angles. 
In sum, the impedance, Z 2 , as estimated by the two parameter model, was accepted 
as the impedance of the plywood surface. 


6.4 Summary 

Apparently good estimates of the surface impedance at grazing angles of incidence 
were found by applying the theory for a point source radiating near an impedance 
plane. The single parameter model of Delany and Bazley provided a working model 
for the impedance of the carpet-on-plywood surface. The single parameter known as 
the effective flow resistance was estimated to be 800 x 10 s and 1600 x 10 s kg/ (m -s) at 
high and low frequencies, repectively, for the carpet-on-plywood surface. However, 
the theory encountered difficulty in predicting the field near the plywood surface 
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when the single parameter model was used. Instead, a two parameter model sug- 
gested by Attenborough and Donato for thin layers atop a rigid back was used to 
estimate the impedance. The two parameters were an effective layer thickness, d e , 
and, again, an effective flow resistance, a e . Predictions of the field, obtained using 
the two parameter model with o e = 120 x 10 6 kg/(m 3 -s) and d e = 7.5 x 10" 6 m, 
improved upon the predictions found using the single parameter model. Specifically, 
the negative excess attenuation measured near the plywood surface was predicted 
when the two parameter model was used. Such a prediction was practically prohib- 
ited by the single parameter model. 
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CHAPTER VII 

MEASUREMENT 
OF THE 

ACOUSTIC FIELD ON THE CURVED SURFACE 


This chapter describes the results of experiments conducted to measure the inser- 
tion loss on the curved surface. Data was obtained at a series of points spanning 
the curved surface from front to back along the ridge centerline. From this data, 
experimental values of insertion loss were deduced at each point. The experimental 
procedure was that as described in chapter 4. Predictions of the insertion loss at 
these same points were obtained by computing (3.24). The frequency content of 
the incident waveforms was such that 200 < kR < 1800. Discrepancies between the 
experimental data and the predicted data were often less than 2 dB in those regions 
where the signal-to-noise ratio was sufficiently large. The following is devoted to 
detailing important aspects of the predicted and measured fields on the surface, and 
the differences between them. 
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7.1 Equations for the Insertion Loss on the Surface 

An implication of the Kirchhoff-Helmholtz Integral theorem and of the Geometrical 
Theory of Diffraction is that if the incident field and the field on the surface of an 
isolated scatterer are known, then the solution for the field surrounding the scatterer 
is uniquely determined. Therefore, the determination of the field (or its associated 
insertion loss) at the surface was of fundamental interest and was a natural starting 
point for an investigation of the surrounding field. 

As presented in chapter 3, the magnitude of the pressure, |p|, at a point in the 
field is proportional to the magnitude of the Fock-vaji-der-Pol-Bremmer function 


= A= [ [v(a — 17 ) Tf~) — ^^Kw l (a-v)]e iat da (3.21) 

y/KJ-oo 1 * K («) “ 9 U>! (a) -I 

which is reexpressed here for convenience. A form of G valid on the surface 
and near the apex was shown to be 
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Further examination of (3.21) revealed that another form of G valid on the 
surface is 


1 f°° c"* 


is 


(7.1) 


Thus, for the case of a completely rigid surface, q = 0, G can be expressed as 


1 f°° e Ui 

G(£,0,0) = -p/ -rrrds 
V*J- CO KW 


(7.2) 


In this form, |G(£,0,0)| can be shown to approach 2 (or 6dB) for £ large and 
negative. Of particular interest was |G(0,0,0)| = 1.399 (or 2.92 dB) which gives 
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the value of the field at the apex of a rigid curved surface. As seen from (7.2), 
this value is independent of frequency. (In general, |G| is a function of frequency 
through the parameters q , £, smd rj . However, at the apex of a rigid surface all 
three parameters are identically zero.) Interestingly, the value at the apex is not 
the —6 dB predicted for a plane wave incident on a flat rigid boundary, nor is it 0 dB 
(i.e. equal to the magnitude of the incident field). In Electromagnetic Diffraction 
and Propagation Problems , Fock lists the values of |£r(£,0,0)| for —4.5 < £ < 4.5. 
Values for G(£,0,g) were computed using Simpson’s integration rule. The values 
of q were chosen to be those deduced in chapter 6 for the plywood and carpet- 
on-plywood surfaces. The computations were limited to dimensionless arc lengths 
corresponding to |£| < 2. 

7.2 Particulars of the Measurements 

Experiments were conducted to measure the insertion loss along the surface of the 
plywood ridge with and without the thin carpet overlay. The experimental method 
and data reduction was the same as that described in chapter 4. In short, the field 
microphone was laid upon the ridge surf ace with the microphone cartridge housing 
touching the surface. In this way, the center of the microphone diaphragm was 
approximately 1/8 in (or « 3 mm) above the surface. In the case of the carpet-on- 
plywood surface, the field microphone was not buried in the carpet but was able to 
rest on top of the carpet fibers. 

Originally, all data was captured with the normal to the microphone diaphragm 
tangent and parallel to the ridge centerline. However, it was found that such an 
orientation introduced significant microphone sensitivity effects at points of large 
and negative £. Therefore, at points of £ < 0, the microphones were oriented such 
that the diaphragm normal was tangential but transverse to the ridge centerline. 
Reference waveforms were captured as described in chapter 4 but always such that 
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the reference microphone shared the same orientation as the field microphone. In 
this way, a reference waveform and field waveform pair was. captured over much 
of the span of the ridge. In fact, for the plywood surface data was captured every 
10 cm from s = -80 cm to s = 100 cm, where s is the arclength from the apex and is 
positive in the shadow zone. At s = 100 cm, the field waveform was not significantly 
larger than some of the background noise. Contamination by noise was an even 
bigger problem for the data captured on the carpet-on-plywood surface since the 
softer surface automatically made for small signals. Thus, for this surface, data was 
captured every 5 cm from s — — 60 cm to s = 60 cm. Past s — 60 cm, the field signal 
was practically indistinguishable from the background noise for this surface. 

The measurement error associated with the arc length was estimated to be less 
than 5 mm including the influence of deviations in the curved from that of a perfect 
circle or parabola. The dimensionless arc length was defined as 



and so d£ a ds. In this way, the error in £ was estimated to be approximately 0.01 
at 10 kHz. 


7.3 Comparison of the Predicted and Measured Insertion Losses 

The computational and experimental results for both types of surface and for fre- 
quencies of / = 10, 20, and 30 kHz are shown in Figs.(6.1-6.3). In regard of the 
plywood surface, the insertion loss curves are approximately the same as those for a 
rigid surface. The curves for the plywood surface also have little dependency upon 
frequency. Mathematically, this is because |^| ^ 1 in these cases, and so the value 
of G is largely determined by w[ (s), which, of course, is independent of frequency. 
Physically, a simple inspection of the plywood surface revealed that the dimensions 
of the surface pores of the plywood were much smaller than any wavelength within 
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Figure 7.1 Insertion loss on the curved plywood and carpet-on-plywood 
surfaces at / = 10 kHz. The discrete symbols (x) and (0) represent 
the measured insertion loss on the two surfaces, respectively. The solid 
lines represent the predicted insertion losses when the effective flow re- 
sistances are o = 1.6 x 10® and a = 120 x 10* MKS. The dashed curve, 

( represents the predicted insertion loss for a rigid surface. £ is a 

dimensionless expression for the arc length, s. 
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Figure 7.2 Insertion loss on the curved plywood and carpet-on-plywood 
surfaces at / = 20kHz. The discrete symbols (x) and (<>) represent 
the measured insertion loss on the two surfaces, respectively. The solid 
lines represent the predicted insertion losses when the effective flow re- 
sistances are o = 0.8 x 10 6 and a = 120 x 10* MRS. The two other 
dashed curves represent the predicted insertion loss for a rigid surface, 

( and for a surface with effective flow resistance of 1.6 x 10® MRS, 

( ). £ is a dimensionless expression for the arc length, s. 
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o Carpet on Plywood , 800 cgs 



Figure 7.3 Insertion loss on the curved plywood and c&rpet-on-plywood 
surfaces at / = 30kHz. The discrete symbols (x) and (0) represent 
the measured insertion loss on the two surfaces, respectively. The solid 
lines represent the predicted insertion losses when the effective flow re- 
sistances are o = 0.8 x 10® and o = 120 x 10* MRS. The two other 
dashed curves represent the predicted insertion loss for a rigid surface, 

( ), and for a surface with effective flow resistance of 1.6 x 10® MRS, 

( ). £ is a dimensionless expression for the arc length, s. 
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the frequency range. In fact, by the methods of chapter 6, the plywood surface was 
found to behave as a thin porous layer backed by a rigid surface, or equivalently, as 
a porous medium with a rapid variation in porosity with depth into the medium. 
The effective flow resistance was found to be very large. This result implied either 
a very small porosity or very small pore dimensions relative to a wavelength. In 
this case, the real part of the impedance, or resistance, was shown to be effectively 
independent of frequency while the imaginary part of the impedance, or reactance, 
was inversely proportional to frequency, although much larger than the resistance 
except at very high frequencies. It is well known that the resistance is related to 
the am ount of acoustic energy absorbed by the surface. For instance, in the case of 
plane wave reflection, the time-averaged acoustic power flowing into the surface per 
unit area is proportional to the resistance. In the case of spherical wave reflection 
one expects a similar relationship. Therefore, the absorption of the plywood surface 
is expected to be relatively independent of frequency for the range of frequencies 
considered. From this perspective, it follows that the insertion loss on the plywood 
surface should be relatively independent of frequency as well. 

In the case of the carpet-on-plywood surface, the insertion loss is seen to be a 
strong function of frequency. This result is in keeping with the fact that \q\ is on 
the order of unity or larger. Also, by inspection, it was clear that the porosity and 
pore dimensions of the carpet were much larger than those of the plywood although 
still smaller than a wavelength. Again, in chapter 6, it was found that the semi- 
empirical impedance model of Delany and Bazley for fibrous absorbent materials 
provided a good match with the experimental data obtained with the carpet-on- 
plywood surface. In their model, the resistance and reactance are relatively strong 
functions of frequency. Therefore, the scaled admittance, |^|, is not small but is 
a strong function of frequency. Thus, the insertion loss on the carpet-on-plywood 
surface should be a strong function of frequency. 
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The general nature of the curves show that the insertion loss of a hard surface 
is approximately — 6dB at £ — —2 and varies relatively slowly with f when ^ < 0. 
However, when £ > 0, that is for points in the shadow zone, the insertion loss 
increases steadily with £. For softer surfaces, and for £ < -1, the value of the 
insertion loss is relatively independent of { but is a strong function of frequency 
and, in fact, increases with frequency. At points of positive £, in the shadow zone, 
the insertion loss increases quickly with £. The slope of the curve in the region £ > 0 
also appears to increase with frequency. This is in keeping with the concepts of the 
Geometrical Theory of Diffraction which predicts that the attenuation of creeping 
waves is a function of frequency. 

Figs. (6.2-6. 3) show the effect on the predicted insertion loss on the carpet- 
on-plywood surface due to a change in the effective flow resistance, o e , from Oq 
1.6 x 10 6 (MKS) to a„ = 0.8 x 10® (MKS). As discussed in chapter 6, a 0 was the 
best fit, low frequency flow resistance while was the adjusted high frequency 
flow resistance. The use of a frequency dependent flow resistance, o j , was justified 
in chapter 6. Since a = a 0 until 14.5 kHz, only one curve of the insertion loss above 
the carpet-on-plywood surface was shown in Fig. (6.1). The use of the frequency 
dependent o resulted in noticeably improved predictions of the measured insertion 
losses. 

An interesting result for the plywood surface is that the predicted insertion loss 
at the apex is 0.42 dB Jess than that for a perfectly rigid surface. This effect is no 
doubt related to the negative excess attenuation effect predicted and measured for 
the flat plywood surface when the source and receiver were close to the surface, that 
is, when the angle of incidence is near grazing. The apex, of course, is defined as 
the point of grazing incidence so the effect appears to be in accord with the theory 
discussed in chapter 6. 

At 10 kHz, there is excellent agreement between the predicted and measured 
insertion losses in the region £ < 0. In fact, the agreement would be good for a 
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wide range of large effective flow resistances, o. However, in the region £ > 0, the 
sensitivity of the predicted insertion loss to the magnitude of o is not negligible. 
Not surprisingly, when £ > 0, the agreement is not quite as good as when £ < 0, 
particularly at the larger values of £. This result was partially explained by the 
decreased signal-to-noise ratio at large positive £. Still, the figures show that, when 
£ > 0, the predicted increases in insertion loss, as £ or / increases, do not keep 
pace with the corresponding measured increases in insertion loss. Lower values of 
a would be required to match the measured data. 

Remember that the results of chapter 6 indicated that the plywood surface is 
not locally reacting and that its best fit a was a factor of ten lower for slightly larger 
angles of incidence (i.e. slightly larger than grazing incidence). Further, the field 
at a point on the surface was certainly influenced by the impedance of all points 
of the surface which lay forward (toward the source) of that point, although the 
impedance of points near the apex probably had more influence on the insertion 
loss than did the impedance at points farther forward. From this perspective, it 
is not surprising that the measured surface data for £ > 0 is more like what one 
would predict if a e was some average of o e ' s found at various angles near grazing. 
One other consideration is that the actual plywood surface of the ridge was three 
times as thin as the plywood surface of the flat table top, which, of course, was the 
surface used to deduce the plywood impedance. This difference also may account 
for part of the indication that the curved surface is softer than the flat table top. 

The data for the carpet-on-plywood surface was perhaps more interesting be- 
cause of its strong dependence on frequency. The general shape of each curve was 
the same and similar to that found for the plywood surface. However, the inser- 
tion loss increased with frequency for all points of £. This result is in accord with 
the predictions of the theory and with what might have been expected by physical 
intuition. The agreement between the measured and predicted insertion losses was 
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again excellent at 10 kHz, but the predictions based upon the low frequency effective 
flow resistance, <j 0 i underestimate the measured insertion loss at the higher frequen- 
cies. However, when the predictions are calculated for , then the agreement with 
measurement was improved. 

In conclusion, when one considers that all measurements were likely to contain 
an inherent error of at most ±2dB, the experimental data is in excellent agreement 
with that predicted by the theory even when both surfaces were modelled as fibrous 
absorbent materials with a = 1.6 x 10* and a = 80 x 10* kg/(m 3 -s) for the carpet- 
on-plywood and plywood surfaces, respectively. More sensible considerations of the 
surface structure and its associated impedance were not detrimental to the match 
between experiment and prediction. In the case of the carpet-on-plywood surface, 
the inclusion of a frequency dependent flow resistance improved the match at higher 
frequencies. In the case of the plywood surface, the inclusion of the thin, porous 
layer model did not improve the match, but did not significantly degrade the match 


either. 
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CHAPTER VIII 

MEASUREMENT OF THE ACOUSTIC FIELD 
ABOVE AND BEHIND THE CURVED SURFACE 

Just as chapter 7 was devoted to an analysis of the theoretical and experimental 
results found for the case of a listener on the curved surface, this chapter is devoted 
to an analysis of the results found when the listener is located above or behind the 
curved surface. However, this chapter is more extensive than chapter 7 because of 
the need to analyze the results at various distances from the apex in the bright 
zone, pen umb ra, and shadow zone. These regions were defined in chapter 3 as 
those regions located well above the curved surface, behind the surface yet near the 
line of sight, and behind the surface but well below the line of sight, respectively. 
The solution valid in all three regions was expressed in terms of the Fock-vaji- 
der-Pol-Bremmer function, G. An alternate form of this function was used to 
calculate the theoretical results at the surface, as presented in chapter 7. In this 
chapter, theoretical results which were obtained from an approximate form of G are 
presented. In addition, theoretical results which were obtained directly from the 
geometrical acoustics solution, the creeping wave series solution, and the knife-edge 
diffraction solution axe presented. 

The analysis and comparisons of the theoretical and experimental data are 
similar to those made in chapter 7. Experimental data was captured along vertical 
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lines located at x = 0, 63, 120, 200, and 400 cm downstream from the apex of 
the curved surface. These distances correspond roughly to 0, Rj 4, R/ 2, 4Rj5, 
and 8/2/5. As before, both origins of the (x,y) and (£,»;) coordinate systems were 
located at the apex. Enough data was captured at each x location such that, at a 
given frequency, a fairly well resolved curve of insertion loss versus distance above 
and below the line of sight was obtained. Typically, an incident acoustic pulse was 
such that reasonable levels of its Fourier transform ranged from 5—40 kHz. The 
data presented here was limited to frequencies of 10 and 20 kHz. At each of these 
frequencies, the magnitude of the Fourier transform of a typical incident pulse was 
near that of the peak magnitude of the entire transform. The intent was to minimize 
errors associated with low signal-to-noise ratios while also providing a fair range or 
ratio of frequencies. For the laboratory value of R — 2.5 m and c = 345.6 m/s, these 
values of frequency corresponded to kR = 454 and 909, respectively. 

In general, the results showed a good comparison between theory and exper- 
iment. This was particularly true of the data captured at x = 0 and x — 63 cm. 
At these locations the discrepancy was less than 2dB in most cases. However, at 
points further away, the discrepancy was found to be as large as 5dB in some cases. 
A more detailed assessment of all the results follows. 

Measurement errors were of the same order as those discussed in chapter 7, 
and the derivation of error ranges was done in the same manner as well. 

8.1 Insertion Loss in the Bright Zone 

The ability of the plane wave geometric acoustics solution to predict the insertion 
loss at points in the bright zone is important because this solution is conceptually 
and computationally simple. In addition, an asymptotic version of it served as the 
radiation boundary condition in the derivation of G. In the bright zone, then, the 
accuracy of G should be comparable to that of the this geometric acoustics solution. 
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Thus, it seems appropriate to compare the plane wave geometric acoustics solution 
to the experimental data before presenting and discussing the results involving 
G. Further, since the spark source simulated a point source in the experiments, 
the geometric acoustic solution for a point source is also presented. This solution 
should yield a better comparison to experiment than does the plane wave version, 
and as such can be used as a test of the reasonableness of the experimental data. 

The insertion loss above a curved surface should be well approximated by both 
geometrical acoustics solution so long as the receiver is in the bright zone; that 
is, not too close to the curved surface or the line of sight, which are caustics. In 
chapter 3, the geometric acoustics solution for an incident plane wave was shown 

to be 

p = Pie ikx + P,[A(0)/A(£)] 1/3 Ke <fc *‘e <k ‘ (3.1) 


When the listener is in the vicinity of the apex, a simpler approximate form of the 


solution was written as 
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where 


Q = [(4/9)z 3 + (2/3 )Ry} 1 ' 7 (3-3) 

= (2 k/R 2 )[Q 3 - (8/27)z 3 - [2jZ)Rxy\ (3.4) 

An asymptotic form of (3.2) was the radiation boundary condition, referred to 
above, used to solve the boundary value problem of chapter 3. In all of the figures, 
(3.2) was referred to as the “plane wave geometric solution” . 

As mentioned above, the experimental data was generated by a point source, 
namely the spark source. The point source solution is developed in appendix C and 
is similar to that for a plane wave. This similarity is expected since the plane wave 
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solution is equivalent to the point source solution when the point source is moved 
to infinity. In short, the point source solution is written as 



where 5 is the source strength. 

One aspect of these solutions which should appear in the calculated data has 
to do with the relative phase of the point source solution compared to the plane 
wave solution. In the case of a plane incident wave, the phase of the direct (non- 
reflected) ray is independent of the listener height, y, while, in the case of a point 
source, the phase of the direct ray increases with an increase in y. In both solutions, 
the phase of the reflected ray increases with an increase in y; but since the point 
source is located relatively far from the apex, the phases of the respective reflected 
rays are nearly equal. Further, in the case of a point source, it is easy to see from 
the ray geometry that an incremental change in y produces a greater net change in 
the reflected ray phase than in the direct ray phase. This fact is obviously true in 
the case of a plane wave. Thus, one can conclude that an incremental change in y 
produces less of a net change in phase in the case of a point source solution than 
in the case of a plane wave. Consequently, the oscillation rate with respect to y of 
the plane wave solution should be higher than that of the point source solution. 

Both geometrical acoustics solutions were derived under the assumption that 
the listener is in the vicinity of the apex but atleast several wavelengths away 
from the surface. Therefore, it is reasonable to expect the best agreement with 
experiment to occur at small rather than large values of x. Consider the data at 
i = Ocm and / = 10 and 20 kHz as shown in Figs.(8.1-8.2). 

As discussed above, the plane wave solution did indeed oscillate more rapidly 
than the point source solution. Another interesting aspect of both geometric solu- 
tions was that the magnitudes of adjacent peaks and troughs are similar but not 
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Figure 8.1 Insertion loss above the curved plywood surface at x = Ocm 
and / = 10 and 20kHz, respectively. The discrete symbols (x) represent 
the insertion loss measured above the laboratory scale ridge. The solid 
and dashed lines represent the point source and plane wave geometric 
solutions, respectively. Y is the dimensionless height (3.29). 
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Figure 8.2 Insertion loss above the curved carpet-on-plywood surface at 
x = Ocm and / = 10 and 20 kHz, respectively. The discrete symbols (0) 
represent the insertion loss measured above the laboratory scale ridge. 
The solid and dashed lines represent the point source and plane wave 
geometric solutions, respectively. Y is the dimensionless height (3.29). 
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identical. In fact, in the case of the carpet-on-plywood surface, the peak and trough 
magnitudes decrease slightly with height, while in the case of the plywood surface, 
the peak and trough magnitudes increase slightly with height. An explanation of 
this result was found in the functional dependence of the reflection coefficient on 
the angle at which the reflected ray struck the surface. As y increased, the angle of 
incidence decreased, and in the case of the plywood surface, the magnitude of the 
reflection coefficient increased, thus yielding larger peak to trough variations as y 
increased. 

In the figures, the abscissa is Y which is related to but distinct from y. As 
was shown in chapter 3, Y is the dimensionless height found in the knife-edge 
diffraction formula (3.28) which was derived from the Fock-van-der-Pol-Bremmer 
integral through the use of the method of steepest descents. That is, Y is meaningful 
for the general solution, G, and was, as will be soon shown, used as the abscissa in 
plots of those results. Y is not necessarily more meaningful than y in the case of 
the geometric acoustics solutions. However, for consistency, Y was also used as the 
abscissa in the plots of the geometric acoustics results. 

Consider the case of the carpet-on-plywood surface. The agreement between 
the experimental data and the point source solution is excellent when Y > 1 . In 
fact, in this region, the experimental data for the carpet-on-plywood surface fit 
the point source curve to within 1/2 dB at each of / = 10 and 20 kHz. The plane 
wave curve also agrees rather well with the experimental data when 1 < Y < 2. 
However, when Y > 2 the agreement becomes progressively worse as Y increases. 
This is as expected because, as discussed above, the plane wave solution does not 
accurately represent the net change of phase associated with the rays emanating 
from the actual source used in the experiments. In the region Y < 1, both the 
point source and plane wave curves fail to maintain the same accurate fit with the 
experimental data. This is not surprising since neither solution is expected to be 
valid for small Y. 
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In the case of the bare plywood surface, the geometric solutions compare well 
with experimental data although not as well as in the case of the carpet-on-plywood 
surface. In fact, an eyeball best fit of the experimental data would appear to 
fall in between the plane wave and point source curves, although for Y > 3, the 
experimental data is clearly better fit by the point source curve. Part of the problem 
with the plywood surface data is that there is not as much of it as there is of 
the carpet-on-plywood surface data. This is particularly evident in the results at 
/ = 20 kHz. However, it also appears that the experimental data is simply not 
as well fit by the theory as in the carpet-on-plywood surface case. An interesting 
result is that both solutions fit the experimental data rather well at small values of 
y. This result is due to the fact that the plywood is hard enough that, at listener 
locations near to the surface, the insertion loss is expected to be around -3dB. In 
chapter 7, the insertion loss measured at the apex of the plywood surface was also 
approximately -3 dB. Thus, the good comparison between the geometric solutions 
and the experimental data for x = 0 and small Y is not surprising for the plywood 
surface. However, such a comparison can only be anticipated when the surface is 
hard. 

At x = 63cm, Fig. (8. 3-8. 4), the results are very similar to those at x = 
0 cm. Again, for the carpet-on-plywood surface, the point source results fit the 
experimental data extremely well when Y > 1 but the fit becomes progressively 
poorer as Y gets small. For the bare plywood surface, the results mirror those at 

x = 0 cm. 

At x = 120 cm, the results are somewhat different from those at x = 0 and 
63 cm. In the case of the plywood surface, Fig.(8.5) shows that the point source so- 
lution follows more closely the experimental data than does the plane wave solution 
although both solutions appear to oscillate faster than what one might envision as 
being described by the data. In the case of the carpet-on-plywood surface, Fig. (8.6), 
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Figure 8.3 Insertion loss above the curved plywood surface at x = 63 cm 
and / = 10 and 20 kHz, respectively. The discrete symbols (x) represent 
the insertion loss measured above the laboratory scale ridge. The solid 
and dashed lines represent the point source and plane wave geometric 
solutions, respectively. Y is the dimensionless height (3.29). 
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Figure 8.4 Insertion loss above the curved carpet-on-plywood surface at 
x = 63 cm and / = 10 and 20 kHz, respectively. The discrete symbols 
(<>) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the point source and plane 
wave geometric solutions, respectively. Y is the dimensionless height 
(3.29). 
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Figure 8.5 Insertion loss above the curved plywood surface at x = 120 cm 
and / = 10 and 20kHz, respectively. The discrete symbols (x) represent 
the insertion loss measured above the laboratory scale ridge. The solid 
and dashed lines represent the point source and plane wave geometric 
solutions, respectively. Y is the dimensionless height (3.29). 
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Figure 8.6 Insertion loss above the curved ca rpet-on-plywood surface at 
x = 120 cm and / = 10 and 20 kHz, respectively. The discrete symbols 
(<>) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the point source and plane 
wave geometric solutions, respectively. Y is the dimensionless height 
( 3 . 29 ). 
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the plane wave solution appears to overestimate the oscillation rate indicated by the 
data. However, the point source solution appears to underestimate the oscillation 
rate of the data. That is, the point source solution reaches its first trough at a Y 
which is larger than the Y of the data’s first apparent trough. 

It is clear that both geometric acoustic solutions compare worse with the x = 
120 cm data than they do with the x = 0 or 63 cm data. This result is not surprising 
since both geometric acoustics solutions assume that the radius of curvature, R — 
2.5 m, is much larger than any of the other relevant dimensions. Obviously, 250 cm 
is not that much larger than 120 m. In addition, intuition leads one to expect 
greater diffraction effects at larger x since for a given y, points of larger x are in 
closer angular proximity to the line of sight than axe points of smaller x. This idea 
is depicted by the well-known penumbra cone emanating from the apex of a curved 
surface. Perhaps, then, the experimental data at x = 120 contains more significant 
diffracted contributions than does the data at x = 0 and 63 cm. (This possibility is 
examined in the next section.) As with the x = 0 and 63 cm results, the results at 
x = 120 cm indicate that, when Y < 1 , both geometric acoustic solutions compare 
poorly with the measurements made over the carpet-on-plywood surface. 

The results at x = 200cm, Fig.(8.7-8.8) are quite similar to those at x = 
120 cm. Again, the point source solution compared better with the experimental 
data than did the plane wave solution. Also, as at x = 120 cm, the apparent oscilla- 
tion rate of the experimental data was between that of the point source and plane 
wave solutions in the case of the carpet-on-plywood surface, and was greater than 
that of both solutions in the case of the plywood surface. Some of the discrepancy 
may be due to diffraction effects which are discussed in the next section. 
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Figure 8.7 Insertion loss above the curved plywood surface at z = 200 cm 
and / = 10 and 20kHz, respectively. The discrete symbols (x) represent 
the ins ertion loss measured above the laboratory scale ridge. The solid 
and dashed lines represent the point source and plane wave geometric 
solutions, respectively. Y is the dimensionless height (3.29). 
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Figure 8.8 Insertion loss above the curved carpet-on-plywood surface at 
x = 200 cm and / = 10 and 20 kHz, respectively. The discrete symbols 
(0) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the point source and plane 
wave geometric solutions, respectively. Y is the dimensionless height 
(3.29). 
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8.2 Insertion Loss in the Penumbra and the function G 

Near the apex and along the line of sight, the basic concepts underlying the geo- 
metric acoustic solutions fail to be valid. So, it was not surprising to find that these 
solutions also failed to consistently yield an accurate prediction of the field for small 
Y. In contrast, the Fock-van-der-PoI-Bremmer function, G, of chapter 3, which is 
restated here as 


i /" [•(— ) - *■ 


was derived for small Y ; that is, at listener positions near the apex or in the not too 
distant penumbra. G was also shown to asymptotically approach the plane wave 
geometric acoustic solution for small Y. Thus, G should yield a good comparison 
with the experimental data when Y < 1 and should approach the plane wave 
geometric acoustic solution as Y increases. The form of G used in the calculations 
was the saddle point approximation 


G = $ + 


(3.26) 


where 


* = _ J- {" ~ 
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and 
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1 


(3.28) 


where 


y = _0) ,/ %*/*(f--? ,/J ) ( 3 - 29 > 

This form was chosen because it offered a ready physical interpretation as well as a 
meaningful dimensionless parameter, Y . As discussed in chapter 3, V is viewed as a 
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dimensionless height and is proportional to where (£,*?) are parabolic 

coordinates. Since ij > 0 everywhere off the surface, Y = 0 represents any and all 
points on the line = rj, which, for increasing z is an upward sloping parabola 
in the (x, y) coordinate space. Although Y is dimensionless, one could not a priori 
expect Y to normalize the results at different x positions. This is true because G 
is an explicit function of the coordinates (£,rj) as well as of Y . Therefore, since 
an infinite number of (£,r)) coordinate pairs exist which correspond to any given 
value of Y, many different magnitudes of the background terms can exist for any 
given value of Y . However, for a given x and /, a one-to-one mapping does exist 
between Y and y. Thus, in those figures where x and / are specified, the results 
would appear the same if y was used in place of Y . 

Figs.( 8.9- 8.10) show the results at x = Ocm when / = 10 and 20kHz, re- 
spectively. The figures show clearly that, for both surfaces, G does, in fact, match 
the plane wave geometric acoustic solution when Y > 1. However, as seen by the 
progressive separation of the curves, the geometric solution has a slightly slower os- 
cillation rate with Y than does G. The discrepancy is not significant for 1 < Y <3, 
and is probably due to the fact that G matches the geometric solution when Y is 
small, but not necessarily when Y is not small, such as when Y > 3. The discrep- 
ancy appears to be slightly dependent on the value of q. For instance, in the case 
of the carpet-on-plywood surface, the match is nearly exact for Y < 3, while, in the 
case of the plywood surface, the match is on the order of 1 dB or less. Of course, 
then, the match with experiment, when Y > 1, is roughly the same as was discussed 
for the plane wave geometric solution. At the larger Y , the curve of G does match 
the experimental data better than does the plane wave geometric acoustic curve. 
However, this result is probably coincidental since the geometric acoustic solution 
is expected to be more valid than G when Y is large. 

The most significant aspect of the figure is that, when Y < 1, G manages to 
match the experimental data to within 2—4 dB, and in most locations to within 2 dB. 
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Figure 8.9 Insertion loss above the curved plywood surface at x = 0 cm and 
f — io and 20kHz, respectively. The discrete symbols (x) represent the 
ins ertion loss measured above the laboratory scale ridge. The solid and 
dashed lines represent the knife-edge diffraction plus background for- 
mula, (3.26-3.28), and the plane wave geometric solution, respectively. 
Y is the dimensionless height (3.29). 
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Figure 8.10 Insertion loss above the curved carpet-on-plywood surface at 
x = Ocm and / = 10 and 20kHz, respectively. The discrete symbols 
(O) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the knife-edge diffraction 
plus background formula, (3.26-3.28), and the plane wave geometric 
solution, respectively. Y is the dimensionless height (3.29). 
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This result holds true for both types of surface and at both frequencies, / = 10 and 
20 kHz. As was noted previously, the geometric acoustics solutions were expected 
to fail in this region although the results for the bare plywood surface were still 
very good until Y < 1. The results for the carpet-on-plywood surface show the 

geometric solution’s failure more clearly. 

In chapter 7, it was shown that an alternate form of G gave very good agreement 
with experiment for most listener locations on the surface. So it was not surprising 
that G also matched the experimental data when z = 0 and Y were small. In fact, 
G was derived under the assumption of z, y R and was expected to be valid near 
the surface. However, it was unclear how valid G remained at x on the order of 
or greater than R. It was expected that G would compare less and less well with 
experiment as z increased, but the extent of this departure could not be ascertained 
a priori. 

Figures( 8.11-8.12) show the results when / = 10 and 20kHz at z = 63cm, 
which was roughly a distance of RjA downstream from the apex. Again, when 
Y > 1, G matched the geometric acoustics solution to within ±2 dB for the plywood 
surface and to within il dB for the carpet-on-plywood surface. Further, for both 
surfaces, G matched the experimental data nearly point-to-point for -1 < Y < 
1. For the carpet-on-plywood surface, there was one experimental data point at 
y « —0.3 which was noticably removed from the theoretical curve and from the 
apparent trend of the other experimental data. Since it stood alone and the other 
data were rather well aligned, it was assumed that this data point was spurious. 

The significant result was that G had produced an excellent prediction of the 
experimental data near a caustic, in this case, the line of sight, as well as at points 
in the lower portion of the bright zone. It should be mentioned that the predicted 
insertion loss would not be significantly effected by assuming the plywood surface 
to be rigid. A similar result was found when the listener was on the surface. Thus, 
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Figure 8.11 Insertion loss above the curved plywood surface at x = 63 cm 
and / = 10 and 20kHz, respectively. The discrete symbols (x) repre- 
sent the insertion loss measured above the laboratory scale ridge. The 
solid and dashed lines represent the knife-edge diffraction plus back- 
ground formula, (3.26-3.28), and the plane wave geometric solution, 
respectively. Y is the dimensionless height (3.29). 
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Figure 8.12 Insertion loss above the curved carpet-on-plywood surface at 
x = 63 cm and / = 10 and 20 kHz, respectively. The discrete symbols 
(<>) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the knife-edge diffraction 
plus background formula, (3.26-3.28), and the plane wave geometric 
solution, respectively. Y is the dimensionless height (3.29). 
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there was a certain amount of liberty in choosing the effective flow resistance pa- 
rameter, a, without fear of jeopardizing the result. However, the same was not 
found to be true for the carpet-on-plywood surface. In that case, the predicted 
insertion loss was quite sensitive to the value of a\ that is, |G| was sensitive to 
a when its associated magnitude of the scaled admittance, q, was on the order of 
unity. Since the predicted insertion loss involved no free parameters, the fit between 
the predicted and measured insertion losses was indeed excellent for the carpet-on- 
plywood surface. Note that the insertion loss below the line of sight was larger 
for the carpet-on-plywood surface than for the bare plywood surface. This result 
is in accord with the co mm on sense notion of carpet as an absorber compared to 
plywood. 

At x = 120 cm, Figs.( 8.13- 8.14) show that the results here were only slightly 
less good than those at x — 63 cm. Aside from the same type of mismatch between 
theory and experiment in the bright zone, as was seen and discussed in the data 
at x = 0 and 63 cm, the predicted insertion loss at x = 120 cm was within 2 dB of 
the measured insertion loss except at points deep in the shadow zone of the carpet- 
on-plywood surface. There, the disagreement was approximately 3-5 dB. For both 
surfaces, and for listener positions in the shadow zone, the theory over predicted the 
measured insertion loss. Near the line of sight, which corresponded to Y — —0.08 
when x = 120 cm and / = 10 kHz, the agreement was within 1-2 dB for both 
surfaces. The figures also show the influence of diffracted contributions in the lower 
bright zone. This influence is seen in the separation between the geometric solution 
and G when Y > 1/2. This same effect was observed in the results at x = 63 cm 
but, there, it was not as significant. Comparisons of both G and the geometric 
solution to the measured data in the bright zone should be viewed with caution 
because the source used in the experiments was a point source, not a plane wave; 
and, as discussed above, the phase effects of these two types of sources are different 
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Figure 8.13 Insertion loss above the curved plywood surface at x = 120 cm 
and / = 10 and 20kHz, respectively. The discrete symbols (x) repre- 
sent the insertion loss measured above the laboratory scale ridge. The 
solid and dashed lines represent the knife-edge diffraction plus back- 
ground formula, (3.26-3.28), and the plane wave geometric solution, 
respectively. Y is the dimensionless height (3.29). 
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Figure 8.14 Insertion loss above the curved carpet-on-plywood surface at 
x = 120 cm and / = 10 and 20 kHz, respectively. The discrete symbols 
(<£>) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the knife-edge diffraction 
plus background formula, (3.26-3.28), and the plane wave geometric 
solution, respectively. Y is the dimensionless height (3.29). 
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and this difference can be important. Results similar to those at x = 120 cm were 
found in the lower bright zone when x = 200 cm. 

The first large disagreements between predicted and measured insertion losses 
were found in the penumbra and shadow zone at z = 200 cm. In the case of the 
plywood surface, Fig. ( 8.15) shows that G over predicted the measured insertion loss 
by 2-3 dB throughout the penumbra and shadow zone. In the case of the carpet- 
on-plywood surface, Fig. (8.16), G over predicted the measured insertion loss by 
3dB near the line of sight, and by 3-5 dB in the shadow zone. The disagreement is 
slightly more pronounced in the data at / = 20 kHz, particularly at listener points 
deep in the shadow zone behind the carpet-on-plywood surface. At the deepest of 
these points, the discrepancy approaches 10 dB. This disagreement does not appear 
to be due to a low signal-to-noise ratio because the data, as a group, does not exhibit 
an unreasonably large scatter. It is more likely that the solution, G, itself is near 
or beyond its range of applicability. An investigation of the knife-edge diffraction 
plus background expression (3.26-3.28) for G proved useful. 

Fig. (7. 19) shows the magnitudes of the knife-edge diffraction and background 
components of the total solution at different x positions. For a given Y , the magni- 
tude of the background component becomes increasingly small as x increases while 
the magnitude of the knife-edge diffraction component, which is a function of Y 
alone, is independent of x. This trend indicated that at some large x, the back- 
ground component of the total solution could be discarded leaving G « $ (3.28). 
However, $ is not a reasonable solution at large x because of the nature of the 
dimensionless height Y. As was stated earlier, Y = 0 is a parabolic curve ema- 
nating from the apex and curving upward as x increases. For example, Y = 0 at 
x = 200 cm represents a point whose height, y, is greater than that at a point of 
Y = 0 and x < 200 cm. The result of these facts is that, at large x, G predicts that 
the insertion loss at a point well above the line of sight is approximately the same 



INSERTION LOSS (dB) INSERTION LOSS (dB) 


149 




Figure 8.15 Insertion loss above the curved plywood surface at x = 200 cm 
and / = 10 and 20kHz, respectively. The discrete symbols (x) repre- 
sent the insertion loss measured above the laboratory scale ridge. The 
solid and dashed lines represent the knife-edge diffraction plus back- 
ground formula, (3.26-3.28), and the plane wave geometric solution, 
respectively. Y is the dimensionless height (3.29). 
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Figure 8.16 Insertion loss above the curved carpet-on-plywood surface at 
x = 200 cm and / = 10 and 20 kHz, respectively. The discrete symbols 
(<*►) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the knife-edge diffraction 
plus background formula, (3.26-3.28), and the plane wave geometric 
solution, respectively. Y is the dimensionless height (3.29). 
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Figure 8.17 Insertion loss above the curved plywood surface at x = 400 cm 
and / = 10 and 20kHz, respectively. The discrete symbols (x) repre- 
sent the insertion loss measured above the laboratory scale ridge. The 
solid and dashed lines represent the knife-edge diffraction plus back- 
ground formula, (3.26-3.28), and the plane wave geometric solution, 
respectively. Y is the dimensionless height (3.29). 
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Figure 8.18 Insertion loss above the curved carpet-on-plywood surface at 
x = 400 cm and / = 10 and 20 kHz, respectively. The discrete symbols 
(0) represent the insertion loss measured above the laboratory scale 
ridge. The solid and dashed lines represent the knife-edge diffraction 
plus background formula, (3.26—3.28), and the plane wave geometric 
solution, respectively. Y is the dimensionless height (3.29). 
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Plywood 



Carpet- on -Plywood 



Figure 8.19 Examples of the magnitude, in dB, of the background terms 
(3.27) relative to that of the knife-edge diffraction term (3.28) for the 
cases of the plywood and ca.rpet-on-plywood surfaces. Shown are the 
results at each of the four x positions when / = 10 and 20 kHz. Y is the 
dimensionless height (3.29). 
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as that predicted by (3.31) at a point on the line of sight. (3.31) is the solution at 
large x for the case of a thin rigid wedge. Therefore, it follows that, at similar points 
along the line of sight of a curved surface and at large x, G will predict the insertion 
loss to be approximately the same as that at a point deep in the shadow of a thin 
rigid wedge. Such a prediction runs counter to intuition in that, at large enough x, 
the field behind any curved surface, independent of the surface impedance, would 
be indentical to that behind a thin rigid wedge. That is, G does not transition 
to (3.31) as required of any solution valid at large x. The conclusion is that, at 
moderate to large x, the background contribution to G is too small to adjust for 
the discrepancy between ${Y) and the thin rigid wedge solution. The result is that, 
independent of the surface impedance, G overestimates the insertion loss behind a 
curved surface when x is large. 

The results at x = 400cm, shown in Figs.(7.17-7.18), support the above rea- 
soning. At this value of x, which corresponds to 1.6 R, the disagreement between 
G and experiment is on the order of 8-12 dB in the penumbra and shadow zone. 

It should not be inferred from the above discussion and results that the back- 
ground contribution at x = 400 cm was negligible. Nor should it be inferred that the 
background contribution always produced a better fit between G and measurement 
than that between $ (F) and measurement. Fig. (8.20) shows all of the experimental 
data at x = 63, 120, 200, and 400 cm, as well as the corresponding curves of G, for 
both surfaces when / = 10 kHz. Also included is a curve of $(30- In the case of 
the plywood surface, the background contribution is seen to reduce the predicted 
insertion loss in the penumbra and shadow zone relative to that predicted by the 
knife-edge diffraction contribution alone. However, as x increases, the measured 
insertion loss tends to shift to the left, while the predicted insertion loss shifts to 
the right toward the knife-edge diffraction curve. The result is that the disparity 
between the measured and predicted insertion loss increases as x increases, as dis- 
cussed above for each individual case. In the case of the carpet-on-plywood surface, 
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Figure 8.20 Insertion loss above both the c&rpet-on-plywood and the ply- 
wood curved surfaces at x = 63, 120, 200, and 400 cm, when / = 10 kHz. 
This is a compilation of previously shown data. The discrete symbols 
( X ), (<>), (□), (+), and (v) represent the insertion losses measured above 
the laboratory scale ridge at the respective x locations. The solid and 
dashed lines represent the knife-edge diffraction formula (3.28), and G 
(3.26), at the same x locations. Y is the dimensionless height (3.29). 
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the background contribution is seen to increase the predicted insertion loss in the 
penumbra and shadow zone relative to that predicted by the knife-edge diffraction 
contribution alone. In this case, as z increases, the measured insertion loss clearly 
shifts to the left, but the predicted insertion loss also shifts to the left, yet still to- 
ward the knife-edge diffraction curve. However, the leftward shift of the predicted 
insertion loss does not keep pace with the leftward shift of the measured insertion 
loss. In fact, the background contribution is such that the predicted total insertion 
loss is always greater than the insertion loss predicted by the knife-edge diffraction 
term alone. A misleading result of this behavior is that the measured insertion 
loss at x = 200 cm happens to be shifted such that it is well fit by the knife-edge 
diffraction curve. At first glance, this result suggested that, at x = 200 cm, the 
carpet-on-plywood curved surface had already begun to effect the field in the same 
manner as a rigid knife-edge. However, the measured data at i = 400 cm did not 
also match the so-called knife-edge diffraction curve; in fact, the match was much 
worse. This result, along with the fact that the functional form, $, of the knife-edge 
diffraction term is identical to that of (3.31), suggested what was hinted and implied 
earlier: Y is not an appropriate dimensionless height when i is not small. 

A good candidate to replace Y is the dimensionless height 



which has the quality of being zero when y is zero. Further, the magnitude of ip is 
related to the angular distance of a point from the line of sight. This feature results 
in predictions of s imil ar insertion losses for points of similar angular distance from 
the line of sight. Such a result is in accord with the concepts of the Geometrical 
Theory of Diffraction. In that theory, diffracted rays are seen as having traveled 
around the curved surface from points of diffraction to points of departure from 
which the diffracted rays are shed tangentially and travel into the shadow zone. 
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Figure 8.21 Measured insertion losses above the curved plywood surface 
at x = 63, 120, 200, and 400 cm when / = 10 and 20 kHz, respectively. 
The discrete symbols ( x ),(□), (+), and (v) represent individual mea- 
surements made above the laboratory scale ridge. The abscissa is no 
longer Y but rather another dimensionless height, 0. As shown, 0 = 0 
is the line of sight. 
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Figure 8.22 Measured insertion losses above the curved ca.rpet-on-plywood 
surface at x = 63, 120, 200, and 400 cm when / = 10 and 20 kHz, 
respectively. The discrete symbols (O), (□), (+), and (v) represent 
individual measurements made above the laboratory scale ridge. The 
abscissa is no longer Y but rather another dimensionless height, rp. As 
shown, ^ = 0 is the line of sight. 
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Thus, for rays having been diffracted near the apex, the shadow zone rays will follow 
straight paths which would nearly intersect with the apex if they were extended 
back toward the apex. In this way, the GTD also predicts similar insertion losses 
for points in the shadow zone of similar angular distance from the line of sight, losses 
due to spreading not withstanding. Finally, the magnitude of $ is proportional to 
the square root of frequency. The result of this relationship is that any described 
width of the penumbra decreases with increasing frequency but the penumbra itself 
is nonetheless always centered about the line of sight. In Fig.( 8.21- 8.22), the 
experimental data at x = 63, 120, 200, and 400 cm is replotted with as the 
abscissa. It is clear from the figure that for both surfaces, t/> gave a reasonably 
good normalization of the data. It is worth noting that at small but non-zero z, 
Y is approximately equal to 0. Thus, at x = 63 cm, for instance, it is appropriate 
to refer to the solution G as being comprised of the superposition of a knife-edge 
diffracted field and the field of the background terms, (3.28). At large z, $(F) 
no longer approximates the knife-edge diffraction formula. 

8.3 Insertion Loss in the Penumbra and the Rigid Wedge Solution 

In chapter 3, and in the previous section, it was suggested that the insertion loss 
at large x behind a curved surface may be approximated by that behind a rigid 
knife-edge or thin wedge. This suggestion arose partly out of the form of (3.26- 
3.28) because the second term, was shown to be small when z is large. At such 
locations, G is approximately equal to $, which has the form of the rigid-knife edge 
solution. 

To investigate this suggestion, the insertion loss behind a wedge was measured 
at x = 63, 120, 200, and 400 cm from the vertex. The wedge used for these mea- 
surements was described in chapter 4. Figs. (8.23-8.26) show the data at the four 
downstream positions, and at / = 10 and 20 kHz. Again, the abscissae are the 
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Figure 8.23 Insertion loss behind a thin steel-tipped wedge at x = 63 cm 
when / = 10 and 20kHz, respectively. The discrete symbols (A), (x), 
and (O) represent the measured insertion loss behind the laboratory 
wedge, and the laboratory scale plywood and carpet-on-plywood ridges, 
respectively. The abscissa is the dimensionless height, rp, where rp = 0 
is the line of sight. 





INSERTION LOSS (dB) INSERTION LOSS (dB) 


161 



a Steel Wedge 



Figure 8.24 Insertion loss behind a thin steel-tipped wedge at x = 120 cm 
when / = 10 and 20kHz, respectively. The discrete symbols (A), (x), 
and (0) represent the measured insertion loss behind the laboratory 
wedge, and the laboratory scale plywood and carpet-on-plywood ridges, 
respectively. The abscissa is the dimensionless height, 0, where 0 = 0 
is the line of sight. 
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Figure 8.25 Insertion loss behind a thin steel-tipped wedge at x = 200 cm 
when / = 10 and 20kHz, respectively. The discrete symbols (A), (x), 
and (A) represent the measured insertion loss behind the laboratory 
wedge, and the laboratory scale plywood and carpet-on-plywood ridges, 
respectively. The abscissa is the dimensionless height, tp, where rp = 0 
is the line of sight. 
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Figure 8.26 Insertion loss behind a thin steel-tipped wedge at x = 400 cm 
when / = 10 and 20kHz, respectively. The discrete symbols (A), (x), 
and (0) represent the measured insertion loss behind the laboratory 
wedge, and the laboratory scale plywood and carpet-on-plywood ridges, 
respectively. The abscissa is the dimensionless height, 4>, where ip = 0 
is the line of sight. 
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dimensionless height, rp. Also shown in each figure is the predicted insertion loss 
behind a rigid wedge as well as the insertion loss measured behind the plywood 
and carpet-on-plywood curved surfaces. The predicted insertion loss was computed 

from 


-ikR 

A C 

f aS ~R- 


H(X)-^ x ’A„(X) 


(3.31) 


where 

X = ^ A<£ an d = ^ — (</>. — *) (3.32) 

which assumes that the listener is near the line of sight and sufficiently far from the 
vertex. In this formula, the dimensionless height, X, is similar to ip. Therefore, the 
actual computations were carried out with X replaced by rp. 

Computations based upon (3.31) compare fairly well with the experimental 
data obtained using the laboratory wedge. At x = 63, 200 and 400 cm, the difference 
is about ldB except at points deep in the shadow zone. At x = 120 cm, the 
predicted insertion loss behind a rigid wedge is approximately 2 dB greater than 
that measured behind the laboratory wedge. At all four x positions, the difference 
tended to decrease as tp increased. At points deep in the shadow zone, the larger 
deviation from the theory might be attributed to the limitation of the theory, namely 


that the listener be near the line of sight. 

Interestingly, at x = 120, 200 and 400 cm, the predicted thin wedge data 
matches the measured plywood curved surface data rather well in the lower bright 
zone; that is, when 1 < tp < 2. However, at x = 63 cm, the first oscillation in the 
plywood data begins at a smaller t p than does that in the predicted wedge solution, 
and so, the match is not quite as good in this case. These results are not surpris- 
ing since, as x gets large, the diffracted rays travel only a short distance along the 
top portion of the curved surface and are, therefore, influenced only by the surface 
properties in the immediate vicinity of the apex. In this way, the diffraction process 
is like that of a knife-edge. Since both surfaces are hard, one might expect similar 
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diffracted transmissions into the bright zone, and therefore, also similar insertion 
losses, when x is large. 

At all four x positions, and for any location in the penumbra or shadow zone, 
the insertion loss measured behind the wedge is greater than that measured behind 
the plywood surface and less than that measured behind the carpet-on-plywood 
surface. In addition, the difference in the measured insertion losses behind the 
wedge and behind the plywood curved surface data is slightly smaller at each x 
than at the next smallest x. The differences are approximately 3, 2, 1, and 1 dB at 
x = 63, 120, 200, and 400 cm, respectively. In the case of the carpet-on-plywood 
surface, the difference in the measured insertion losses is in the range of 4-6, 3-4, 
2-3, and 2-3 dB at x = 63, 120, 200, and 400 cm, respectively. 

In sum, the experimental data indicated that a hard curved surface, such as 
the plywood ridge, scatters sound to locations sufficiently far from the apex in a 
manner s imil ar to that of a rigid wedge. The data also indicated that a soft surface 
continues to give insertion losses greater than those of a rigid wedge even when x 
is not small. Finally, the data also suggests that the insertion loss in the penumbra 
and shadow zone of a curved surface whose finite impedance is softer than the 
plywood but harder than the carpet-on-plywood would be roughly equivalent to 
that of a thin rigid wedge, independent of x. 


8.4 Insertion Loss in the Shadow Zone and the Creeping Wave Series 
In chapter 3, it was shown that a residue series expression for G of the form 

» 2*' "if: , Wl{a ;, < 3 - 40 > 

u -<n «>?(<*») 

exists when rj is sufficiently small. It was suggested that reasonable convergence 
of the series can be expected for £ > r) 1 ^ 2 . In terms of the dimensionless height, 
Y, this convergence criterion corresponds to listener locations of Y <0; that is, 
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to listener locations in the shadow zone or penumbra. Each term of the series was 
identified with a creeping wave such that the field in the shadow zone and penumbra 
was seen as a superposition of an infinite number of creeping waves or of the rays 
which they shed. 

In the shadow zone, the series may require only a few terms to converge. How- 
ever, at listener locations near the line of sight, the number of terms required for 
convergence may be large. For example, Fig. (8.27) shows the number of terms 
required for convergence at each of 35 points spanning —0.58 < Y < 3.35 when 
x = 63 cm and / = 10 kHz for the plywood surface. The series was deemed to have 
converged when two consecutive terms in the series and their sum did not alter 
the total sum by 1 part in 10 6 . Nearly identical numbers were found at the same 
frequency and x position for the carpet-on-plywood surface. 

The insertion loss predicted by the creeping wave series when x = 63 cm and 
/ = 10 and 20kHz is shown in Fig.(8.28). As expected, the match with exper- 
imental data when Y < 0 was excellent, just as it was when G was computed 
as the sum of knife-edge diffraction plus background terms. However, the match 
continued to be excellent well into the bright zone. A plot of the creeping wave 
series solution compared to the knife-edge diffraction plus background terms solu- 
tion revealed that the two solutions were practically identical until Y = 2 and only 
diverged slowly thereafter. However, at some point, call it Y * , the creeping wave 
solution failed completely and abruptly. An example of this behavior is shown in 
the results at / = 20kHz. In those results, Y* was approximately 2.3 and 2.6 for 
the carpet-on-plywood and plywood surfaces, respectively. 

The reason for the abrupt failure at some Y* was unclear. However, it was 
significant that the creeping wave series converged at all when Y > 0. As stated 
in chapter 3, Pierce has suggested that the creeping wave series expression for G is 
valid when £ > rj 1/3 . 36 From the results of the figure, this statement is clearly 
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Figure 8.27 Number of terms required for convergence of the creeping 
wave series solution at various values of Y behind the curved plywood 
surface, when x = 63 cm and / = 10 kHz. Similar results were found 
when x = 120, 200, and 400 cm, and for the c&rpet-on-plywood surface. 
At more negative values of Y than those shown, the number of terms 
required for convergence continued to decline, although slowly, as might 
be extrapolated from the above curve. 
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Figure 8.28 Insertion loss behind the curved plywood and carpet-on- 
plywood surfaces at x — 63 cm when f — 10 and 20 kHz, respectively. 
The discrete symbols (x) and (0) represent the insertion loss measured 
behind the laboratory scale ridges. The solid lines represent the creeping 
wave series solution (3.40). Y is the dimensionless height (3.29). 
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true. However, the creeping wave series continued to yield approximately the same 
result as G even when £ < i7*^ a (i.e. Y > 0). The explanation for the apparent 
convergence of the series when Y > 0 is that the contribution from the contour 
integration at infinity in the upper half of the complex plane is non-zero but still 
quite small. If the contribution is small enough to be negligible, then the residue 
series would still approximate the integral solution, although not exactly. Until the 
computations were carried out, it was unclear just how far the { > ^ criterion 

could be exceeded. The results shown in the figure indicate that the criterion is 
quite conservative from a practical point of view. 

Not surprisingly, at x = 120 cm, the insertion loss predicted by the creeping 
wave series also compared rather well with the measured insertion loss (Fig. 8.29). 
Here, the good comparison is as expected since it has already been shown that 
the insertion loss predicted by the knife-edge diffraction plus background formula 
compares rather well with the measured insertion loss when x = 120 cm. As before, 
Y * ss 2.3 for the carpet-on-plywood case at / = 20 kHz. However, for the plywood 
case, y* was not found but was shown to be at least 3. 

At x = 200 and 400cm, Figs.( 8.30-8.31), the results are also similar to those 
found before. Again, the measured insertion loss was over predicted by the theory 
for both types of surface, particularly at x = 400 cm. At x = 200 cm and / = 
20kHz, the value of Y * for the carpet-on-plywood was approximately 2.4 while it 
was approximately 2.6 for the plywood. It would appear that the value of Y is a 
function of £, r}, and q, although no expression for Y‘ is known. 
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Figure 8.29 Insertion loss behind the plywood and the carpet-on-plywood 
curved surfaces at x = 120 cm when / = 10 and 20 kHz, respectively. 
The discrete symbols (x) and (<>) represent the insertion loss measured 
behind the laboratory scale ridges. The solid lines represent the creeping 
wave series solution (3.40). Y is the dimensionless height (3.29). 



171 



* Carpet on Plywood 



Figure 8.30 Insertion loss behind the curved plywood and carpet-on- 
plywood surfaces at z = 200 cm when / = 10 and 20 kHz, respectively. 
The discrete symbols (x) and (0) represent the insertion loss measured 
behind the laboratory scale ridges. The solid lines represent the creeping 
wave series solution (3.40). Y is the dimensionless height (3.29). 
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Figure 8.31 Insertion loss behind the curved plywood and carpet-on- 
plywood surfaces at x = 400 cm when / = 10 and 20 kHz, respectively. 
The discrete symbols (x) and (0) represent the insertion loss measured 
behind the laboratory scale ridges. The solid lines represent the creeping 
wave series solution (3.40). Y is the dimensionless height (3.29). 
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8.5 Summary 

In the first part of this chapter, the results of two geometric acoustic solutions, 
one for an incident plane wave, another for a point source, were compared to the 
experimental data. For listener locations in the bright zone, a comparison of the 
point source solution with the data proved to be excellent at x = 0 and 63 cm, 
particularly for the carpet-on-plywood surface. At the larger values of x = 120 
and 200cm, the comparison was not quite as good but was, in all cases, better than 
that produced by the plane wave solution. These results were satisfying because the 
actual experimental source, the spark source, was shown in chapter 5 to behave as 
a point source. For listener locations in the penumbra, neither geometric acoustic 
solution compared very well with the experimental data, although, for the plywood 
surface, the comparison was good through small but positive values of Y . In the 
case of the carpet-on-plywood surface, the comparison was poor for Y < 1. 

Next, results based upon the Fock-van-der-Pol-Bremmer function, G, were 
compared with the plane wave geometrical acoustics results and the experimental 
data. The comparison showed that G mimicked the geometric acoustics results in 
the bright zone, although not exactly. Disparities between the two solutions in 
the lower bright zone were attributed to diffraction effects implicit in G but not 
included in the geometric solution. These effects were particularly significant at 
the larger x locations. In the penumbra, G also accurately predicted the measured 
insertion loss at x — 0 and 63 cm. The same held true at 120 cm, although a 
small but consistent over prediction of the measured insertion loss was also found. 
At the larger z positions, x = 200 and 400 cm, G consistently over predicted the 
measured insertion loss in the penumbra and shadow zone. In the plywood surface 
case, the over prediction was approximately 2 dB and 8—10 dB at x = 200 and 
400 cm, respectively. In the carpet-on-plywood surface case, the over prediction 
was approximately 3— 5dB and 10 dB at the same respective x locations. It was 
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argued that, since the theory was derived by assuming (x,y) < R, the disagreement 
at the larger values of x was not surprising. An explanation for the consistent 
over prediction was found by examining G in terms of the knife-edge diffraction 
plus background formula. It was shown that the background component became 
progressively smaller as x increased, leaving G « $(Y). However, $(Y) was shown 
not to be a good solution at any x because of the nature of the dimensionless height, 
Y . In a cartesian plane, Y = 0 was shown to represent an upward sloping parabola 
with the origin as its vertex. The result was that the theory, minus the background 
component, predicted a line of 6 dB loss which lay higher and higher above the line 
of sight as x increased, independent of the type of surface. This result ran counter 
to the intuitive notion of a curved surface acting as a wedge when x is large. In the 
wedge solution, the 6dB loss line is the line of sight, independent of x. Further, it 
was shown that the experimental data at all four x locations was normalized rather 
well by another dimensionless height, t/>, found in the wedge solution. Since Y ~ xp 
when x is small, it was concluded that ip is a better suited dimensionless height 
than Y, and that the knife-edge diffraction plus background formula would yield 
better results everywhere if Y was replaced by ip. In sum, G was shown to provide 
excellent predictions at listener locations near to the curved surface and that its 
difficulties at large x could be corrected by substituting ip for Y . 

An attempt was made to check the notion that a curved surface acts as a 
wedge when the listener is at large x. The insertion loss behind a fairly rigid wedge 
was measured at z - 63, 120, 200, and 400 cm, just as was done for the curved 
surfaces. A comparison of the wedge data with that of the two types of curved 
surfaces revealed that the insertion loss of the wedge lay between that of the plywood 
and carpet-on- plywood surfaces. In fact, the measured insertion loss behind the 
wedge was only slightly greater than that behind the curved plywood surface. This 
difference decreased as x increased, so much so that the two sets of data practically 
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coincided in the penumbra at x = 400 cm. The difference between the wedge data 
and the carpet-on-plywood data also decreased as i increased, but not to the extent 
seen with the plywood data. Furthermore, the theoretical wedge solution over 
predicted the experimental data by approximately 1-2 dB. The conclusion drawn 
from these results was that a curved surface is fairly well approximated by a wedge 
when x is large, and particularly so when the curved surface is hard. Another 
conclusion was that for some curved surface, whose impedance is less than that 
of the plywood but greater than that of the carpet-on-plywood, the insertion loss 
in the penumbra and shadow zone would approximate that at the same locations 
behind a thin rigid wedge, independent of x. 

Finally, the creeping wave series expression for G was found to give equally good 
results as the knife-edge plus background formula, even when Y > 0. In general, 
the creeping wave series was shown to yield reasonable results for Y <2, and 
perhaps slightly higher. However, the computations also showed that the creeping 
wave series tends to fail rather abruptly at some height, designated Y . Preliminary 
indications are that Y' > 2 in most, if not all, cases. In many previous publications, 
it was suggested that the creeping wave series only converges for Y < 0. From a 
practical point of view, this has been shown to be not true. 

The creeping wave series results are important because the computation time 
involved is considerably less than that associated with the knife-edge plus back- 
ground expression. This fact is seen in that only two Fock functions need to be 
calculated for each term of the series while eight Fock functions are involved in the 
integrands of the background terms. In addition, the integration involved with the 
background terms required these eight Fock functions to be evaluated roughly 75 
times at each spatial point, (£, tj). There is a certain amount of time involved in 
calculating the roots or singularities needed by the series, but this time is short and 
the roots need only be calculated once for a given q. Thus, when the number of se- 
ries terms needed for convergence is small, the series solution is significantly faster. 
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Even at Y = 2, the creeping wave series required only about 100 terms. Thus, the 
creeping wave series offers no worse than a 3-to-l time savings over the expressions 
involving integration. When Y < 0, the time savings can be 20-to-l or greater. 
However, when the background terms are negligible, both solutions are inferior to 
a simple knife-edge solution, which is considerably faster than either expression of 
G. 

In s umma tion, the theory was quite good at predicting the measured insertion 
loss near caustics such as the surface and the not too distant line of sight. Further, 
the theory performed equally well for both the hard plywood surface and the fairly 
soft carpet-on-plywood surface. An important point is that the values of g used 
in the theory were deduced from measurements made over flat surfaces of similar 
impedances as the curved surfaces. In this way, no free parameters were used in the 
calculation of the predicted insertion losses. Thus, the MAE theory of diffraction 
over a curved surface of finite impedance appeared to be an excellent description 
of the field near such a surface, at least when kR is large. There is some prospect 
that the results may also be good for smaller kR on the order of unity since Keller 
previously found reasonable results using the GTD in this range. 
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CHAPTER IX 

CONCLUSIONS AND RECOMMENDATIONS 


Phenomena associated with long range propagation of sound over irregular topog- 
raphy motivated the research work which was described in this thesis. Specifically, 
the goal of the work was to analyze the diffraction effects which would occur near 
the tops of hills and ridges. From this particular goal, the research work evolved 
into a study of the diffraction of a high frequency plane wave due to its grazing of a 
two-dimensional curved surface of finite impedance. Laboratory scale models were 
constructed and measurements were made of the field on, above, and behind either of 
two curved surfaces possessing distinctly different impedances; that is, one was soft 
while the other was hard. The experimental technique consisted of simultaneously 
measuring the pressure at a reference point and at a field point due to a transient 
pulse generated by an electric spark. The pressure waveforms were digitized and 
processed. As described in chapter 4, the ratio of the discrete Fourier transforms 
of the two waveforms provided an estimate of the insertion loss between them. The 
results of the measurements were compared with the predictions of a theory which 
was derived by Pierce using the method of Matched Asymptotic Expansions (MAE). 
The predictions relied upon the experimental evaluation of the impedance of each 
surface at grazing angles of incidence. This evaluation was achieved by a fairly stan- 
dard technique involving empirical models of various generic types of surfaces. An 
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example was shown of the important role that the structural intricacies of a surface 
play in the determination of an appropriate model. The comparison between the 
mea s urements and predictions clearly indicated that the theory gives an excellent 
description of the field anywhere near a curved surface. Further, with a simple 
modification, the theory was also shown to give nearly as good of a description of 
the field surrounding a curved surface even at distances far behind the surface yet 
near the line of sight. 

The elementary yet fundamental nature of this study should make the results 
obtained therein useful in other studies or applications. For instance, the MAE 
theory, when used in conjunction with the versatile Geometrical Theory of Diffrac- 
tion (GTD), should prove to be a powerful tool for predicting the field in more 
complicated geometries. This expectation is based upon the fact that the GTD has 
provided good predictions of fields associated with canonical type problems, such 
as the propagation of a plane wave past a sphere. The weakness of the GTD is 
that it relies upon canonical solutions for the evaluation of its so-called diffraction 
coefficients. Thus, when the curved surfaces involved in a propagation problem are 
not of simple geometric shapes, such that a canonical solution is not available, the 
MAE theory could fill the breach since it requires only that the IocbJ curvature of 
a surface be smooth. An interesting extension of the current MAE theory would 
be to those cases where the local curvature is not smooth but rather is continually 
varying or discontinuous. In a related vein, problems involving surfaces of continu- 
ously varying impedance might be well handled by a variation on the current MAE 
theory as well. Further, previous work by Fock suggests that a generalization of the 
solution for a plane wave to the solution for a point source located near the surface 
could also be found by the methods used in the plane wave case. 

Another interesting and related problem is that of echoes occuring from the 
passage of a wave over curved surface. The work in this thesis dealt almost ex- 
clusively with the forward scattering problem. However, there may be cases where 
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diffracted echoes play an important role in the description of a field. An example 
of such a case would be when one curved surface sits directly in front of another . 
Then, echoes from the latter curved surface may make a significant contribution to 
the field in the shadow zone of its forward sitting partner, particularly when the 
forward sitting surface is soft. The same methods used in the derivation of the 
theory for the forward field should be applicable to the problem of solving for the 
backward field. 

In chapter 1, it was mentioned that a clear analogy exists between the dif- 
fraction of sound around a curved surface and the refraction of sound due to wind 
or temperature gradients above a flat surface. In fact, for the case of a linear 
sound speed profile, the forms of the respective solutions are identical. Therefore, 
an extension of the MAE theory to the case of a surface with continually varying 
curvature might also lead to the solution of more complicated refraction problems. 

Clearly, each of the above mentioned problems constitutes a separate aspect of 
the more general problem of long range sound propagation in the atmosphere above 
irregular terrain. Thus, the MAE theory and the method used to obtain it could 
be effective tools in the construction of a more comprehensive solution for the total 
field. 

In an even more general context, the MAE theory should be useful in suggesting 
solutions to other types diffraction problems such as the diffraction of pulses or 
bounded beams from curved surfaces. 
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APPENDIX A 


MATCHED ASYMPTOTIC EXPANSIONS 
AND THE THEORY OF 

SOUND PROPAGATION OVER A CURVED SURFACE 


This appendix contains many of the details involved in the derivation of the theory 
as presented in chapter 3. Most of what follows pertains to the derivations of (3.2), 
(3.21), and (3.26-3.28), and is taken directly from a collection of notes compiled by 
Pierce. Part of the derivation in section A.4 is due to Zhou. 


A.l Geometrical Acoustics Field near a Curved Surface 


Consider a horizontal ray which is incident upon a curved surface and which reflects 
according to the law of mirrors (Fig.3.2). Let R be the radius of curvature of the 
surface at the reflection point (z 0 ,!/o). The reflected ray path is given by 


x = x 0 + lcos(ir - 2&i) 
y = y 0 + £sin(7r - 20.) 


(A.l) 


where l is the distance along the reflected ray. 

If the incident field is a plane wave of unit amplitude, then the acoustic field 
complex amplitude in the geometric acoustics approximation is 


c -") 


p = + \A(Q) I c ikt 


(A.2) 
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where JR is the reflection coefficient, and A(l) is the so-called ray tube area associated 
with a small bundle of adjacent reflected rays. In two-dimensions, A[l) is the 
spreading of two adjacent reflected rays after propagation a distance i from the 
surface. If one considers the geometry associated with the reflection of two rays 
which strike the surface a distance As apart, then one can show that 


A(0) _ R cos Qj 
A{t) R cos 0i + 2 1 


(A.3) 


The reflection coefficient JR is that for plane wave reflection at a locally reacting 
surface such that 


5R = 


Z, cos &i — pc 


(A.4) 


Z, cos Oi -I- pc 

where Z, is the surface impedance, and pc is the characteristic impedance of the 
surrounding fluid medium. When these expressions are substituted into (A.2), a 
new expression for the acoustic field complex amplitude is 


* ikx . 

p = e + 


R cos 0i 

1/2 

Z, cos 0i — pc 

Rcos0 { + 21 


Z, cos 0i + pc 


p ik{* o + l) 


(A.5) 


The preceding analysis is valid for any surface in the geometrical acoustics 
limit. Of particular interest here is the form of (A.5) in the vicinity of the top of a 
curved surface. As shown in Fig. ( A.l), near the top one can write 

i 0 = — RsinOo 

y 0 = -/E(l - costfo) (A.6) 

$i = jt /2 - 0 O 

where 0 O is small and x, y, and l are all much less than R. Upon substitution of 
(A.6) into (A.l), one finds to “leading order”, 

x « —R0 o + t 


y « -R0% /2 + 2 10 O 


(A.7) 
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or, after elimination of £ 


-(ROoY +2{R9 o )x-Ry = 0 
& 


(A.8) 


This quadratic equation in R9 0 can be solved such that as x — * 0” and y — » 0 + , 
one has R9 0 > — i as dictated by the geometry. The solution is 

11/3 


nA -2x + [4x 3 + 6Ry] 

R0 ° = i 


(A.9) 


where the result of the square root is considered positive. Then (A. 7) yields 

, 1/2 


x + [4x 3 + 6Ry] 
l= - 


(A.10) 


(A.ll) 


If one adds £ to both sides of (A.l), then one has 

£ + x 0 = x + £(1 - cos 2 9 0 ) 

« x + 2191 

The higher order term must be kept because kR is assumed large, and so kl9* may 
not be negligible. Then, upon substitution of (A.9) and (A. 10) into (A.ll), one 
finds 

r- T 2 

(A.12) 



1 

1 

'2 

£ + x 0 « x + 2 

3 X + Q 

~R? 

[i 1 - Q \ 


where 


Q = 


4 2 2 _ 
-x 2 + -Ry 
9 3 


• 11/2 


(A. 13) 


With further algebra, one can derive 


l + x,=z+- 


.Ix’-f Rxy + Q > 


(A.14) 


Other factors in (A. 5) can be approximated in a similar fashion. In particular, 


R cos 9 t 


R9 a 


R cos 9i + 2 1 R9 0 + 21 


fss 


-i* + Q 


3 Q 


(A.15) 
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and 

Z, cos 0, - pc ~ 3 X + Q - / A 16 \ 

Z t cos di + pc — j x + Q + jrR 

Thus, in the limit of an incident ray reflecting as a plane wave near the top of a 
curved surface, the net result for the acoustic complex pressure-amplitude is 

i 1/1 r-j x + q - 


p si e 


■ ikx 


1 + 


where 


and again 




zh±Q 

3 Q 

2k 


-}x + Q + 


z R \ \ 


R 2 


—x 3 + -Rxy — Q 9 
27 3 


Q = 


4 2 2 

9 X + t Ry 


-. 1/3 


(A.17) 


(A. 18) 


(A. 19) 


which are identical to (3.2-3.4). As a check of (A.17-A.19), consider x < 0 and let 
2Ry = —x 2 such that (x,y) is on the surface. In this case, since y < 0, Q = — x/3, 
and ip = 0, so 

i — ’ — R 1 

Jkx ) i i ^ * \ 


p « C 




(A. 20) 


which is just 2e <fca in the limit of Z, -*■ oo for finite x. This result conforms to the 
expectation of pressure doubling at a rigid surface. 

Next, one can introduce scaling parameters L x and L y such that pe~ tkx of 
(A.17) is independent of k and R (although not of their product). Let 


x = xL x y = yL y 


Q = QL X 


If one substitutes for x and y in (A.18-A.19), and requires that ip be expressed 
solely in terms of dimensionless quantities, then the scaling parameters are found 
to be 

L x = R{kR)~ 1 ' 3 L y = J?(*fl)- a/3 (A.21) 
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With this terminology, another form of (A. 17) is 

a ikx 

p « e 


f 

1 + 

V. 


1/2 


•-A 

3 Q 


-\z + Q+?:(kRyi’ 

) 


where 


rp = xp = 2 

and the local shape of the surface is given by 

x 3 

y = —r 


+ \a - <?* 


(A.22) 


(A. 23) 


(A.24) 


A.2 Approximate Wave Equation for Diffraction by a Curved Surface 


To derive a solution valid on and near the surface, one attempts to find an approx- 
imate form of the Helmholtz (reduced wave) equation which can be solved readily. 
To begin, as was done in deriving the geometric acoustics solution, the local shape of 
the surface near the apex is approximated by the parabola y = — x 3 /2 R. Consider 
the Helmholtz equation 


f^P 

dx 2 


+ 


d h 

dy 2 


+ k 7 p = 


0 


when x < R and y < R and with the boundary condition 


(A.25) 


Vp • n H — —p = 0 at the surface, (A.26) 

Z, 

where n is the outward normal at a point on the surface. 

To facilitate the solution, one adopts a parabolic coordinate system which is 
a natural coordinate system for the given geometry. The relationship between the 
cartesian coordinates and the new parabolic coordinates is chosen to be 


x = u [1 + (v/f?)] 
y = v [1 + (v/2 R)] - tt 3 /2f2 


(A.27) 
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such that 

dx 
dy 

With this choice of coordinates, the diffracting surface corresponds to v = 0. The 
(u,v) coordinate system is orthogonal in that 


1 + v/R 

u/R 

du 

-u/R 

1 + v/R 

dv 


(A. 28) 


Vu • Vt; = 0 


(A.29) 


The determinant of the 2x2 matrix in (A.28) is 

D(u, v) = (1 + v/R) 2 + (u/R) 2 (A. 30) 

The basic rules of orthogonal curvilinear coordinate transformations can be used to 
express the Helmholtz equation (A. 25) in the parabolic coordinates as 

Likewise, the boundary condition (A.26) cam be reexpressed as 

+ o) p = o at t> = 0. (A.32) 

ov Z g 

Since the incident wave is planar, let 

p = e iku F(u,v) (A. 33) 


such that F satisfies 


n .,dF &F 9*F l2 

2lk JZ + + as + * 


/ 2u v 2 u 2 \ 

\R + R? + R?) 


F = 0 


(A. 34) 


with the same boundary condition (A.32). The second and fifth terms in (A. 34) 
are expected to be much smaller than the other terms when |u| C R and v C R. 
If the smaller terms are omitted, and if the (u,v) coordinates are scaled as in the 
geometrical acoustics solution such that 


v = L v v 


u = L x u 


(A. 35) 
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then approximate versions of the Helmholtz equation and its boundary condition 
are found to be 

2 *H + fr£ + (25 + « 2 ) F = 0 (A. 36) 

and 

|£ + !^p = 0 at v = 0 (A. 37) 

ov Z, 

where it is assumed that (fci2) 1/3 1. 

The next step was to derive an outer boundary condition on (A. 36) from the 
requirement that the solution match the geometrical acoustics solution at points 
well above the curved surface. Note that 

kx = ku + kuv/R 

= ku+ (kL,L y /R) uv (A. 38) 

= ku + uv 


and 


and 


Q 2 = \x 2 + | Ry = iu J (l + v/R) 2 + ? [Rv(l + v/2R) - u 2 / 2 ] 


9 

4 2 1 , 

« TtX - ~U 

9 3 


« Ll 


l„ a 
-U + -v 
9 3 


2 2 r-i 

2. 


8 2 2 8 2 1 2 2 
—i + -Ry w — u - -u + -Rv 
27 3 27 3 3 


« X 3 


2. 1 .2 

-V u 

3 27 


(A.39) 


(A.40) 


Thus, 


xp « 2 


1-2 2 . 
-U 2 + -V 

9 3 


- 11/2 


2 __ 1 
i" - 27“ - « 


(A.41) 

(A.42) 



Hence, the outer boundary condition was posed as a requirement that 




1 + 


„ - \ 1/2 


ZQ 


-!* + <?- ft ( fc *) 1/3 

-i u + Q+£{kRy/> 




(A.43) 


where 


Q = 


1-2 

5“ + 


H 


1/2 


(A.44) 


j, = tuv- — u 3 - 2 Q 3 (A.45) 

3 27 

This limit is expected to hold when u < 0 and v 1. 

If (A.43) is reasonable, then uv and uv — tp should be solutions of the eikonal 
equation 1 

n / dPhase\ ( dPhase\ 7 , 2 , A 

+ = 2v + u (A46) 

The phase uv is obviously a solution to (A.46). To confirm that the other phase 
uv — tp satisfies (A.46), note that 


and 


dtp 

du 



(A.47) 


dtp 4 . - 

-JL = - u - 2Q 
dv 3 


(A.48) 


Then, for Phase = uv — tp, one finds that 

I 2 


v — 


dtp 

du 


+ 


u — 


dtp 

dv 


2, 4_ 3 4 

-r+r + 5 tt<J 

-| { ’ + I“ a + 4 « J 

3 9 


1 

’ 1- A 

+ 

— 3 U ^ 


(A.49) 
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Simplified forms of (A. 36— A.37) result when the following notation is intro- 
duced. Let 


and 


.fkR\ 1/3 pc 

9 A 2 ) z. 

(A.50) 

t = — = - 5 - 

5 Re 2 1 / 3 

(A.51) 

2v 25 

V ~ Re 2 ~ 2 3 / 3 

(A. 52) 

-(A)'" 

(A. 53) 

G = e -«’/ 3 iT 

(A.54) 


such that G satisfies 



d 2 G 



+ i)G = 0 


with the boundary conditions 


— 1- qG — 0 at rj = 0 

dr\ 


and 


G 


where 




\ \ 3(?+3,)‘/» ) 

( -2t+(e+w 7 +x<i \ \ 

V-2£+(f J +Sl y/’-iii) I 


0 = - j£’> + {£ s } 


(A. 55) 


(A. 56) 


(A.57) 


(A. 58) 


Again, the above asymptotic limit is expected to hold when |£| is small and rj » 1. 
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In order to solve the boundary value problem (A.55— A.58), one can introduce 
a function T(»7, a) and assume a solution for G in terms of a Fourier type transform 
of T. That is, write 

&{£>*)) = f e’ 0< T {r),a)da (A.59) 

Jc 

where the contour C is independent of £ and tj. The introduction of this transform 
into (A.55) separates the variables and results in the requirement that 

T 

— + (r7 -a)T = 0 (A.60) 


which is the Airy differential equation. Two linearly independent solutions to the 
Airy differential equation are 


v(a - tj) = 7r 1/3 Ai(a - >7) 
tt;, (a — rj) = e'-^Tr^Ai (e‘ 2 ' /3 (a - »?)) 


(A.61) 


where v[z) and Wi ( z ) are so-called Fock functions and Ai(z) is the Airy function of 
complex argument. The Airy function is an entire function and can be expressed 
as a contour integral of the form 


Ai(a - tj) = J- [ «*!••/*+<— »M ds (A.62) 

2n Jc Ai 

where the integration contour C Ai can be any path originating at infinity in the 
sector 2?r/3 < arg(s) < n and ending at infinity in the sector 0 < arg(s) < tt/3. 

The boundary condition (A.56) is easily satisfied if one substitutes T(fj,a) = 
v(a — tj) + c(a)u>j (a — tj) and solves for c(a). The solution is found to be 


T = ff(«) 


, \ tF'(a) -g*(a) , x 

v(a-rj) 77-r 7 ~\ w i ( a ~ r? ) 

w[ (a) - qw x (a) 


(A.63) 


such that 


G{£,rj,q) = f e iai g{a) 
Jc 


v(ot-rj) - 


v'(a) — gt>(a) 
w[ (a) - qw x (a) 


w x (a - tj) 


da (A.64) 



190 


where it remains to fix the contour C and determine the function 5 ( 0 ) . One takes 
a clue from (A.57) and tries setting 


= f g* a< y(a)v(a — 17 ) da 

Jc 


(A.65) 


when Z is large and negative. One can substitute (A.61— A. 62) for v(o: 77 ) and 

obtain 

c -«*/3 c «n = _L [ e <.’/3 «-<•»* f e^+Og^dads (A.66) 

2 >A Jc Ai Jc 


Since 


f e ia{,+t) da = 2w6(s + Z) 

J — 00 

one is tentatively led to the choice 


(A.67) 



(A.68) 


with C traversing the real axis. Thus, a tentative solution is 


= ^ /__ [»(“ - ») - J' (a) - qwj^a) Wl (a — ^ ia < A ' 69 > 


What remains is to show that (A.69) has the asymptotic limit specified by 
(A.57-A.58). In order to show that such a limit is obtained, the asymptotic limits 
of the integrand of (A.69) can be analyzed. Let arg(a) = <p. Then one can write 


Ai(a) 


1 [<.-*«*'* /3 + 1 2?r/3 < <f> < 4tt/3 (A.70) 

27r 1 / 2 a 1 / 4 L J 


,- 3 a* /J /3 




-2tt/3 < <t> < 2tt/3 (A.71) 
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In fact, Ai(a) changes its behavior as one crosses rays of constant phase situated at 
arg(a) = titt / 3 in the complex a plane, where n is any integer. In particular, when 
a= |a|e*' 3 '/ 3 , 


Ai(|a|e <3 ' /3 ) 

When a = |a|e‘ 4 ’/ 3 , 

Ai(|a|e’‘ 4 ' /3 ) 

When a = \a\e iu /3 , 

Aidale’^ 3 ) 

When a = \a\e~ iK /\ 

Ai(|a|e-^/ 3 ) 




► 

2 ^i/ 3 | a |x/ 4 


e**/ 6 


2^/5 | a |l/ 4 




27T 1 / 3 lal 1 / 4 


e .>/is 


27T 1 / 3 |a| 1 / 4 


€ 3|a|*/*/3 




e -.'2|a| »'*/* 


e » 2 l a l* / */ 3 


Finally, when a = |a|e*" , 

Ai(|a|e ,,r ) 


cos(2|a| 3 / 3 /3-*r/4) 

jt 1 / 3 !^ 1 / 4 


(A.72) 

(A.73) 

(A.74) 

(A.75) 

(A.76) 


Thus, the asymptotic nature of the complex Airy function Ai(a) continuously 
changes from a decaying exponential on the positive real axis to an oscillator on 
the line <f> = tt/ 3, to a growing exponential on the line <f> = 27r/3, and back to an 
oscillator on the negative real axis. The behavior in the lower half-plane is a mirror 
image of the upper half-plane bacause Ai(a) = Ai(a). 

Consider the second term in the integrand of (A.69). The saddle point of this 
term is essentially the same as that of 


(A.77) 

(ot) 

and is most likely on the negative real axis. If one writes /? 2 = — a, then the saddle 
point is found from 

|f)-o <*•«> 



192 


where the multiplicity of signs is due to the fact that v(a) has two terms in its 
asymptotic limit for a real and negative. After the differentiation in (A. 78) is 
carried out, the equation of the saddle point becomes 

-2PZ =F 2/? 3 - 2/? 3 + (/? 3 + 17) 1/3 2/? = 0 (A.79) 


Two distinct forms for /? can be found depending upon which sign is chosen. These 
forms are 

j3 = -!e±i [£ 1 + 3 »] 1/3 or (A.80) 

The saddle point at /? = (£ 3 — rj) 1 / 2 exists only if £ > 0 and £ 3 > rj. However, 
these conditions correspond to points in the shadow zone. Thus, this saddle point 
is ignored. The other two possible saddle points can be checked to see if either is a 
proper solution of 

(/?= + T7) 1/3 = e-+-2>3 (A.81) 


which can be rewritten as 


-11/2 


= - j{ ± j(£ j +3|))‘ /5 


(A.82) 


When the upper signs are in place, the right side of the equation is always positive; 
however, when the lower signs are in place, the right hand side is always negative, 
which is a contradiction. Hence, the proper saddle point (denoted by subscript Jp ) 
is 


Ar = (— “.r) l '’ , = -f«+f(e’+3>l) ,/ ’ 

(A.83) 

(/?’+,)■/> = -i( +§(£»+ 3, )■'’ 

(A.84) 
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Further, the other terms in the argument of (A.78) can be expressed in terms of 
0, p with the result that 


-/3 s e-j0 1 + 5(/5 s + ») s/ ’ 


- je’ + (n (A ' 85) 

+ ^ («» + («’ +3,)^) 


which compares correctly with the exponent in (A.57-A.58). 

Near the saddle point, an approximation of the left side of (A. 85) can be found 
from the first two non-trivial terms of its Taylor series expansion. Let 


m = 


-0 , (-\F + l(P , +’i) a ' 2 


(A.86) 


Then 


m .... «/(/?.,) + row 






+ 






+ 


e 3 +4»7 


(A.87) 


(0 ~ A,)’ 


Thus, if one evaluates (A.69) by the saddle point method, the result is 



where 

0 = -^£>? + ^{e 3 + (e 2 +3»/) 3/3 } 


(A.88) 


(A. 89) 
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and 


^ FTTv 


Further analysis yields 

r dp « ~^ e ~ iW, *^lw 

J — oo 

Then, the second term in (A.88) can be rewritten as 

i (i±h\ ( jHP+M V 

</F+3 V \/» - I (/»’ + l) 1 " [*> - C(f +3-0"’] ) 


ab 


Finally, it can be shown that 


gMg ± 4 l) = -Hf + i(f» + = A, 

(/? 2 + *?) 1/3 [ 2>7 - £(£ 2 + 317 ) 1 / 3 ] 3 3 


(A.90) 


(A.91) 


(A.92) 


(A.93) 


Thus, upon comparison with (A.57-A.58), it becomes clear that (A.69) has the 
desired asymptotic solution and is the appropriate solution of the boundary value 
problem (A.55-A.58). 


A.3 A Study of a Definite Integral occurring in the Theory of High Fre- 
quency Diffraction by Curved Surfaces — Field on the Surface 


In the previous section, the integral solution 


■ ’ »■ - •> - 7? £ [«- - ’> - (a - ">] e ‘“ ! *• (A - 94> 


was identified, where 


cPv(a) 

da 7 

<Pw x (a) 
da 2 


—a v(a) = 0 
—otw x (a) = 0 


v(a — tj) = tt 1 / 2 Ai(a — rj) 
w x (a ~v) = e' ,r/e 2ir 1/2 Ai ^e* 2,r/3 (a — rj)^j 


(A.95) 


(A.96) 
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and 


Ai(a) 


27r 1 / 3 a 1 / 4 L 

2n i/3 a 1 / 4 




+ te 


iS*'"'*] 


2x/Z <4> < 47T/3 
- 2tt/3 < <j> < 2?r/3 


In the limit of the listener being on the surface rj = 0, one has 

G(f,0, ? ) = 1 r <(«)«(«) da 

y/ir J- oo «»i («) - 9^1 («) 

However, it follows from (A.95) that 

Wx v" — vw" = 0 
(ttf! v' — vw[) = 0 

Wx v' — w'x v = constant 

But as |a| — * oo for <j> = 0, one has from (A.96) and (A.97) 


2a 1 /* 


Wx 


,-2«* /3 /3 


,2a> /3 /3 


a 


1/4 


u>. 


—a lf2 v 


a'^Wx 


and so, 


Hence, when tj = 0 


Wx 


v' - w[v = —2a 112 — 7— = -1 
1 2 a 1 / 2 


1 /“ e* a 4 


In a 1945 article entitled “The Distribution of Currents Induced by 
Wave on the Surface of a Conductor” 2 , V.A. Fock defines 


?(e) - j?/ ri t»V) 


(A.97) 

(A.98) 

(A.99) 
(A. 100) 

(A.101) 

(A.102) 

(A.103) 
a Plane 


(A. 104) 
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where Tx is a contour which runs from ooe l2w/ 3 —►()—> oo, and 

tu(a) = - 4 = f e°‘e _, * /3 dt (A.105) 

V 7r J r, 

where T 3 is the mirror image of Tx in the complex plane. Now let 


s = e i6w/e t 


(A.106) 


so that 

at = iase' 2 * /3 and — i 3 /3 = is 3 / 3 

Then 

w(a) = -^e''/ 6 f t u ' /3 ds 

yj'n J Tt 

where r 3 is the contour e ,5,t / 6 oo — ► 0 — ► e* w / 8 oo. Thus, 

to (a) = 4=«‘ >/6 2?rAi(ae <3,r/3 ) 

V* 

= (a) 


(A. 107) 


(A. 108) 


(A.109) 


Further, it is well known that the zeros of w\ (a) are on the ray where arg(a) = x/3, 
so the contour r x can be deformed to coincide with the real axis. Therefore, 


G(f, 0,0) = ff (£) 


(A. 110) 


where g(£) is the function derived and tabulated by Fock. 

It is possible to deform the integration contour (— oo — *■ oo) in (A.95) to Fock’s 
contour r x , which will now be referred to simply as C. To prove that this deforma- 
tion is valid, it is helpful to introduce the connection formulae for the Airy function 
and its derivative 


Ai(t) = e in/3 A\ {tr i2 ”' 3 ) + e~ iw/3 A\ (te i2 * /3 ) 
Ai'(t) = e~ iw/3 Ai'(te* <3,r/3 ) + e iw/3 Ai' {te i2 " /3 ) 


(A.lll) 
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and define a third Fock function w 2 as a linear combination of u/j and v such that 
u> 2 (t) = uii (0 — 2iv{t) 

= 2n 1 ' a e i * / *Ai{tc i:,n '*) - 2*V 1/a Ai(t) (A.112) 

= 2v 1/7 e- iw/6 Ai(te- i3w/3 ) 

which may be compared with 

w, (t) = 2* 1 'V'' e Ai {te i3 *' 3 ) (A.113) 


From (A.61) and (A.70-A.71), it can be shown for \t\ large that 


u; 2 (f) 


1 if*/ 


t */ 4 


g 2 1* / */3 


S—l 0 < ^ < 4tt/3 

l *| l/4 


(A. 114) 


where <f> = arg(f). Similarly, it can be shown that 


Wl (0 — 2tt/3 < 4 < 2tt (A. 115) 


l *| 1/4 


Hence, 


tu 2 (|f|c^)| — » 0 as |i| — ► 00 and ir/Z<<f><v (A.116) 


Now, write 

v(f) = Yi W ” w * (01 (A. 117) 

and substitute for v(a) wherever it occurs in (A.95). The portion of the integrand 
associated with w x j2i drops out identically. Then, one is left with 





ta;(a) -gu> 2 (a) 
(<*) - gw 1 (a) 



(A.118) 
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In this form, it is clear that the contour can be deformed to C. Thus, it has been 
shown that 

G(£, n, 9) = —7= [ [u>2 (a — n) — -77-I — ^7-rWi (<* “ »?)] «** e da 

2 vW c L SH-gwxCa) J ( An9 ) 

*/[,(«_,)- - «"<«U (a _ ,)] da 

\/Wc L w i ( a ) - 9 w i ( a ) J 

Since [11/3 1 — * 0 as joe | — * 00 for <f> — 2x/3, and since |v| — ► 0 as |qe| ♦ 00 for <f> — 0, 

the first version of (A. 119) is more appropriate for the first leg of C, while the 

second version of (A. 119) is more appropriate for the second leg of C. Therefore, 

one can write 

G({, f) . i ( U (a - ,)]«•'■ < da 

K ’ 20F J legl l J (A120) 

+ -L/ f v(a-T))- 

vW le , 3 l «;;(«) -9tM«) J 

The integral over leg 1 in (A. 120) can be reexpressed by letting a = se‘ 2 ' /3 , 
and by noting the following identities 


w , 


1 (se ,2,r/3 ) = e ,,r/3 ti; a (s) 


and 


u;;(se‘ 2,r/3 ) = «-‘' /: X(s) 
w 3 (sc i7 *t 3 ) = 2e-‘ >/8 v(s) 


(A. 121) 


w' 2 {se i7w/3 ) = 2e- <6 ' /6 t;'(s) 

Substitution of these identities leads to another expression for G in the form 


(A.122) 


G(£ t T},q) = i J [*>(*-*?) 


w’ 7 (s) - qw 7 (s) 

«*(*)- H*) Wi 

w[ (s) - qw x (s) 1 


(s — »7)j t xii ds 
(s-» 7)]e’* f ds 


(A. 123) 


where fj = rje ~ i2w/ 3 , q = gc’ 2 * /3 , and s = sc i2 " /s . This form is preferable to (A.95) 
because the integrands become vanishingly small for values of s greater than s « 8. 

On the surface, one can simplify (A. 123) by multiplying the numerator and 
denominator of the first term of the integrand by the denominator of the second 
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term and combining the two terms. Then one can make use of the Wronskian 
relations W(v,iy 1 ) = W(v,u> 2 ) = 1, which were derived above to yield 


G(£,0, q) 



g-i.e/Sg-.ev's/a 

w 2 ( 5 ) “ e i2w / 3 qtv 2 (s) 
e‘*< 

rr T 7 ** 

< ( 3 ) - qw x (s) 


(A. 124) 


A.4 Fock’s Method for treating Transition from Light to Shadow in High 
Frequency Diffraction by a Curved Surface 

This section contains the results of an attempt to apply the analysis in Fock’s 
1948 paper “Fresnel Diffraction from Convex Bodies” to the approximate evaluation 
of 


c(e, i /_" [»(« - -1) - da < A12S > 


which is similar to an integral obtained by Fock except that Fock’s integration 
contour, C , is ooe‘ 2,r/3 — * 0 — ► 00 . 

It was shown in the previous section that an improved expression for G is 




(a) ~ quh (<*) 

<(<*) -<pi/ x (a) 
v'(a) — gv(a) 


w x 


(a 


da 


+ -V/ [.(«-«)- ° 

y/*Ju a 2 L wKa) -qwx (a) 1 


(A. 126) 


This expression for G can be rearranged by letting 


G = $ + 


(A. 127) 


where 



e*° ( tv 3 1 


(a-rj)da + -^= f e ia( v(a - rj)da 
V* Jteg2 


(A. 128) 
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and 


„ if i*t to, (a) — qtv 2 (a) . .. 

tf = — — 7= / « “tH ■rr\tv 1 {a-T))da 

Jiegi w i ( a ) ~ Q w i ( a ) 


(A. 129) 


1 f iat v'(a) - qv( a ) t \j 

7= / e 4 -77T --ru^a-rjjda. 

\/*J (ej , 3 ^(a) -gu;;(a) 

In the expression (A. 128) for $, leg 1 can be deformed to the contour (— oo — ► 0), 

and w 2 can be replaced by Wi — 2iv. Then, a new expression for $ is 


$ = 


(A.130) 


where 


and 


*•“ ijs/.. 1 

<$ 2 = — f c ,at v(a — rj)da = e"'* ^ 3 e ,<n 

V 7r » — oo 


(A. 131) 


(A.132) 


such that $ 2 is related to the incident wave. In order to evaluate , it is convenient 
to express the Fock function Wj (a — rj) in the integral form 

.iw/6 


w 1 (a - rj) = f 

v 71 " Jc Ai 


(A.133) 


where the contour C Ai = (ooe* 5 *’^ 6 — ► 0 — ► ooe*'/ 8 ) (previously r 3 ). Consequently, 
one has 

2 = te^_J j J da ds (A.134) 


or 




g»»/e r git* /3giti*e~ ,w> * 

> = / — — ds 

2?r Jr' £ — st~ x "! 3 

A i 

where the contour C' Ai is similar to C A , but passes above the pole at s = £e x *t 3 . 
The change of variable z = se XK ^ 3 allows the pole to be moved to the real axis, and 
the new expression for $i is 

-ti|i 


(A. 135) 


, If t~ % * ^e" 1 * 
1 ~ 2iV J Cr z - £ 


dz 


(A. 136) 
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where C r runs from ooe 1 *^ 2 , intersects the real axis to the right of the pole z = 
and then proceeds towards ooe’**/ 6 . Since the residue of the above integral at 
z = £ equals $2 j one can deform the contour C r to a similar contour Cj which, 
however, passes to the left of the pole, so that one can write 

$ = i, + i, = -L / ^ <f*. (A.137) 

The above integral can be evaluated by the saddle point technique. The saddle 
point occurs at z sp = ry 1 ^ 2 . Construct a new contour Co that is similar to both Ci 
or C ri except that it cuts the real axis at the saddle point. If r? 1/2 > £, the contour 
C 0 is equivalent to the contour C P , and $ reduces to $ 1 . If, however, rj 1 ^ 2 < £, the 
contour C 0 is equivalent to Cj . Near the saddle point z = z ip , 

+,*« In’ 12 (A. 138) 

o O 


so that 


±[ 

2:* J Ci 


and 


-L ( 

2l7r ./c, 






dz = $1 if t? l/2 > £ 


dz = $ if < £. 


(A. 139) 


By changing the variable of integration from z to u — i/ 1 ^ 4 (z — z, p )e*'^ 4 , and by 
deforming the integration contour to the real axis, one can actually perform the 
integration with the result that 

!^r *» , r = * if ,•/•<« (A. 140) 

L m (£ - »7 1/a ) - tte- ,,r / 4 Jj- 1 / 4 

that same integral being equal to if r? 1/2 > It follows that 




for X > 0 

+ 1 e (‘f+}'"’lyi I) (x) for X<0 
V2 


(A. 141) 
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where A D ( X ) is defined by 

L 


+ oo 


e— du 


(A.142) 


V2J-00 ( 5r / 2 )‘/ 3 X - e --/ 4 tt 

The parameter X acts as a dimensionless height in the normalized parabolic coor- 
dinate system (£,»?) and is deduced from (A.140-A.142) to be 




(A. 143) 


The properties and asymptotic limits of the diffraction integral Ad (X) can be found 
elsewhere. 3 

The function ^ can be reexpressed in the same manner as was (A. 116) in the 
previous section. Again, let a = se i2n ^ 3 for the leg 1 integral, and since 

Wi (se i2w/3 ) = c l " ,3 w 2 (s) 

(A. 144) 

w[{se i2 "' 3 ) = (s) 

and 

u; 2 (se ,2 ' /3 ) = 2c- in/ *v{s) 

(A.145) 


tv; 


^(se* 2 "' 3 ) - 2e~ iSK/fi v , (s) 


another expression for ^ similar to (A. 125) can be obtained. The final result is 

- -r f < A - 146 > 

•A-/0 «£(«) - qe' 2 w/ 3 w 3 [s) 

1 f°° i.t t>'(s)-gv(s) , \j 

— 7= / c—rri (\ w i (s-ri)ds 

V* Jo W 1 ( 5 ) - 9^1 ( 5 ) 

where H(Y) is the Heaviside step function which is zero for negative values of 
Y and unity for positive values of Y. The argument Y = - X , which represents 
a dimensionless height defined from (A.143), is positive in the bright zone and 
negative in the shadow zone. 
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APPENDIX B 

COMPUTATION OF AIRY FUNCTIONS 


This appendix contains a summary of an algorithm developed and coded by Pierce 
for the calculation of the Airy function and its derivatives for arbitrary complex 
argument. 1 This algorithm was applied in the computation of the solutions of 
chapter 3. 

In general, the Airy function of complex argument may be expressed as the 
difference of two functions of the form 


Ai(z) = Cj f{z) - Czg(z) 


(B.l) 


where 

Cl = 3“ 2/3 /T(2/3) = 0.355028053887817... 

c 2 = 3' 1/3 /r(l/3) = 0.258819403792807... 
and where the functions /(z) and g(z) are power series of the form 


/W = i+jV + 

g(z) = z + + 


2 

6 • 5 • 3 • 2 

z 7 

7 • 6 • 4 • 3 


+ 

+ 


z 9 + • • 

9 • 8 • 6 • 5 • 3 • 2 

z 10 + 

10 • 9 • 7 • 6 • 4 • 3 


(B.2) 


(B.3) 

(B.4) 


These series converge absolutely for all complex z. However, significant errors in 
the numerical results can arise for |z| > 3. 
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For very large \z\, an asymptotic series may be used to represent Ai(z). How- 
ever, the series does not converge absolutely because the terms in the series even- 
tually reach a minimum but then rise without limit thereafter. If one keeps all of 
the terms up to and including the minimum, then the error is on the order of the 
first excluded term. As a result of this behavior of the asymptotic series, another 
approach is needed for intermediate values of \z\. 

In order to derive this approach, the Airy function is expressed as a contour 
integral of the form 

Ai(z) = ~ / (B.5) 

Ztt Jc Ai 

where C A i is a contour in the complex s-plane going from ooe* 5 "/ 6 to the origin 
and from the origin to ooe’"^ 6 . With the help of the connection formulae of the 
Airy function, this contour can always be deformed to pass through a saddle point 
along a path of steepest ascent to the saddle point and then along a path of steepest 
descent away from the saddle point. Letting s = iz 1 ^ 2 + u and changing variables 
by letting f? — z 1/2 u 2 - (t'/3)u 3 , the Airy function may be rewritten, after further 
manipulation, as 


Ai(z) 


J_ c -(2/3)i’ /1 

27T 



21 

2 uz 1 / 2 — iu 2 


dt 


(B.6) 


Of the possible solutions of the cubic equation for t i(e//), only one fits the require- 
ments that u(0) = 0, that u(£) be continuous, and that the contour deformation 
onto the real l axis of the mapping of C A i be admissable. For small £, the solution 
is easily seen as u(t) = Ijz 11 * . Thus, it is appropriate to write 


u = 


Kl 

z i/< 


(B.7) 


Then, with further algebra, the Airy function can be expressed as 


Ai(z) 


e -(2/3)«* /a 

27T 1 / 2 Z 1 / 4 


1 




J — CO 




,3/4 



(B.8) 



205 


where 


and 



(B.9) 

(B.10) 


The integrand in the new expression for Ai(z) is no longer oscillatory and, thus, 
the integral has better convergence properties. An Hermite integration scheme was 
used to calculate the numerical value of the integral. A similar approach can be 
applied to solve for the derivative of the Airy function when |z| is moderate. 
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APPENDIX C 

GEOMETRIC ACOUSTICS SOLUTION 
FOR A POINT SOURCE 
RADIATING OVER A CURVED SURFACE 


The geometric acoustics solution derived here for a point source radiating over a 
curved surface of finite impedance follows the same approach as that used to derive 
the geometric acoustics solution for a plane incident wave. The two solutions are 
similar as they should be; in fact, the point source should yield the plane wave 
solution when the source position is moved to infinity. The point source solution is 
not presented in as neat and compact a form as was done for the plane wave solution 
in chapter 3 and appendix A. However, the goal here was to simply formulate a 
point source solution that contained the important approximations used to derive 
the plane wave solution. No attempt was made to present the solution in a form 
suitable for use as an asymptotic limit, as was done in the plane wave case. 

Consider two rays which are incident upon a curved surface and which reflect 
according to the law of mirrors (Fig. C.l). Again, let R be the radius of curvature 
of the surface at the reflection point {x 0 ,yo ) ■ <f > o is the angle between the first leg 
of the reflected ray and the line of sight. Let the position of the point source be 
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Figure C.l Geometric description of the acoustic field due to a point source 
in the vicinity of a curved surface. The reflection occurs according to 
the law of mirrors. Although the incident wave is spherical (or cylin- 
drical), the effects at the point of reflection are described by the plane 
wave reflection coefficient. However, the spreading of adjacent rays after 
reflection does include the effect due to the fact that the rays are not 
parallel when incident. 
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designated by the coordinates (x, ,0). Then, the second leg 
is given by 

X — Xq + /cOs(7T - 2 &i — 4> o) 


of the reflected ray path 


(c.i) 


y = y 0 + £sin(7r - 20, - <f> 0 ) 
where t is the length of the second leg. 

If the source strength is designated by S , then the acoustic field complex am- 
plitude in the geometric acoustics approximation is 

'*( 0)1 1/J 


£ 

A 

5 


.ikr. 


+ 


r o [A(l)\ 


&c ik r °c ikt 


(C.2) 


where 9? is the reflection coefficient, rj is the length of the direct path between the 
source at (x,,0) and the listener at (x,y), and A(l) is the so-called ray tube area 
associated with a small bundle of adjacent reflected rays. In two-dimensions, A(£) 
is the spreading of two adjacent reflected rays after propagation a distance l from 
the surface. If one considers the geometry associated with the reflection of two rays 
which strike the surface a distance As apart, then one can show that 


A(0) Rr% cos 0, 

A(t) ~ Rr l cos + 2r 0 £(i2 cos 0< + r 0 ) + 1? (R cos 0< + 2r 0 ) 


(C.3) 


where 0< is the angle of incidence normal to the surface. When r 0 — ► oo, this 
expression reduces to that given in Appendix A for an incident plane wave. The 
reflection coefficient 9i is that for plane wave reflection at a locally reacting surface 
such that 


a _ Z, cos 6j - pc 
Z, cos Oi + pc 

where Z. is the surface impedance, and pc is the characteristic impedance of the 
surrounding fluid medium. 


(C.4) 
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The preceding analysis is valid for any surface in the geometrical acoustics 
limit. Of particular interest here is the form of of the solution in the vicinity of the 
top of a curved surface. As in appendix A, near the top one can write 

x Q = — JZsinflo 

I/O = -i2(l - cos0 o ) ( c - 5 ) 

Of = ?r/2 — 9 0 — 4>o 

where $ 0 and <f> 0 are small and x, y, and l are all much less than R. Upon substi- 
tution of (C.5) into (C.0), one finds to “leading order”, 

x « —R6 o + 1 

(C.6) 

y a -ROll 2 + H2$ 0 +M 

or, after elimination of l 


^(iZ0 o ) 2 + {R<t > o + 2x)(i20 o ) + R<f> 0 x - Ry = 0 


(C.7) 


As before, the quadratic equation in R0 O can he solved so that as x — ♦ 0 and 
y — ¥ 0 + , one has R8 0 > — x, as required by the geometry. The solution is 


R6 o = 


■{R<f>o + 2x) + [(2x + R<f> o) 2 - 6iZ(^ 0 x - y)] 


1/2 


5 -W° + 2 *) +Wo) 


(C.8) 

(C.9) 


where the result of the square root is considered positive. When r 0 — ► oo, then 
♦ 0 and the above result again reduces to that found for the case of a plane 
wave. Substitution of (C.9) into (C.6) yields 


<= (x-^) +Wo) 


(C.10) 


An expression for ro can also be found in terms of R9q and R<j> o , but the algebra is 
messy. However, the expression for r 0 does reduce to Xo when <j>o — + 0, as it should. 
The other terms in (C.2) are also functions of R0 o and R<{> 0 . Again, they are not 
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written here, but when <t> 0 -* 0, they all reduce to their corresponding version in 
the plane wave solution. 

What remains is to express R<f > o in terms of ROq , or vice versa, so that the 
solution can be computed. A suitable expression is found from the geometry shown 


in the figure. The law of sines gives 


sin(<£+^ 0 ) _ sin 2 Bj _ sin(2^ 0 + 2<ft 0 ) 
l r, Ti 


(c.u) 


where <f> is the angle between the direct ray and the line of sight. When Bq and 4>o 
are small enough, then sin(20 o + 2<f>o) 2sin(0 o + ^o)> and 30 


sin(0 + <£ o ) 2sin(^ 0 +^ 0 ) 

t ~ r. 


(C.12) 


Further, when the listener is not too far above the curved surface, so that <f> is also 


small, then 


Further algebra yields 


<j> + 4>o 2(9 0 + <t>o) 

t ~ r. 


2xR0p +2iZ 3 gg - rM 
r x - 2(x — RB 0 ) 


(C.13) 


(C.14) 


This expression along with (C.8) can be solved iteratively. Experience indicates 
that a small initial guess for either of 6 0 or <f > o will quickly lead to a solution. When 
the height of the listener above the curved surface is not small compared to R, the 
above approximate formula for <f> 0 in terms of 0 O should be used with caution. In 
those cases, one can simply use the law of sines relation directly. Of course, the 
validity of the entire solution is questionable at too large listener heights, so the 
above approximate formula may suffice. 
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